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FOREWORD 


This  report  was  prepared  by  McDonnell  Aircraft  Company 
(MCAIR),  St.  Louis,  Missouri,  for  the  Structures  Research  Branch 
of  the  Naval  Air  Development  Center  ' (NADC ) ,  Warminster, 
Pennsylvania,  under  Contract  No.  N62269-79-C-0275.  Mr.  M.  S. 
Rosenfeld  was  the  Project  Engineer. 


The  Structural  Research  Department  of  McDonnell  Aircraft 
Company  had  responsibility  for  performance  of  this  program.  The 
program  manager  for  MCAIR  was  H.  D.  Dill.  Principal  author  of 
this  report  is  C.  R.  Saff.  W.  D.  Nowak,  MCAIR,  was  a  consultant 
for  development  of  the  landing  gear  load  spectrum.  K.  C.  Garland 
and  F.  J.  Coffey,  MCAIR  Metallurgical  Laboratories,  performed  the 
test  program. 


This  report  covers  work  accomplished  during  the  period  August 
1979  through  August  1980. 

This  report  was  released  by  the  author  in  October  1980  for 
publication. 


NADC-79095-60 


TABLE  OF  CONTENTS 


Section  Page 

I  INTRODUCTION  .  1 

II  MODEL  CALIBRATION  AND  EVALUATION  PROGRAM  .  3 

1.  TEST  PROGRAM  SUMMARY .  3 

2.  TEST  MATERIALS .  3 

3.  TEST  SPECIMENS .  5 

4.  TEST  PROCEDURES .  7 

5.  SUMMARY  OF  CONSTANT  AMPLITUDE  TEST  RESULTS.  ...  10 

6.  CONSTANT  AMPLITUDE  FATIGUE  STRESS  RATIO 

EVALUATIONS .  11 

7.  CONSTANT  AMPLITUDE  FATIGUE  FREQUENCY  EVALUATIONS.  16 

8.  CONSTANT  AMPLITUDE  FATIGUE  WAVE  SHAPE 

EVALUATIONS  .  .  .  . .  2  5 

9.  SPECTRUM  AND  DUPLICATE  TESTS .  45 

10.  SURFACE  FLAW  TESTS .  47 

III  ENVIRONMENT-LOAD  INTERACTION  MODEL  .  52 

1.  SUMMARY .  52 

2.  STRESS  RATIO  EFFECTS .  52 

3.  LINEAR  SUPERPOSITION  APPROACH  FOR  PREDICTION  OF 

ENVIRONMENTAL  ACCELERATION .  53 

4.  WILLENBORG  MODEL .  62 

5.  MODEL  CALIBRATION  .  66 

6.  SEMI -ELLIPTICAL  SURFACE  FLAW  ANALYSIS  .  68 

IV  LANDING  GEAR  STRESS  HISTORY .  71 

1.  F-4J  LANDING  GEAR  DESCRIPTION .  71 

2.  LOADS  SPECTRUM  SUMMARY.  '. .  74 

3.  FLIGHT-BY-FLIGHT  STRESS  HISTORY  DEVELOPMENT  ...  74 

4.  TIMES  FOR  STRESS  APPLICATION .  81 

5.  CYCLE-BY -CYCLE  STRESS  SPECTRA  .  85 

V  VERIFICATION  TEST  PROGRAM .  90 

1.  TEST  PROGRAM  SUMMARY .  90 

2.  TEST  SPECIMENS  AND  TEST  CONDITIONS .  90 

3.  PRECRACKING  PROCEDURES .  91 

4.  TEST  PROCEDURES,  CRACK  GROWTH  MONITORING, 

INSTRUMENTATION,  AND  ENVIRONMENTAL  CONTROL.  ...  92 

VI  COMPARISON  OF  SPECTRUM  LIFE  PREDICTIONS  AND  TEST 

RESULTS .  94 

1.  SUMMARY .  94 

2.  CENTER  CRACK  PANEL  SPECIMEN  RESULTS  .  94 

3.  SURFACE  FLAW  SPECIMEN  TEST  RESULTS .  99 

4.  CRACK  GROWTH  UNDER  NAVY  LANDING  GEAR  SPECTRUM  .  .  102 


v 


NA DC-79095-60 


TABLE  OF  CONTENTS  (Concluded) 

Section  Page 

VII  CONCLUSIONS .  108 

APPENDIX  A  CRACK  GROWTH  ANALYSIS  ROUTINE  INCLUDING 

ENVIRONMENT-LOAD  INTERACTION  EFFECTS  .  Ill 

APPENDIX  B  STRESS  INTENSITY  SOLUTIONS  FOR  ELLIPTIC 

SURFACE  FLAWS .  128 

REFERENCES .  139 


vi 


NADC-7 909 5-60 


LIST  OF  ILLUSTRATIONS 

FIGURE  PAGE 

1  Center  Cracked  Panel  Specimen  .  5 

2  Surface  Flawed  Panel  Specimen  .  6 

3  Containment  Method  for  Low  Humidity  Air  Tests  ...  8 

4  Container  Used  for  Alternate  Immersion  Tests.  ...  9 

5  Algorithm  Test  Program  Summary  Procedure .  10 

6  Effect  of  Stress  Ratio  on  Crack  Growth  in  300M 

Steel .  11 

7  Effect  of  Stress  Ratio  on  Crack  Growth  in  30 0M 

Steel .  12 

8  Effect  of  Stress  Ratio  on  Crack  Growth  Rate  in 

300M  Steel  (Plotted  versus  Kmax) .  13 

9  Effect  of  Stress  Ratio  on  Crack  Growth  in  HP 

9-4-0.30  Steel .  14 

10  Effect  of  Stress  Ratio  on  Crack  Growth  Rate  in  HP 

9-4-0.30  Steel .  15 

11  Effect  of  Stress  Ratio  on  Crack  Growth  Rate  in  HP 

9-4-0.30  Steel  (Plotted  versus  Kmax) .  16 

12  Effect  of  Sine  Wave  Loading  Frequency  on  Crack 

Growth  in  300M  Steel .  17 

13  Effect  on  Sine  Wave  Loading  at  10  cps  on  Crack 

Growth  Rate  in  300M  Steel  in  Sea  Water .  18 

14  Effect  of  Sine  Wave  Loading  at  1  cps  on  Crack 

Growth  Rate  in  300M  Steel  in  Sea  Water .  19 

15  Effect  of  Sine  Wave  Loading  at  0.1  cps  on  Crack 

Growth  Rate  in  300M  Steel  in  Sea  Water .  20 

16  Effect  of  Sine  Wave  Loading  Frequency  on  Crack 

Growth  in  HP  9-4-0.30  Steel .  21 

17  Effect  of  Sine  Wave  Loading  at  10  cps  on  Crack 

Growth  Rate  in  HP  9-4-0.30  Steel  in  Sea  Water  ...  22 

18  Effect  of  Sine  Wave  Loading  at  1  cps  on  Crack 

Growth  Rate  in  HP  9-4-0.30  in  Sea  Water .  23 


vii 


NADC-79095-60 


LIST  OF  ILLUSTRATIONS  (Continued) 

FIGURE  PAGE 


19  Effect  of  Sine  Wave  Loading  at  0.1  cps  on  Crack 

Growth  Rate  in  HP  9-4-0.30  Steel  in  Sea  Water  ...  24 

20  Trapezoidal  Wave  Shape .  25 

21  Effect  of  Trapezoidal  Wave  Loading  Frequency  on 

Crack  Growth  in  300M  Steel .  26 

22  Effect  of  Alternate  Immersion  Environment  on  Crack 
Growth  Rate  on  300M  Steel  -  10  cps  Trapezoidal 

Wave  Loading .  27 

23  Effect  of  Alternate  Immersion  Environment  on  Crack 

Growth  Rate  in  300M  Steel  -  1  cps  Trapezoidal  Wave 
Loading .  28 

24  Effect  of  Alternate  Immersion  Environment  on  Crack 

Growth  in  30 0M  Steel  -  0.1  cps  Trapezoidal  Wave 
Loading .  29 

25  Effect  of  Trapezoidal  Wave  Loading  Frequency  on 

Crack  Growth  on  HP  9-4-0.30  Steel .  30 

26  Effect  of  Trapezoidal  Wave  Loading  at  10  cps  on 

Crack  Growth  Rate  in  HP  9-4-0.30  Steel .  31 

27  Effect  of  Trapezoidal  Wave  Loading  at  1  cps  on 

Crack  Growth  Rate  in  HP  9-4-0.30  Steel .  32 

28  Crack  Growth  Rate  in  HP  9-4-0.30  Steel  in  Alternate 
Immersion  Environment  Using  0.1  cps  Trapezoidal 

Wave  Loading .  33 

29  Effect  of  Stress  Ratio  on  Crack  Growth  in  Alternate 

Immersion  Environment  -  300M  Steel .  35 

30  Crack  Growth  Rate  for  300M  Steel  in  Alternate 

Immersion  Environment,  R  =  0.5 .  36 

31  Crack  Growth  Rate  for  300M  Steel  in  Alternate 

Immersion  Environment,  R  =  -1.0 .  37 

32  Effect  of  Stress  Ratio  on  Crack  Growth  in  Alternate 

Immersion  Environment  -  HP  9-4-0.30  .  38 

33  Crack  Growth  Rate  in  HP  9-4-0.33  Steel  in  Alternate 

Immersion  Environment,  R  =  0.5 .  39 

34  Crack  Growth  Rate  in  HP  9-4-0.30  Steel  in  Alternate 

Immersion  Environment,  R  =  -1.0 . .  .  .  .  40 


viii 


NA DC-79095-60 


4 

LIST  OF  ILLUSTRATIONS  (Continued) 

FIGURE  PAGE 

35  Effect  of  Wave  Shape  and  Environment  on  Crack 

Growth  in  300M  Steel  at  10  cps .  41 

36  Effect  of  Environment  and  Wave  Shape  on  Crack 

Growth  in  HP  9-4-0.30  Steel  at  10  cps .  42 

37  Effect  of  Wave  Shape  on  Crack  Growth  in  300M 

Steel  in  Sea  Water  at  0.1  cps .  43 

38  Effect  of  Wave  Shape  on  Crack  Growth  in  HP 

9-4-0.30  Steel  in  Sea  Water  at  0.1  cps .  44 

39  Spectrum  Test  Results  for  FP  9-4-0.30  in  Dry  Air.  .  45 

40  Spectrum  Test  Results  for  300M  in  Dry  Air .  46 

41  Comparison  of  Surface  Flaw  Test  and  Analysis 

Results  in  300M  Steel .  48 

42  Comparison  of  Surface  Flaw  Test  and  Analysis 

Results  in  HP  9-4-0.30  Steel .  49 

43  Comparison  of  Surface  Flaw  and  Center  Crack 

Growth  Rates  for  300M  .  . .  50 

44  Comparison  of  Surface  Flaw  and  Center  Crack 

Growth  Rates  for  HP  9-4-0.30 .  51 

45  Crack  Growth  Behavior  of  300M  Steel  in  Alternate 

Immersion  in  Sea  Water  at  Low  Cycle  Rates .  55 

46  Comparison  of  Sustained  Load  Crack  Growth  Rates 
in  300M  Steel  for  Alternate  and  Continuous 

Immersion  Environments .  56 

47  Comparison  of  Results  of  Simple  Load-Environment 

Analysis  With  Crack  Growth  Behavior  of  300M  Steel  .  59 

48  Comparison  of  Average  da/dN  Curve  with  Crack 

Growth  Rate  Test  Data .  60 

49  Spectrum  Load  Definitions  Used  for  Analysis  ....  61 

50  Improved  Analysis  of  Compression  Effects  with 

Residual  Stress  Intensity  Model  .  65 

51  Effect  of  Overload  Ratio  on  Delay  in  HP  9-4-0.30 

Steel .  67 


ix 


k 


NA DC-79095-60 


LIST  OF  ILLUSTRATIONS  (Concluded) 

FIGURE  PAGE 

52  Comparison  of  Surface  Flaw  Stress  Intensity 

Solutions .  69 

53  F-4J  Main  Landing  Gear  Selected  as  Study  Base  ...  72 

54  Areas  of  F-4J  Main  Gear  Predicted  to  Sustain 

Greatest  Fatigue  Damage  .  73 

55  Loading  Diagrams  for  F-4J  Main  Gear  Fatigue 

Spectrum . . .  77 

56  Distribution  of  Landing  Stress  Cycles  in 

Stress  History .  79 

57  Comparison  of  FI ight-by-Fl ight  Spectrum  and 

Design  Spectrum  for  F-4J  Main  Landing  Gear .  80 

58  Stress-Time  Profiles  for  Each  Loading  Condition  .  .  86 

59  Real  Time  Stress  History  for  Longest  Flight  ....  88 

60  Real  Time  Stress  History  for  Shortest  Flight.  ...  89 

61  Development  of  Initial  Flaws .  92 

62  Comparison  of  Analysis  and  Verification  Test 

Lives  for  Center  Cracked  Panels .  96 

63  Verification  Test  Results  for  300M  Center 

Cracked  Panels .  97 

64  Verification  Test  Results  for  HP  9-4-. 30  Center 

Cracked  Panels .  98 

65  Comparison  of  Analysis  and  Verification  Test 

Lives  for  Surface  Flawed  Panels .  100 

66  Verification  Test  Results  for  300M  Surface 

Flawed  Panels  .  101 

67  Verification  Test  Results  for  HP  9-4-0.30 

Surface  Flawed  Panels  .  102 

68  Predicted  Growth  of  Small  Surface  Flaws  Under 

Navy  Landing  Gear  Spectrum .  103 

69  Comparison  of  Air  Force  and  Navy  Main  Landing 

Gear  Stress  Spectra .  104 

70  Comparison  of  Navy  Landing  Gear  Spectrum  Without 
Carrier  Landings,  and  Air  Force  Landing  Gear 

Spectrum .  105 

71  Comparison  of  Air  Force  and  Navy  Main  Landing  Gear 

Stress  Spectra .  106 

B-l  Surface  Flaw  Stress  Intensity  Solution  of  Depth  .  .  129 

B-2  Surface  Flaw  Stress  Intensity  Solution  at  Surface  .  130 


x 

1 


NADC-79095-60 


TABLE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


LIST  OF  TABLES 


PAGE 


Model  Calibration  and  Evaluation  Test  Program 

Summary  .  4 

Tensile  Test  Results .  4 

Summary  of  Constant  Amplitude  Lives  .  10 

Parameters  Used  for  Stress  Ratio  Corrections.  ...  53 

Comparison  of  Test  and  Analysis  Results  in  300M 

Steel .  57 

Comparison  of  Test  and  Analysis  Results  in  HP 

9-4-0.30  Steel .  58 

Willenborg  Model  Parameters  Used  for  Analysis  ...  66 

F-4J  Landing  Gear  Repeated  Load  Conditions .  75 

F-4J  Main  Landing  Gear  Design  Fatigue  Spectrum.  .  .  76 

Stress  Spectrum  for  Axle-Piston  Fork  Area .  7  8 

Flight  Stress  Sequence  for  Axle-Piston  Fork  Area.  .  79 

Comparison  of  Design  and  Flight-By-Flight  Stress 
Spectra  F-4J  Main  Landing  Gear  Axle-Piston  Fork 
Spectra  8,000  Takeoff  and  Landings .  81 

Summary  of  Stress-Time-History  for  Baseline 

Flight .  87 

Verification  Test  Program  .  90 

Summary  of  Spectrum  Test  Lines . . .  95 

Comparison  of  Test  Lives  Under  Air  Force  and  Navy 
Landing  Gear  Spectra .  107 


xi 


NADC-79095-60 


SECTION  I 
INTRODUCTION 


Aircraft  are  subjected  to  a  combination  of  environmental 
attack  and  varying  loads.  Structural  integrity  can  be  impaired  by 
surface  degradation  due  to  corrosive  action  or  when  crack  damage 
is  developed  or  aggravated  by  the  environment. 

Landing  gear  of  Navy  aircraft  are  subjected  to  an  extreme 
combination  of  aggressive  chemical  environment  and  high  load 
levels,  particularly  during  carrier  operations.  Landing  gear 
structure  usually  consists  of  a  series  of  nonredundant  elements. 
Weight  and  geometrical  constraints  require  that  landing  gear  be 
fabricated  from  high  strength  materials,  even  though  these 
materials  may  be  susceptible  to  environmentally  accelerated  crack 
growth. 

Great  care  is  taken  during  manufacturing  and  processing  to 
ensure  that  these  nonredundant  structural  elements  are  as  flaw- 
free  as  possible.  However,  even  with  good  quality  control,  minute 
cracks  can  be  developed,  for  example  during  grinding  and  plating 
of  nigh  strength  steel  parts.  Also,  as  with  any  aircraft  struc¬ 
tural  subsystem,  small  cracks  can  be  initiated  during  service  life 
from  latent  damage  sites  developed  by  fretting,  pitting,  intergran¬ 
ular  stress  corrosion  cracking  (for  aluminum  alloys)  and  fatigue. 

In  landing  gear  structural  components,  the  cracks  of  principal 
concern  initiate  on  the  surface  of  the  component.  Cracks,  whether 
initially  present  or  service  developed,  can  propagate  to  failure 
in  these  highly-loaded  nonredundant  structural  subsystems  and 
cause  system  failure. 

The  objective  of  this  program  was  to  systematically  investi¬ 
gate  crack  growth  behavior  under  chemical  environment  and  loading 
combinations  typical  of  Navy  aircraft.  Both  constant  amplitude 
and  spectrum  tests  were  performed.  Stress  histories  were  selected 
to  represent  usage  conditions  for  a  landing  gear  component  of  a 
multi-mission,  carrier-based  Navy  aircraft.  Experimental  results 
were  used  to  evaluate  the  ability  of  available  analysis  procedures 
to  predict  crack  growth  behavior  under  aggressive  environments. 

In  Task  I,  Model  Calibration  and  Evaluation  Test  Program, 
constant  amplitude  and  spectrum  crack  growth  data  were  obtained  to 
calibrate  and  evaluate  the  current  analysis  procedure.  Twenty- 
eight  specimens  were  tested  in  dry  air  and  alternate  immersion  in 
sea  water  environments  under  a  variety  of  loading  frequencies. 
Results  of  the  alternate  immersion  tests  were  used  to  revise  the 
crack  growth  analysis  methodology  to  include  analysis  of  the 
decaying  rate  of  crack  growth  which  occurs  during  drying  in  the 
alternate  immersion  tests  of  300M  steel.  Materials  studied  were 
300M  and  HP-9Ni-4Co-. 30  steels. 
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In  Task  II,  Landing  Gear  Stress  History,  a  flight-by-flight 
test  stress  history  was  prepared  for  a  landing  gear  component. 

The  F-4J  main  landing  gear  was  selected  to  determine  the  stress 
values  and  to  estimate  times  associated  with  each  stress 
condition.  This  stress  history  was  developed  based  on  design 
loads  and  not  on  field  measurements. 

In  Task  III,  Crack  Growth  Prediction,  the  methodology  devel¬ 
oped  and  revised  in  Task  I  was  used  to  predict  spectrum  crack 
growth  in  dry  air  and  alternate  immersion  in  sea  water 
environments.  Predictions  were  prepared  for  the  flight-by-flight 
stress  history  both  including  and  ignoring  sustained  loads. 

In  Task  IV,  Verification  Test  Program,  ten  tests  were  per¬ 
formed  to  verify  the  crack  growth  prediction  methodology. 

Spectrum  fatigue  tests  were  performed  on  center  cracked  and 
surface  flawed  panels  in  both  dry  air  and  synthetic  sea  water 
environments  using  both  accelerated  and  sustained  load  stress 
histories. 

In  Task  V,  Comparison  of  Analysis  and  Test  Results,  results  of 
the  verification  test  program  were  compared  with  predicted  crack 
growth  from  Task  III.  These  comparisons  indicate  that  the  predic¬ 
tion  methodology  provides  accurate  assessment  of  crack  growth 
under  aggressive  environments  under  varying  loading  conditions. 
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SECTION  II 

MODEL  CALIBRATION  AND  EVALUATION  PROGRAM 


1.  TEST  PROGRAM  SUMMARY  -  The  purpose  of  this  program  was  to 
obtain  experimental  data  necessary  for  calibrating  and  evaluating 
a  set  of  crack  growth  prediction  algorithms.  Table  1  summarizes 
the  crack  growth  program.  Two  steels  (300M  and  HP-9-4-.30)  were 
each  subjected  to  14  tests.  Data  were  obtained  to  determine 
material  behavior,  verify  stress  intensity  solutions  for  part- 
through  elliptical  flaw  geometries,  and  develop  crack  growth 
prediction  algorithms.  The  test  series  was  identical  for  each 
material  with  the  exception  of  Specimen  13  for  which  test 
conditions  were  selected  independently  in  each  case.  The  tests 
had  four  objectives: 

o  Characterize  each  material  by  developing  constant  amplitude 
fatigue  crack  growth  rate  data  in  a  dry  air  environment. 
(Specimens  1-3). 

o  Evaluate  a  linear  summation  model  of  environmentally 

accelerated  crack  growth.  (Specimens  4-11).  Specimens  4-6 
were  used  to  evaluate  the  effects  of  frequency  on  environ¬ 
mentally  accelerated  crack  growth  rates.  Specimens  7-11 
were  used  to  evaluate  wave  shape  effects  under  an  aggres¬ 
sive  environment  and  the  interaction  of  wave  shape  with 
frequency  and  stress  ratio. 

o  Calibrate  the  load  interaction  model  for  the  landing  gear 
stress  history  used  in  the  verification  test  program. 

Center  cracked  panel  specimens  were  tested  in  dry  air  to 
the  landing  gear  stress  history  without  sustained  loads 
(Specimen  12).  Specimen  13  was  used  to  duplicate  any  one 
test  in  each  material  which  appeared  to  give  results  which 
are  inconsistent  with  expected  trends. 

o  Evaluate  the  stress  intensity  solution  for,  semi-elliptical 
surface  flaw  specimens  used  in  verification  testing 
(Specimen  14). 

2.  TEST  MATERIALS  -  Materials  used  in  this  program  were  300M  and 
HP-9-4-.30  steels.  300M  steel  is  used  in  landing  gear  of  several 
Navy  aircraft.  These  materials  were  selected  to  represent  the 
range  of  environmental  sensitivity  found  in  landing  gear  steels. 
Crack  growth  in  300M  steel  is  significantly  accelerated  by 
exposure  to  salt  water  or  sea  water  under  loading  cycles  involving 
sustained  loads  (Reference  1).  Crack  growth  in  HP-9-4-.30  steel 
has  been  shown  to  be  relatively  insensitive  to  immersion  in  sea  or 
salt  water. 


NADC-79095-60 


TABLE  1.  MODEL  CALIBRATION  AND  EVALUATION  TEST  PROGRAM  SUMMARY 


Specimen 

Number 

Specimen  Type 

Environment 

Stress 

Ratio 

Frequency 

(cps) 

Wave  Shape 

Test  Type 

Objectives 

H 

Ce 

Cracks 

nter 
>d  Panel 

<10% 

RH  Air 

0 

0.5 

-1 

n 

Const 

-Amp 

Develop  da/dn  and 
Evaluate  Stress  Ratio 

4 

5 

6 

■ 

■a 

0 

0 

0 

10 

i 

0.1 

I 

I 

■ 

Evaluate  Frequency 

7 

8 

9 

10 

11 

0 

0 

0 

0.5 

-1 

10 

1 

0.1 

0.1 

0.1 

Evaluate  Wave  Shape 

12 

<10% 

RH  Air 

- 

10 

Sinusoidal 

Spec 

trum 

Calibrate  Model 

13 

- 

““ 1^ 

| 

Duplicate  Test 

14 

Elliptical  Flaw 

<10%  RH  Air 

0 

10 

Sinusoidal 

Const-Ainp 

Evaluate  K 

Test  series  is  identical  for  both  materials  GPO3-O0M-2 


Tensile  test  results  for  the  steel  materials  used  in  this 
program  are  presented  in  Table  2.  Three  tests  were  performed  for 
each  material.  The  strength  levels  found  are  typical  of  those 
used  in  these  materials  for  landing  gear  applications. 


TABLE  2.  TENSILE  TEST  RESULTS 


Alloy 

Specimen 

Yield 

Strength 

(kai) 

Ultimate 

Strength 

(kail 

Percent 

Elongation 

Percent 
Reduction 
of  Area 

Modulus  of 
Elaaticitv 
(pei  x  10&) 

1 

248.0 

296.0 

11 

37 

29.5 

300M 

2 

296.0 

11 

37 

29.2 

3 

248.0 

295.0 

11 

38 

28.5 

Average 

248.0 

295  5 

11 

37 

29.1 

1 

210.0 

228.0 

16 

59 

28.5 

HP  9-4-0.30 

2 

210.0 

229.0 

14 

56 

28.3 

3 

211.0 

229.0 

15 

62 

28.4 

Average 

210.5 

228.5 

15 

59 

28.4 

QP034UM 
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3.  TEST  SPECIMENS  -  The  specimens  were  of  two  types:  center 
cracked  panels  and  surface  flawed  panels.  Test  specimens  are 
shown  in  Figures  1  and  2. 


Dimensions  in  inches 


Figure  1.  Center  Cracked  Panel  Specimen 
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Figure  2.  Surface  Flawed  Panel  Specimen 


Surface  flaw  specimens  were  used  for  stress  intensity  calibra 
tion  of  part-through  cracks  and  to  determine  environment-load 
interaction  effects  on  part-through  flaws  representative  of  flaws 
found  in  landing  gear  structure  (Reference  2).  The  reduced 
section  (Figure  2)  was  required  to  properly  simulate  stresses 
experienced  in  landing  gear  components,  and  maintain  load  levels 
within  the  capacity  of  available  fatigue  test  equipment.  The 
predicted  finite  width  effect  on  crack  growth  was  small. 
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4.  TEST  PROCEDURES 

a.  Precracking  Procedures  -  Center  cracked  panel  specimens 
were  precracked  at  zero  stress  ratio  until  the  total  crack  length 
was  approximately  0.20  inch.  The  final  0.05  inch  of  crack  growth 
was  performed  at  a  stress  level  equal  to  or  less  than  that  at 
which  the  subsequent  test  was  performed.  For  alternate  immersion 
tests,  the  cracked  section  of  the  specimen  was  immersed  during  the 
final  0.05  inch  of  precracking.  Specimens  were  immersed  during 
the  final  stages  of  precracking  to  avoid  an  artificial  increase  in 
crack  growth  rate  found  to  occur  when  specimens  precracked  in  a 
laboratory  air  environment  are  immersed  in  salt  water  for  test 
(Reference  1). 

Semi-elliptic  surface  flaw  specimens  had  EDM  notches  having 
0.05  inch  surface  length  and  0.025  inch  depth  to  insure  that  the 
surface  flaw  would  initiate  at  all  points  along  the  EDM  slot 
periphery  during  precracking  to  0.1  inch  surface  length.  The 
initial  flaw  geometry  was  selected  to  insure  that  crack  growth  in 
the  depth  dimension  is  accompanied  by  growth  on  the  surface. 
Specimens  were  precracked  at  zero  stress  ratio.  The  final  0.025 
inch  of  precrack  growth  was  obtained  using  a  load  level  equal  to 
or  less  than  that  used  for  subsequent  testing.  Final  precracking 
was  performed  in  sea  water,  if  the  test  was  to  involve  the 
alternate  immersion  environment. 

b.  Specimen  Loading  and  Instrumentation  -  Cyclic  testing  was 
performed  in  an  MTS  test  system,  consisting  of  a  hydraulic  power 
supply,  load  frame  assembly,  electronic  control  console,  and 
tele-printer.  The  specimens  were  loaded  through  self-aligning 
hydraulic  grips.  Teflon  roller  guides  were  installed  against  the 
specimen  surface  to  prevent  buckling  during  application  of 
compression  loads. 


During  testing,  crack  lengths  were  optically  monitored  using  a 
linear  displacement  transducer  with  a  30X  microscope.  During 
constant  amplitude  tests  in  dry  air,  crack  growth  measurements 
were  made  after  every  .05  inches  of  crack  growth  (approximately). 
During  alternate  immersion  tests,  crack  lengths  were  measured 
before  immersion,  after  immersion,  and  at  the  midpoint  of  the  dry¬ 
ing  cycle  (25  minutes).  During  spectrum  tests,  crack  lengths  were 
measured  every  100  landings  until  2000  landings  were  obtained, 
measurements  were  made  every  500  landings  thereafter. 

c.  Environmental  Control  -  Two  environments  were  used  for 
tests  in  this  program:  low  humidity  dry  air  (10%  R.H.)  and  alter¬ 
nate  immersion  in  ASTM  standard  synthetic  sea  water  (10  minutes 
immersed  50  minutes  blown  dry  by  filtered  shop  air).  The  contain¬ 
ment  method  used  for  the  low  humidity  air  environment  is  shown  in 
Figure  3.  Low  humidity  air  was  obtained  using  silica  jell 
dessicant . 
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The  container  used  for  alternate  immersion  tests  is  shown  in 
Figure  4.  The  plexiglass  container  is  wide  enough  to  accept  the 
specimen  and  a  nylon  tube  which  encircles  the  specimen  to  provide 
air  for  the  drying  cycle.  Synthetic  sea  water  was  introduced 
through  a  plastic  tube  and  stainless  steel  port.  The  synthetic 
sea  water  used  was  ASTM  D  1141-75  substitute  ocean  water  without 
heavy  metals.  This  water  is  commonly  used  for  salt  water  spray 
tests . 


Figure  4.  Container  Used  for  Alternate  Immersion  Tests 


The  crack  was  totally  immersed  during  the  10  minute  immersion 
cycle.  The  sea  water  was  released  through  a  second  stainless 
steel  drain  port.  After  immersion  filtered  shop  air  at  10  psia 
was  used  to  blow  the  specimen  dry.  This  procedure  rapidly  dried 
the  specimen  surface  and  container.  Rapid  drying  provided  better 
control  of  the  environment  than  simply  allowing  the  specimen  and 
container  to  dry  in  laboratory  air. 
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5.  SUMMARY  OF  CONSTANT  AMPLITUDE  TEST  RESULTS  -  The  results  of 
constant  amplitude  crack  growth  testing  were  summarized  using  the 
procedure  shown  in  Figure  5.  They  are  presented  in  Table  3  in 
terms  of  number  of  cycles  to  grow  a  center  crack  over  the  range  of 
crack  lengths  shown  for  each  material. 


Cycles  - •- 

GP03-0838-6 


Figure  5.  Algorithm  Test  Program  Summary  Procedure 


TABLE  3.  SUMMARY  OF  CONSTANT  AMPLITUDE  LIVES 

Results  of  Model  Calibration  and  Evaluation  Tests 


Specimen 

Number 

Specimen  Type 

Environment 

Stress 

Ratio 

Frequency 

(cps) 

Wave  Shape 

Test  Type 

300M 

Cycles  to  Grow  Crack 
from  0.5  to  2.6  in. 

HP  9-4-0.30 
Cycles  to  Grow  Crack 
from  0.5  to  2.6  in. 

1 

2 

3 

Cer 

Cracke 

1 

ter 

d  Panel 

<  10% 

RH  Air 

0 

0.5 

-1 

10 

10 

10 

soidal 

Const 

-Amp 

140,000 

370,000 

122,000 

137,000 

390,000 

115,800 

4 

5 

6 

■ 

Alte 
lmm< 
in  Sea 

mate 

:rsion 

Water 

0 

0 

0 

10 

1 

0.1 

I 

. 

■ 

155,000 

94,000 

22,900 

131,000 

101,900 

22,500  A 

7 

8 

9 

10 

11 

1 

1 

1 

10 

1 

0.1 

0.1 

0.1 

Trape 

zoidal 

1 

140,000 

62,700 

6,210 

7,380 

4,010 

108,500 

149,400 

35,000  A 
76,500  A 
32,000  A 

Note:  These  flam  did  not  grow  until  precracked  beyond  1.25  in.  Lives  shown  are  based  on  1.70  in.  to  2.6  in.  of  growth  QPM4W7S 


2.  Maximum  gross  stress  in  each  test  was  20  ksi 
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6 .  CONSTANT  AMPLITUDE  FATIGUE  STRESS  RATIO  EVALUATIONS  - 
(Specimens  1-3)  -  Constant  amplitude  fatigue  tests  were  performed 
on  6  center  crack  panels  (three  per  material)  to  determine  the 
effect  of  stress  ratio  on  crack  growth  rate  and  to  verify  an 
analytical  stress  ratio  correction.  Tests  were  performed  at 
stress  ratios  of  0,  0.5,  and  -1  for  both  of  the  materials  in  a 
low-humidity  air  environment  (10%  R.H.).  The  wave  shape  used  for 
these  tests  was  a  sinusoid,  applied  at  a  frequency  of  10  Hz. 
During  the  tests,  crack  growth  was  visually  monitored  and  crack 
length  recorded  at  approximately  0.05  inch  intervals. 


The  results  of  the  tests  used  to  evaluate  the  effects  of 
stress  ratio  are  summarized  in  Table  3,  and  Figures  6  through  11. 
The  figures  also  present  predictions  based  on  an  analysis  proce¬ 
dure  discussed  in  Section  III.  Plots  of  crack  growth  rate  versus 
AK  (stress  intensity  factor  range)  are  presented  in  Figures  7  and 
10  to  show  the  effect  of  stress  ratio  on  growth  rate  as  commonly 
presented.  Figures  8  and  11  show  the  same  data  plotted  versus 
Kmax  so  that  crack  growth  rates  show  the  same  behavior  as  that 
shown  in  the  plots  of  crack  growth  measurements  (Figures  6  and  9) 
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Figure  6.  Effect  of  Stress  Ratio  on  Crack  Growth  in  300M  Steel 


da/dN  -  microinches/cycle 
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Figure  9.  Effect  of  Stress  Ratio  on  Crack  Growth  in  HP  9-4-0.30  Steel 
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Figure  11.  Effect  of  Stress  Ratio  on  Crack  Growth  Rate  in  HP  9-4-0.30  Steel 

(Plotted  vs  Kmax) 

7.  CONSTANT  AMPLITUDE  FATIGUE  FREQUENCY  EVALUATIONS  -  (Speci- 
mens  4-6)  -  Constant  amplitude  fatigue  tests  were  performed  on  6 
center  crack  panels  (3  per  material)  to  determine  the  effect  of 
cyclic  frequency  on  crack  growth  rate  and  to  verify  crack  growth 
predictions  obtained  with  a  linear  superposition  model  (Section 
III).  Specimens  were  tested  at  frequencies  of  10  cps,  1  cps,  and 
0.1  cps  for  each  material.  Each  specimen  was  subjected  to  an 
aggressive  environment  (alternate  immersion  in  sea  water)  during 
test.  The  wave  shape  was  sinusoidal.  During  the  tests,  crack 
growth  was  visually  monitored  and  crack  length  recorded  before 
immersion,  after  immersion,  and  25  minutes  into  the  drying  cycle. 
In  addition,  crack  length  measurements  were  recorded  at 
approximately  0.05  inch  intervals. 
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The  results  of  this  series  of  tests  are  summarized  in  Table  3 
(Page  10)  and  Figures  12  through  19.  300M  steel  shows  a  much 

larger  effect  of  test  frequency  in  the  sea  water  environment  than 
HP-9-4-.30.  The  300M  shows  a  factor  of  7  decrease  in  crack  growth 
life  as  the  frequency  is  reduced  from  10  cps  to  0.1  cps. 

HP-9-4-.30  shows  little  effect  of  either  environment  or  frequency. 
Crack  growth  life  for  HP-9-4-.30  varies  less  than  50  percent  as 
frequency  decreases  from  10  cps  to  0.1  cps.  At  0.1  cps  and  stress 
intensity  factors  less  than  25  ksi  / in ,  crack  growth  in  HP-9-4-.30 
stopped  during  the  drying  portion  of  the  alternate  immersion 
cycle.  This  phenomenon  occurred  in  trapezoidal  wave  tests  also 
and  is  discussed  in  detail  in  the  next  section. 


Figure  12.  Effect  of  Sine  Wave  Loading  Frequency  on  Crack  Growth  in  300M  Steel 
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da/dN  -  microinches/cycle 
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A  K  -  ksi  V  in. 

QP03-043A-74 

Figure  14.  Effect  of  Sine  Wave  Loading  at  1  cps  on  Crack  Growth 
Rate  in  300M  Steel  in  Sea  Water 
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Figure  15.  Effect  of  Sine  Wave  Loading  at  0.1  cps  on  Crack  Growth 
Rate  in  300M  Steel  in  Sea  Water 
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Figure  17.  Effect  of  Sine  Wave  Loading  at  10  cps  on  Crack  Growth 
Rate  in  HP  9-4-0.30  Steel  in  Sea  Water 
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Figure  19.  Effect  of  Sine  Wave  Loading  at  0.1  cps  on  Crack  Growth  Rate  in 
HP  9-4-0.30  Steel  in  Sea  Water 


Because  of  the  slow  crack  growth  rate  in  HP-9-4-.30  at  0.1 
cps,  many  immersion  cycles  are  applied  prior  to  obtaining  apprec 
able  crack  growth.  Plots  of  crack  growth  and  growth  rate  are 
based  on  data  obtained  after  each  0.05  inch  of  crack  growth. 
Since  each  point  represents  many  environmental  cycles,  crack 
growth  rates  during  immersion  cannot  be  separated  from  total 
growth. 
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8.  CONSTANT  AMPLITUDE  FATIGUE  WAVE  SHAPE  EVALUATIONS  -  (Specimens 
7-11)  -  Ten  constant  amplitude  fatigue  tests  were  performed  on 
center  crack  panels  to  assess  the  interaction  of  environment, 
stress  ratio,  cyclic  frequency,  and  wave  shape  on  crack  growth. 
Five  tests  were  performed  on  specimens  from  each  material  using  a 
trapezoidal  wave  form,  Figure  20,  in  which  the  load  rate  allowed  a 
great  portion  of  the  cycle  to  be  held  at  a  constant  peak  stress. 
Tests  with  the  trapezoidal  wave  shape  were  run  at  frequencies 
ranging  from  0.1  cps  to  10  cps.  Tests  at  0.1  cps  were  run  at 
stress  ratios  of  0,  0.5,  and  -1.  All  specimens  were  subjected  to 
the  aggressive  environment  (alternate  immersion  in  sea  water) 
during  the  test. 


Maximum  Load 


Minimum  Load 


GP03-08 38-28 

Figure  20.  Trapezoidal  Wave  Shape 


The  effect  of  sustained  load  duration  are  presented  in  Figures 
21-28.  In  300M  steel,  increasing  the  sustained  load  time  in  the 
alternate  immersion  environment  shortened  life  considerably. 

Crack  growth  life  under  trapezoidal  wave  loading  at  0.1  cps  is 
about  five  percent  of  that  for  10  cps  loading. 
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Figure  21.  Effect  of  Trapezoidal  Wave  Loading  Frequency  on  Crack  Growth  in  300M  Steel 


26 


da/dN  -  microinches/cycle 


4 


N ADC -79095-60 


5  10  20  30  50  100 


A  K  -  ksi  V  in.  gpo3<»3«-69 


Figure  22.  Effect  of  Alternate  Immersion  Environment  on  Crack  Growth  Rate 
on  300M  Steel  ■  10  cps  Trapezoidal  Wave  Loading 
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Figure  23.  Effect  of  Alternate  Immersion  Environment  on  Crack  Growth  Rate 
in  300M  Steel  ■  1  cps  Trapezoidal  Wave  Loading 
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Figure  24.  Effect  of  Alternate  Immersion  Environment  on  Crack  Growth 
in  300M  Steel  ■  0.1  cps  Trapezoidal  Wave  Loading 
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Stress  ratio  -  0 
Cyclic  rate  -  1  cps 
Wave  shape  -  trapezoidal 
Environment  -  alternate  immersion 
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Figure  27.  Effect  of  Trapezoidal  Wave  Loading  at  1  cps  on 
Crack  Growth  Rate  in  HP  9-4-0.30  Steel 
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Figure  28.  Crack  Growth  Rate  in  HP  9-4-0.30  Steel  in  Alternate  Immersion 
Environment  Using  0.1  cps  Trapezoidal  Wave  Loading 


HP-9-4-.30  steel  does  not  show  the  environmental  acceleration 
found  in  300M  steel,  and  consequently  does  not  show  much  effect  of 
frequency  or  environment.  However,  under  the  0.1  cps  trapezoidal 
wave  test  of  HP-9-4-.30  specimen  9  (R=0),  the  threshold  stress 
intensity  factor  in  the  alternate  immersion  environment  was  found 
to  increase  to  about  25  ksi  /In.,  Figure  28.  Below  this  thresh¬ 
old  value,  the  specimen  was  precracked  in  sea  water  in  0.05  inch 
increments  and  then  subjected  to  the  alternate  immersion  environ¬ 
ment  and  the  0.1  cps  trapezoidal  wave  loading;  the  crack  stopped 
growing  once  the  specimen  was  dried  and  did  not  start  growing 
during  immersion.  Thousands  of  precracking  cycles  were  required 
to  reinitiate  the  crack  growth  when  precracking  in  salt  water. 
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Currently,  our  conjecture  is  that  corrosion  products,  created 
on  the  crack  faces  during  immersion,  solidify  during  drying. 

These  corrosion  products  reduce  the  stress  intensity  factor  ranges 
applied  to  the  crack  tip  (i.e.  increasing  crack  closure)  causing 
the  crack  growth  to  slow  or  stop.  During  immersion,  these 
corrosion  products  are  formed  but  are  kept  in  solution  in  the  sea 
water  and  do  not  significantly  affect  the  crack  tip  behavior. 

This  would  explain  the  negligible  impact  of  continuous  immersion 
environment  on  crack  growth  in  HP-9-4-.30  steel  found  in  previous 
work.  Reference  1. 

In  the  alternate  immersion  environment,  high  cycle  rates  may 
tend  to  reduce  the  buildup  of  corrosion  products  near  the  crack 
tip  during  drying,  producing  little  impact  on  crack  growth  at  low 
stress  intensity  factors.  At  low  cycle  rates,  the  crack  tip  is 
held  open  long  enough  during  drying  that  corrosion  products 
solidify  and,  acting  like  a  wedge  in  the  crack,  prevent  growth. 
Evidently,  the  ten  minute  immersion  period  used  in  these  tests  is 
not  sufficient  to  dissolve  the  corrosion  products.  At  the  thresh¬ 
old  stress  intensity  factor  range,  the  crack  begins  to  grow,  in 
spite  of  the  corrosion  product  wedge.  Once  crack  growth  begins 
the  crack  tip  grows  away  from  the  corrosion  products  during  the 
drying  cycle  and  crack  growth  rates  rapidly  approach  those  found 
in  dry  air. 

The  combined  effects  of  wave  shape,  stress  ratio,  and  alter¬ 
nate  immersion  environment  are  demonstrated  in  the  crack  growth 
results  presented  in  Figures  29  through  34.  The  data  presented  in 
Figures  29  and  32  show  that  the  effects  of  stress  ratio  under 
trapezoidal  wave  shape  loading  are  somewhat  less  for  300M  than  for 
HP  9-4-. 30. 

Stress  ratio  effects  for  HP  9-4-. 30  steel  are  relatively 
insensitive  to  environment  and  load  frequency.  This  is  demon¬ 
strated  by  the  agreement  between  the  crack  growth  rate  data  and 
analysis  shown  in  Figures  33  and  34.  The  analysis  results  shown 
are  based  on  dry  air  data.  Also  shown  in  Figures  33  and  34  is  the 
variation  of  threshold  stress  intensity  factor  range  (AK),  with 
stress  ratio.  The  little  data  taken  at  lowAK's  indicates  that 
crack  growth  rates  in  the  alternate  immersion  environment  are  very 
close  to  those  in  dry  air  whenever  Kmax  exceeds  about  25  ksi  /  in 
( AK  =  12.5  ksi  /  in  at  R  =  0.5,  AK  =  50  ksi  / in  at  -1.0). 

High  frequency  trapezoidal  wave  loading  has  little  effect  on 
crack  growth  life  in  either  steel  in  an  alternate  immersion 
environment  (Figures  35  and  36).  At  low  frequencies,  both  sine 
and  trapezoidal  wave  loadings  have  significant  effect  on  crack 
growth  behavior  (Figures  37  and  38).  In  HP  9-4-. 30  low  frequency 
loadings  of  either  wave  shape  in  alternate  immersion  tests  raised 
the  threshold  stress  intensity  factor  above  the  initial  stress 
intensity  factors  used  in  the  high  frequency  tests  (Figures  25  and 
26,  pages  30  and  31). 
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Test  data 
A  Immersed 
D  1st  half  drying  cycle 
O  2nd  half  drying  cycle 
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Environment  -  alternate  immersion 
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Figure  30.  Crack  Growth  Rate  for  300M  Steel  in 
Alternate  Immersion  Environment,  R  =  0.5 
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Figure  31.  Crack  Growth  Rate  for  300M  Steel  in 
Alternate  Immersion  Environment,  R  =  -  1.0 
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Figure  34*  Crack  Growth  Rate  in  HP  9-4-0.30  Steel  in  Alternate 
Immersion  Environment,  R=  -1.0 
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Figure  36.  Effect  of  Environment  and  Wave  Shape  on 
Crack  Growth  in  HP  9-4-0.30  Steel  at  10  cps 
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Half  Crack  Length  -  in. 


Figure  38.  Effect  of  Wave  Shape  on  Crack  Growth  in  HP  9-4-0.30  Steel 
in  Sea  Water  at  0.1  cps 
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In  300M  steel,  at  low  frequency,  a  trapezoidal  wave  causes 
faster  crack  growth  than  a  sine  wave  at  a  similar  frequency.  As 
shown  in  Table  3  (page  10)  the  crack  growth  life  for  a  0.1  cps 
trapezoidal  wave  is  about  one  fourth  that  under  an  0.1  cps 
sinusoidal  wave.  This  is  expected  since  300M  exhibits  large 
sustained  load  cracking  rates  in  a  salt  water  environment 
(References  1,  3,  and  4). 

9.  SPECTRUM  AND  DUPLICATE  TESTS  -  (Specimens  12-13)  -  Specimens 
12  and  13  were  used  to  evaluate  crack  growth  retardation  under 
landing  gear  spectrum  loads  in  center  crack  panels.  Two  tests  in 
each  material  were  performed  in  the  dry  air  environment  using  the 
spectrum  defined  in  Section  V.  The  test  frequency  was  approxi¬ 
mately  10  cps. 

In  HP-9-4-.30,  spectrum  tests  were  performed  at  40  ksi  and  48 
ksi  maximum  stress  levels.  The  results  presented  in  Figure  39 
indicate  the  sensitivity  of  crack  growth  lives  to  stress  level 
variations.  Analysis  results  indicate  that  the  effect  on  life  of 
such  variations  is  adequately  predicted.  Some  accelerated  growth 
can  be  noticed  in  both  results  prior  to  the  occurrence  of  the 
maximum  stress  level  at  about  3800  landings.  This  is  probably  due 
to  the  low  precracking  stress  level  (12  ksi)  used  for  these  tests. 


Figure  39.  Spectrum  Test  Results  for  HP  9-4-0.30  in  Dry  Air 
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In  300M  steel,  both  spectrum  tests  were  performed  at  48  ksi 
and  two  initial  crack  lengths  were  used  (1.0  inch  and  1.14  inch). 
Results  from  both  tests  are  presented  in  Figure  40.  The  first 
test  took  20,000  landings  to  reach  a  crack  of  1.14  inches.  Data 
from  the  second  test  are  plotted  starting  at  this  point  to  sim¬ 
plify  comparison  of  the  two  test  results.  The  results  presented 
in  Figure  40  indicate  the  specimen-to-specimen  variation  in  crack 
growth  under  the  landing  gear  spectrum.  Reasonably  good  correla¬ 
tion  between  analysis  and  both  test  results  was  obtained.  Initial 
acceleration  is  not  evident  in  these  results. 


Figure  40.  Spectrum  Test  Results  for  300M  in  Dry  Air 
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These  test  results  permitted  selection  of  overload  interac-  j 

tion  zone  sizes  for  both  materials,  which  were  used  to  predict  the  i 

crack  growth  behavior  expected  in  the  verification  test.  The  j 

prediction  methodology  is  discussed  in  Section  III.  ’j 

10.  SURFACE  FLAW  TESTS  (Specimen  14)  -  Two  tests  of  surface  j 

flawed  panels  were  used  to  verify  the  crack  growth  analysis  I 

procedure  for  surface  flaws.  One  constant  amplitude  test  in  each  ; 

material  was  performed  in  the  dry  air  environment.  The  test 

frequency  was  approximately  10  cps.  Tests  were  intended  to  last  i 

approximately  10,000  cycles  at  100  ksi  and  crack  lengths  were  ’ 

measured  after  500  cycle  increments.  ’ 

Comparisons  of  predicted  and  measured  crack  growth  on  the  j 

specimen  surface  are  presented  in  Figures  41  and  42.  Predicted 

and  measured  growth  rates  are  compared  with  center  cracked  panel  I 

growth  rates  in  Figures  43  and  44.  Measured  growth  rates  in  both  j 

tests  are  initially  lower  than  predicted.  In  300M  the  growth  rate 
rapidly  approaches  predicted  values  and  good  agreement  between 

predicted  and  measured  growth  is  obtained  (Figure  41).  In  HP  I 

9-4-. 30  the  initial  growth  rates  are  low  enough  that  the  crack 
growth  curve  comparison  of  Figure  42  is  not  as  good  as  for  30 0M. 

In  both  cases  the  initial  low  growth  rates  measured  may  be  due  to  i 

differences  between  assumed  and  actual  initial  flaw  shapes.  These 

differences  rapidly  disappear  as  the  flaw  shapes  approach  natural 

flaw  shapes.  The  generally  good  agreement  between  predicted  and 

measured  results  indicates  the  accuracy  of  the  surface  flaw 

analysis  procedure. 


NADC-79095-60 


a aaa aaaaa  a aaaa a aaaa  aaaaa aaaaa  *■ 

aaaa  aaaaaaaaaa  INHUIli  . 

aaaaaaaaaaaaaaaaaaa  aaaaaaaaaf  •  -a a ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

aaaaaaaaa aaaaaaaaaa aaaaaaaa'  aaaaa 

aaaaaaaaa  aaaaaaaaaa  aaaa  ».»-  •■■■■■■■  aa aaaaa  aaaaa aaaaaaaaaa  aaaaaaaaaa  aaaaa 


!■■■■■■■ ■■■■■■■■■•■■■•■■■■■■■■■■■■■■■■■■■■■ III . ■«■•■■■ 

.  . . ■■■■ . «■■■■■  aaaaa aaaaa  aaaaa 

■  aaa . . . . . . 

aaaaaaiaaaaiaaaaiaaaaaaaaaiaaaaaiBaaiaaaaaaaaat •aaaiiaaaaaaaaa 

. . .  . . . 

. . . . .  aaBBBaaaaa  matataai  . . a  aaaaa 

■■■■■■•■a  aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa  aaBaaaaaaaiaaBa>aaaa aaaaa aaaaa  aaaaa 
aaaa  aaaaa  aaaaa aaaaa  aaaaa aaaaa  aataaaaaaaiBaBBaaaai aaaaa  aaaaa aaaiiaaaaa  aaaaa 
aaaaaiaaa ■■■■■■aaaa  aaaaa aaaaa  aaaaa aaaaa  aaaaa aaaaa  aaaaa aaaaa  aaaaa aaaaa  aaaaa 
aaaa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaaaaaaa aaaaa aaaaa aaaaaaaaaa aaaaa ■ 

aaaa . aaaaa aaaaa ■aaaSaaaSa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa aaaaa  aaaaa 

aaaa aaaaa aaaaa aaaaa ■aaaifaaai aaaaa fjfaaiaaaaaaaaa aaaaa aaaaa aaaaaaaaaa aaaaa 


aaaa aaaaa laaiaiaaaa aaaaa aaa 
aaaa  aaaaa  aaaaauMUMauU 
aaaa  aaaaa aaaaa 
aaaaaaaaa  aaaaa  MMIIMai 
■aaa aaaaa aaaaa aaaaa aaaaa aaa 
aaaa  aaaaa  aaaaa  aaaMMMAaaaj 


NADC-79095-60 


Thousands  of  Cycles  qpo3-o»m-2i 


Figure  42.  Comparison  of  Surface  Flaw  Test  and  Analysis 
Results  in  HP  9-4-0.30  Steel 
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da/dN  -  microinches/cycle 


A  K  -  ksi 
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Figure  44.  Comparison  of  Surface  Flaw  and  Center  Crack  Growth  Rates  for  HP  9-4-0.30 
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SECTION  III 

ENVIRONMENT-LOAD  INTERACTION  MODEL 


1.  SUMMARY  -  The  analysis  of  environmental  acceleration  of  crack 
growth  is  based  on  a  linear  superposition  of  cyclic  load  and  sus¬ 
tained  load  crack  growth.  Stress  ratio  effects  are  accounted  for 
by  using  an  effective  stress  intensity  range  concept  based  on 
closure.  Crack  growth  analysis  using  the  combination  of  these  two 
algorithms  shows  good  correlation  with  the  results  of  the  model 
evaluation  and  calibration  tests. 

The  Willenborg  model  was  the  basis  for  extension  of  the  crack 
growth  analysis  to  spectrum  loadings.  Spectrum  test  results  from 
center  crack  panels  were  used  to  calibrate  the  model  for  predic¬ 
tions  of  surface  flaw  and  center  crack  growth  under  spectrum 
loadings  for  the  verification  tests.  Stress  intensity  factor 
solutions  for  semi-elliptical  surface  flaws  were  determined  using 
a  slice  synthesis  model,  validated  by  comparisons  with  finite 
element  solutions. 

2.  STRESS  RATIO  EFFECTS  -  The  method  used  herein  to  predict  the 
effect  of  stress  ratio  on  crack  growth  rate  is  based  on  analysis 
of  crack  closure  caused  by  crack  surface  displacements.  Crack 
closure  determines  the  stress  intensity  range  which  is  effective 
in  propagating  the  crack.  This  effective  stress  intensity  range 
is  obtained  from  a  curve  fit  to  the  analyses  of  Reference  5. 

4Ke££  =  <£f°>  11  -  <l-R)  *•]  for  R  >  0  (1) 


iKe££  ‘  !1  *  K°e°'1R)  Kmax  £or  R  ‘  0  '2> 


where  R  and  Kmax  are  defined  by  the  remote  loading  conditions  and 
K 0  is  the  ratio  of  closure  stress  intensity  to  Kmax  at  R=0, 
generalized  as 

K°  =  (0.33045  +  0.15164  r  -  0.01476  rp2 )  (3) 

[1  +  0.6  (fc/fm  "  1)  -  0.156  (fc/fm  -  l)2] 


where 

fm  =  monotonic  yield  stress 
f c  =  cyclic  yield  stress 
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The  parameter,  rp,  was  used  to  "tune"  Equations  (1)  and  (3)  for 
each  material  analyzed  so  that  constant  amplitude  test  lives  for 
R=0  and  R=0.5  were  matched.  In  the  closure  analysis  of  Reference 
5,  rp,  was  related  to  plastic  zone  conditions,  rp=l  for  plane 
stress,  rp=0  for  plane  strain.  As  presented  in  Table  4,  rp=0 
(Ko=0.33)  was  found  to  provide  good  correlation  with  the  R=0  and 
R=0.5  crack  growth  results  in  both  300M  and  HP  9-4-. 30.  This 
differs  somewhat  from  the  values  found  in  Reference  1  where  Kq  for 
300M  was  0  and  for  HP  9-4-. 30  was  0.242.  However,  because  of  heat 
treating  difficulties  noted  in  Reference  1,  the  material  used  in 
that  program  was  less  well  behaved  than  that  used  in  this  program: 
This  probably  accounts  for  the  differences  in  behavior  at  various 
stress  ratios. 


TABLE  4.  PARAMETERS  USED  FOR  STRESS  RATIO  CORRECTIONS 


Material 

fpA 

A 

ksi 

k°  A 

HP  9-4-0.30 

0 

210 

210 

0.330 

300M 

0 

248 

228 

0.330 

Notes: 


A 


Plastic  zone  size  correction  factor 
Monotonic  yield  stress 
Cyclic  yield  stress 

Ratio  of  Kmax-eff  -  ^eff  ,0  Kmax-applied 


for  R  =  Q,  constant  amplitude  loading 
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The  closure  based  approach,  as  empirically  used  in  this  pro¬ 
gram,  has  two  advantages  over  the  Forman  equation,  Reference  6: 

(1)  predicted  effects  of  stress  ratio  on  crack  growth  rate  can  be 
varied  to  account  for  variations  in  material  behavior,  and  (2) 
negative  stress  ratios  generally  are  predicted  to  accelerate  crack 
growth.  As  shown  in  Figures  6-9  (Pages  11-14),  this  analysis 
agrees  well  with  the  test  data  from  this  program. 

3.  LINEAR  SUPERPOSITION  APPROACH  FOR  PREDICTION  OF  ENVIRONMENTAL 
ACCELERATION 


a.  Modification  of  Wei-Landes  Approach  -  The  Wei-Landes 
linear  summation  hypothesis,  Reference  7,  suggests  that  environ¬ 
mentally  accelerated  crack  growth  rates  can  be  predicted  by  adding 
the  crack  growth  due  to  the  individual  mechanisms  of  environmental 
attack  and  cyclic  loading.  This  hypothesis,  when  used  to  account 
for  frequency  effects,  can  be  expressed: 

da  _  1  da  +  da  ^ 

^total  ^  ^environment  ^Nfatigue 
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The  da/dtenvironment  is  obtained  from  sustained  load  or  low 
frequency  cyclic  tests  in  the  environment,  and  da/dN£at^gue  is 
obtained  from  cyclic  load  tests  conducted  at  high  frequencies  in 
air  or  a  non-aggressive  environment. 


The  modification  of  the  linear  summation  hypothesis  used  in 
this  program  is  that  for  any  environment  the  total  crack  growth 
rate  for  any  cycle  is  the  sum  of  two  components;  one  due  to 
mechanical  loading  in  the  reference  environment  (less  than  ten 
percent  relative  humidity),  and  one  due  to  sustained  load  in  an 
aggressive  environment.  This  model  is  formulated  as: 


da  _  da  I 
dN  dN 


+ 

ref 


da 

dN 


env 


(5) 


where , 


da 

ft2  da 

dN 

j 

env 

t  dt 

immersed 

immersed 

1 

da 

f2  da 

-at 

dN 

env 

1.  dt 

'“I 

e 

immersed 

drying 

1 

(during  immersion) 


(during  drying) 


(6a) 


(6b) 


Where  the  load  period  is  t2~tj_. 


The  sustained  load  crack  growth  rate  during  immersion  is 
determined  from  0.1  cps  trapezoidal  wave  test  data  by: 


da 

dt 


% 

immersed 


1 

10 


da 

immersed 
0.1  cps 


da 

dN 


ref  I 
10  cpsj 


(7) 


Equation  (6a)  was  numerically  integrated  for  a  1  cps,  R=0,  sine 
wave  loading  to  determine  the  crack  growth  rate  for  sine  waves  of 
various  Kmax  values.  Only  crack  growth  rates  during  the  loading 
portion  of  the  sine  wave  are  integrated  because  it  is  assumed  that 
the  unloading  portion  of  the  cycle  contributes  no  growth. 
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During  the  drying  cycle,  the  environmental  acceleration  decays 
exponentially  with  time  after  immersion.  This  results  in  a  stair¬ 
step  crack  growth  curve  at  low  frequencies.  Figure  45.  Using  the 
exponential  decay  function,  Equation  (6b),  the  crack  growth  rate 
for  a  cycle  during  the  drying  period  having  sustained  load  from 
time,  ti  to  t2 /  can  be  approximated  from  integration  of  Equation 
(6b)  as 


da 

dN. 

drying 


da  1  da 

dN  ,  a  f  d  t . 

ret  immersed 

10  cps 


-at 


(e 


(8) 


ere  dt  immersed  is  assumed  not  to  change  significantly  during  a 
single  cycle.  The  exponent,  a  ,  in  Equation  (6b)  was  determined 
empirically  by  correlation  with  the  crack  growth  rates  measured 
during  the  drying  cycle  of  the  trapezoidal  wave  test  using  R=0  at 
0.1  cps  (Specimen  9).  In  300M  steel  ct  was  found  to  be  0.00175. 

In  HP-9-4- . 30 ,  no  environmental  acceleration  was  found  and  a  is 
zero. 


Number  of  Cycles  N 
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Figure  45.  Crack  Growth  Behavior  of  300M  Steel  in  Alternate  Immersion 
in  Sea  Water  at  Low  Cycle  Rates 
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The  prediction  of  sustained  load  crack  growth  by  Equation  (6a) 
was  found  to  adequately  describe  crack  growth  behavior  in 
continuous  immersion  tests  (Reference  1).  However,  in  alternate 
immersion  tests,  the  sustained  load  crack  growth  rate  during 
immersion  was  less  than  that  predicted  by  Equation  (6a)  for  1  cps 
and  10  cps  load  frequencies.  Sustained  load  crack  growth  rates 
for  both  alternate  immersion  and  continuous  immersion  tests  are 
compared  in  Figure  46.  An  empirical  correction  is  applied  to 
Equation  (6a)  to  improve  correlation  with  the  test  results.  The 
corrected  expression  is: 


da 

dt 


immersed 


(0.1th) 


0.25  da 
dt 


immersed 
0.1  cps 


(9) 


where  t^  is  the  hold  time  in  a  given  load  cycle. 


Frequency  •  cps 
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Figure  46.  Comparison  of  Sustained  Load  Crack  Growth  Rates  in 
300M  Steel  for  Alternate  and  Continuous  immersion  Environments 


Sustained  load  cracking  during  immersion  is  assumed  to  occur 
only  during  the  loading  portion  of  the  cycle.  When  a  trapezoidal 
wave  is  applied,  predicted  crack  growth  is  the  sum  of  three 
components;  cyclic  growth  during  loading,  sustained  load  growth 
during  loading  (analyzed  as  a  sine  wave),  and  sustained  load 
growth  during  the  hold  time.  The  unloading  portion  of  a  cycle  is 
assumed  to  produce  no  crack  growth  and  is  ignored  in  the  analysis. 
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This  model  was  used  for  all  analyses  presented  herein.  Compar¬ 
ison  of  analysis  and  test  results  for  300M  steel.  Figures  21-24, 
and  45,  indicate  the  accuracy  of  the  model  to  describe  the  crack 
growth  behavior  of  300M  steel  in  the  alternate  immersion  environ¬ 
ment  at  several  loading  frequencies.  Predicted  and  measured  lives 
for  all  constant  amplitude  crack  growth  tests  in  300M  steel  are 
summarized  in  Table  5.  Good  correlation  is  evident. 


TABLE  5.  COMPARISON  OF  TEST  AND  ANALYSIS  RESULTS  IN  300M  STEEL 


Specimen 

Number 

Specimen  Typt 

Environment 

Stress 

Ratio 

Frequency 

(cps) 

Wave  Shape 

Test  Type 

3 DOM  Test 

Cycles  to  Grow  Crack 
from  O.S  to  2.1  in. 

mu  •  lelysis 
Cycles  te  .  e  Crack 
train  15  t*  '  1  in. 

1 

Center 

<  10% 

RH  Air 

0 

10 

Sinusoidal 

Const-Amp 

140.000 

139  000 

2 

Cracked  Panel 

0.5 

10 

370.000 

378  000 

3 

-1 

10 

122.000 

126.000 

4 

Alternate 

0 

10 

155,000 

135.000 

5 

Immersion 

0 

1 

94,000 

104.000 

6 

In  Sea 

Water 

0 

0.1 

22,900 

22.500 

7 

0 

10 

Trapezoidal 

140,000 

132,000 

8 

0 

1 

62.700 

59.300 

9 

0 

0.1 

6,210 

6.350 

10 

0.5 

0.1 

7,380 

6,630 

11 

-1 

0.1 

4.010 

6,320 

Note  Maximum  stress  level  was  20  ksi  for  each  test. 


The  effect  of  stress  ratio  on  life  in  300M  steel  in  the 
alternate  immersion  environment  is  indicated  in  Figures  29-31. 

The  correlation  of  predicted  and  measured  lives  in  the  alternate 
immersion  environment,  Figure  29,  is  not  as  good  as  that  of  the 
other  comparisons  shown.  Evidently  stress  ratio  has  more  impact 
on  crack  growth  in  the  aggressive  environment  than  is  computed 
from  Equation  (5).  The  correlation  is  considered  adequate  since 
the  analysis  does  predict  the  more  significant  reduction  in  life 
produced  by  the  aggressive  environment. 

Predicted  and  measured  crack  growth  lives  for  HP-9-4-.30  steel 
are  compared  in  Table  6.  These  predictions  assume  that  there  is 
no  impact  of  the  aggressive  environment.  This  assumption  is  con¬ 
firmed  by  the  test  data.  However,  the  0.1  cps  trapezoidal  wave 
shape  test  results  show  a  significant  increase  in  threshold  stress 
intensity  factor.  Plots  summarizing  the  effects  of  frequency  and 
stress  ratio  on  crack  growth  life  in  HP-9-4-.30  steel  are 
presented  in  Figures  25  and  32.  These  plots  indicate  the  accuracy 
of  the  technique. 
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TABLE  6.  COMPARISON  OF  TEST  AND  ANALYSIS  RESULTS  IN  HP  9-4-0.30  STEEL 


Specimen 

Number 

Sptcimen  Type 

Environment 

Stress 

Ratio 

Frequency 

(cps) 

Wave  Shape 

Test  Type 

HP  9-4-0.30  Test 
Cycles  to  Grow  Crack 
from  0.5  to  2.5  in. 

HP  9-4  0.30  Analysis 
Cycles  to  Grow  Crack 
from  0.5  to  2.6  in. 

i 

Center 

<  10% 

RH  Air 

0 

10 

Sinusoidal 

Const  Amp 

137,000 

139,000 

2 

Cracked  Panel 

0.5 

10 

390,000 

370.000 

3 

-1 

10 

115,800 

127,000 

4 

HU 

Alternate 

0 

10 

131,000 

5 

■ 

Immersion 

0 

1 

101,900 

13^,000 

6 

in  Sea 

Water 

0 

0.1 

H 

22,500 

32,000  A 

7 

0 

10 

Trapezoidal 

108,500 

139,000 

8 

0 

1 

149,400 

139,000 

9 

0 

0.1 

35,000 

32,000  A 

10 

0.5 

0.1 

76,500 

67,000  A 

11 

-1 

0.1 

32,000 

29,000  A 

Notes:  /l\  Cracks  did  not  grow  until  precracked  beyond  1 .25  inches.  Lives  presented  here 
A  are  from  1.7  to  2.6  inches. 

/2\  Maximum  stress  level  was  20  ksi  for  each  test. 


b.  Simple  Approach  to  Load-Environment  Interaction  Analysis  - 
A  less  complicated  analysis  technique,  presented  in  Reference  1, 
which  does  not  account  for  the  decay  in  crack  growth  rates  after 
immersion  in  an  aggressive  environment,  was  also  evaluated  using 
the  test  data  presented  in  Section  II.  This  second  technique  is 
similar  to  that  presented  previously  but  uses  average  crack  growth 
rate  data  from  the  0.1  cps,  alternate  immersion  tests  to  determine 
the  sustained  load  crack  growth  rate  via  Equation  (7).  in  this 
case,  the  alternate  immersion  cycle  in  total  is  considered  the 
aggressive  environment  and  computation  of  the  acceleration  due  to 
the  environment  is  made  using  only  Equation  (6a) . 

Comparison  of  analysis  results  using  this  technique  with  alter¬ 
nate  immersion  test  data  for  300M  steel  at  three  different  frequen¬ 
cies  is  presented  in  Figure  47.  Correlation  is  even  better  than 
that  for  the  more  detailed  analysis  shown  in  Figure  21,  Page  26. 
Correlation  with  the  remainder  of  constant  amplitude  tests  in  300M 
and  HP-9-4-.30  is  as  good  as  that  shown  for  the  more  detailed 
analysis  technique  previously  presented.  Thus,  the  simple  load- 
environment  interaction  analysis  presented  in  Reference  1,  while 
it  is  not  capable  of  predicting  detailed  crack  growth  behavior  in 
alternate  immersion  tests,  is  an  accurate  method  for  predicting 
the  gross  behavior  in  constant  amplitude  tests. 
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Figure  47.  Comparison  of  Results  of  Simple  Load-Environment 
Analysis  With  Crack  Growth  Behavior  of  300M  Steel 


The  average  crack  growth  rate  curve  used  for  these  analyses  is 
compared  with  the  test  data  for  0.1  cps,  alternate  immersion  test 
of  300M  steel  in  Figure  48.  This  comparison  indicates  the  inabil¬ 
ity  of  this  technique  to  account  for  the  transitory  acceleration 
which  occurs  in  alternate  immersion  tests.  This  is  not  a  major 
handicap  for  analysis  of  constant  amplitude  behavior,  where  many 
load  cycles  and  immersion  cycles  are  applied  during  the  specimen 
life.  However,  in  spectrum  crack  growth  analyses,  where  the  coup¬ 
ling  of  high  sustained  loads  and  immersion  in  salt  or  sea  water 
provide  significant  acceleration  of  crack  growth  rates  in  300M 
steel,  this  simple  technique  could  provide  erroneous  predictions 
of  crack  growth  behavior.  Therefore,  the  simple  technique  was  not 
used  for  prediction  of  crack  growth  behavior  under  spectrum 
loadings  in  this  program. 
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A  K  -  ksi  y/in^ 

QP034M3S-56 

Figure  48.  Comparison  of  Average  da  /dN  Curve  with  Crack  Growth  Rate  Test  Data 


c.  Spectrum  Analysis  with  Environmental  Acceleration  - 
Application  of  the  crack  growth  algorithms  to  spectrum  analysis 
was  straightforward.  The  spectrum  was  input  as  a  series  of 
stresses  and  corresponding  time  increments  to  peak  stress  (Section 
IV)  .  Stress  levels  were  joined  by  haversine  waves  as  shown  in 
Figure  49.  Sustained  stress  levels  are  represented  by  two  adjoin¬ 
ing  peaks  having  the  same  stress  level,  the  peak  to  peak  time 
being  the  hold  time. 
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Figure  49.  Spectrum  Load  Definitions  Used  for  Analysis 


The  stress  history  was  simplified  so  that  it  contains  only 
peaks  and  valleys  and  the  times  at  which  they  are  applied.  Thus, 
during  simplification,  the  stress  level  sQ  is  the  first  valley, 
stress  level  s^  is  dropped  and  S2  is  the  peak  for  the  first  half 
cycle.  The  time  increment  for  the  first  cycle  is  the  sum  of  ti 
and  t2- 


Under  spectrum  loadings,  da/dNref  and  da/dNenv  are  as.c  _d  to 
depend  on  the  effective  loading  parameters,  Kmax_eff  and  Kmin_eff, 
computed  by  the  Willenborg  model. 

The  following  assumptions  are  also  inherent  in  the  analysis: 

(a)  The  unloading  portion  of  the  stress  cycle  produces  no 
crack  growth  and  is  ignored  in  the  analysis.  The  stress 
intensity  range  and  stress  ratio  used  for  analysis  are 
determined  by  a  valley  stress  and  the  subsequent  peak 
stress. 

(b)  Any  half  cycle  having  a  stress  ratio  greater  than  0.9  is 
treated  as  a  sustained  load. 

(c)  In  the  computation  of  sustained  load  crack  growth,  a  cycle 
has  a  peak  stress  intensity  factor  of  Kmax-ef£»  as 
determined  by  the  Willenborg  model,  and  Reff=0. 

(d)  Fracture  is  assumed  to  occur  whenever  the  stress  intensity 
factor  exceeds  Kc. 
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Computation  of  cyclic  crack  growth  considers  only  the  loading 
portion  of  the  cycle,  i.e.,  S5  to  S6  is  a  load  cycle,  sg  to  S7  or 
S7  to  sg  are  assumed  to  produce  no  cyclic  growth.  Similarly, 
computation  of  sustained  load  crack  growth  uses  only  the  loading 
half  cycle  and  sustained  load,  the  growth  is  predicted  by 
Equations  (5)  through  (9).  The  cyclic  frequency  is  determined 
from  the  loading  period,  i.e.. 


f 


1 

2t4 


for  s3 


to  s4. 


(10) 


For  sustained  load,  i.e.,  half  cycles  having  R  >_  0.90, 


da  _  t  .  da 
dN  ~  C  dN 

env 


(ID 


where  for  sg  to  S7,  At  =  t7 . 

4.  WILLENBORG  MODEL 

a.  Generalized  Willenborg  Model  -  The  Willenborg,  et  al, 
model.  Reference  8,  as  generalized  by  Gallagher  and  Hughes 
(Reference  9)  is  the  load  interaction  model  extended  to  the 
analysis  of  environmental  effects.  The  original  Willenborg  Model 
was  developed  to  describe  crack  growth  retardation  following 
high-low  block  loadings.  It  is  based  on  observations  of  the 
following  phenomena: 


(a)  Retarded  crack  growth  occurs  whenever  the  maximum  applied 
stress  intensity  is  reduced. 

(b)  Such  retardation  is  directly  related  to  the  reduction  in 
maximum  stress  intensity. 

(c)  The  length  over  which  crack  growth  is  retarded,  i.e.,  the 
load  interaction  zone,  is  proportional  to  the  plastic  zone 
created  by  the  maximum  stress  intensity. 

(d)  There  is  no  retardation  of  growth  if  the  current  maximum 
stress  creates  a  load  interaction  zone  which  extends  out 
to  or  beyond  a  previously  established  interaction  zone. 

Based  on  these  observations,  Willenborg,  et  al,  assumed  that  the 
load  interaction  effects  were  caused  by  variations  in  local  stress 
intensity  as  the  crack  grows  through  the  residual  stress  field 
produced  by  the  overload(s). 


62 


NADC-79095-60 


Mathematically,  the  effective  stress  intensity  for  the 
Willenborg  model  is  expressed  in  Reference  8  as: 

Kef f  K  KRED 
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where  K  is  the  applied  stress  intensity,  Kred  Is  the  additional 
stress  intensity  required  to  extend  the  current  interaction  zone 
to  that  created  by  the  overload,  Aa  is  the  growth  following  the 
overload,  and  zql  is  the  overload  interaction  zone  size.  The 
effective  stress  intensity  range  and  stress  ratio  are  computed  as. 
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Thus,  as  noted  in  Reference  8,  the  Willenborg  model  predicts 
retardation  by  depressing  the  effective  stress  ratio  below  that 
remotely  applied  while  leaving  the  stress  intensity  range  intact. 
Since  Kred  decreases  as  the  crack  grows  through  the  overload 
interaction  zone,  the  Willenborg  model  predicts  that  maximum  retar¬ 
dation  will  occur  just  after  the  overload  and  that  the  growth  rate 
will  return  to  constant  amplitude  when  the  current  interaction 
zone  extends  to  the  end  of  the  overload  interaction  zone. 


Due  to  the  dependence  of  the  Willenborg  retardation  on  effec¬ 
tive  stress  ratio,  a  crack  growth  rate  relationship  which  inter¬ 
relates  the  influence  of  stress  ratio  with  stress  intensity  range 
must  be  used.  This  relationship  is  discussed  in  Paragraph  2. 
Overload  ratio  is  defined  as  the  ratio  of  overload  stress 
intensity  to  maximum  stress  intensity,  for  the  current  load,  and 
shut-off  ratio  as  the  overload  ratio  that  prevents  subsequent 
crack  growth. 

The  Willenborg  model  predicts  zero  value  for  the  maximum 
effective  stress  intensity,  and  thus  no  growth  should  gccur  when 
the  overload  ratio  is  two,  that  is,  when  Kred  equals  Kraax  in 
Equation  (12).  This  can  be  shown  by  rewriting  (12)  for  the 
maximum  effective  stress  intensity  as: 


Keff  =  K  -  [K0L  (1  -  ^A-j1/2  -  K" 
max  max  max  '  z„  '  max 

Uli 


(16) 
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where  is  the  maximum  overload  stress  intensity,  Aa  is  growth 

following  overload  and  zql  i-s  overload  interaction  zone  size. 

Immediately  following  the  overload,  Aa  is  usually  very  close  to 

zero  so  that  is  zero  when  K0L  is  twice  Km=y. 

max  max  max 

Test  results  obtained  by  several  investigators,  References  10 
through  13,  show  that  the  actual  crack  growth  shut-off  ratio  can 
be  greater  than  two.  Gallagher  and  Hughes,  Reference  9,  general¬ 
ized  the  Willenborg  model  to  correct  prediction  of  the  overload  to 
maximum  load  ratio  required  to  produce  cessation  of  crack  growth. 
They  proposed  modifying  Equation  (16)  so  that  for  R=0: 
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For  the  condition  of  no  growth,  a=0,  =  KmaxTfl  (threshold 

stress  intensity)  so  that,  max 
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where  the  shut-off  overload  ratio  (K0L  /Kmax)  must  be  obtained 

max 

from  tests  for  the  given  material,  thickness,  and  stress  ratio 
(underload  condition). 


Gallagher  and  Hughes  used  the  generalized  model  quite  success¬ 
fully  to  predict  the  number  of  cycles  required  to  return  to 
constant  amplitude  growth  rate  following  an  overload  in  two  steels 
having  different  yield  strengths.  Gallagher  and  Stalnaker, 
Reference  14,  also  used  the  generalized  model  to  predict  magnitude 
and  trends  of  crack  growth  rate  data  generated  under  transport¬ 
wing  simulation  loading.  The  correlation  of  test  and  analysis  was 
significantly  improved  over  that  of  the  original  model. 

The  computer  routine  for  the  generalized  Willenborg  Model  was 
extracted  from  the  Air  Force's  CRACKS  II  computer  program. 
Reference  15. 
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b.  Compression  Loads  Analysis  -  There  are  two  effects  of  com¬ 
pression  loading.  One  effect  is  to  increase  the  stress  intensity 
amplitude  of  the  cycle  following  compression;  analysis  of  this 
stress  ratio  effect  was  discussed  in  Paragraph  2.  The  second,  and 
frequently  greater  effect  is  to  accelerate  the  crack  growth  caused 
by  subsequent  load  cycles.  In  the  program  "Effects  of  Fighter 
Attack  Spectrum  on  Crack  Growth",  Reference  16,  modifications  were 
made  to  the  Willenborg  model  to  better  account  for  the  accelerated 
growth  following  compression.  The  approach  is  to  reduce  the 
overload  plastic  zone  size  based  on  the  minimum  stress  intensity 
applied  prior  to  the  overload: 
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where  Kmin  is  the  minimum  stress  intensity  applied  prior  to  the 
overload,  Fty  is  the  material  yield  strength,  and  Y  is  the  plane 
stress-plane  strain  coefficient;  Y  =  1  for  plane  stress  and 
=  1/2/2  for  plane  strain.  The  multiplier  3/32  was  empirically 
selected  to  correlate  predictions  with  the  compression  loads  test 
data  presented  in  Reference  17. 


The  improvement  in  prediction  accuracy  obtained  by  including 
compression  in  spectrum  loading  analysis  is  shown  in  Figure  50. 
The  test  data  was  obtained  from  the  program  "Effects  of  Fighter 
Attack  Spectrum  on  Crack  Growth",  and  in  this  series  of  tests  the 
ground  load  was  varied  to  investigate  the  effects  of  compression. 
The  earlier  Willenborg  model  predicted  small  impact  of  this 
variation,  the  improved  model  more  accurately  predicted  the 
acceleration  observed  in  test. 
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Figure  50.  Improved  Analysis  of  Compression  Effects  With  Residual  Stress  Intensity  Model 
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The  complete  environment-load  interaction  analysis  routine  is 
presented  in  the  form  of  a  user's  manual  in  Appendix  A. 

5.  MODEL  CALIBRATION  -  The  generalized  Willenborg  model  was  used 
to  account  for  the  effect  of  high  loads  in  a  spectrum  to  retard 
the  crack,  growth  produced  by  subsequent  low  load  cycles.  The 
Willenborg  model  accounts  for  this  retardation  by  using  two 
parameters;  the  overload  interaction  zone  size,  zql*  and  the 
shut-off  overload  ratio.  The  overload  interaction  zone  size 
determines  the  amount  of  crack  growth  affected  by  an  overload. 
Increasing  the  interaction  zone  size  retards  the  predicted  crack. 

Values  of  both  of  these  parameters  were  selected  by  correla¬ 
tion  with  results  obtained  with  tests  of  center  crack  panels.  The 
shut-off  overload  ratio  selection  was  based  on  results  of  single 
overload  tests  reported  in  Reference  1.  The  overload  interaction 
zone  size  was  selected  to  correlate  with  the  spectrum  tests  in  dry 
air.  The  premise  of  the  Willenborg  model  application  in  this 
program  was  that  once  parameter  values  were  found  which  correlated 
spectrum  analysis  and  test  results,  they  could  be  used  with  confi¬ 
dence  to  predict  the  results  of  similar  spectrum  tests.  The 
Willenborg  model  parameters  used  to  make  all  spectrum  crack  growth 
predictions  presented  herein  are  shown  in  Table  7. 


TABLE  7.  WILLENBORG  MODEL  PARAMETERS  USED  FOR  ANALYSIS 


Material 

R00T2 A 

OLMAX  A 

HP  9-4-0.30 

0.100 

3.50 

300 M 

0.028 

3.50 

Notes: 

A  R00T2  is  an  overload  interaction  zone  size  factor 
The  interaction  zone  size  is  defined  as 


OL ' 


R00T2 

27T 


A 


,0L 


max 

OLMAX  is  the  overload  ratio  - which  is  determined 


from  tests  to  shut-off  subsequent  crack  growth 


QPO3-O0M-36 


a.  Determination  of  Shut-off  Ratio  Using  Overloads  Tests  -  To 
determine  the  overload  ratio  which  shuts-off  subsequent  constant 
amplitude  crack  growth,  a  series  of  single  overload  tests  was 
performed  and  reported  in  Reference  1.  The  constant  amplitude 
crack  growth  between  overload  applications  was  sufficient  to 
preclude  interactions  of  the  overload  effects. 
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The  correlation  of  predicted  and  measured  delay  cycles  for 
overload  tests  in  both  300M  and  HP  9-4-. 30  were  used  to  select 
shut-off  ratio  of  3.5  for  both  steels.  The  correlation  of 
measured  and  predicted  delay  cycles  for  HP  9-4-. 30  tests  in  both 
air  and  salt  water  is  shown  in  Figure  51.  As  shown,  the  shut-off 
ratio  of  3.5  used  for  analysis  of  HP  9-4-. 30  steel  correlates  well 
with  the  higher  overload  ratios  but  predicts  much  less  delay  than 
measured  for  lower  overload  ratios.  This  should  result  in 
conservative  analysis  when  applied  to  spectrum  loadings. 


gpm-mm-w 


Figure  51.  Effect  of  Overload  Ratio  on  Delay  in  HP  0-4-0.30  Steel 
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b .  Determination  of  Overload  Interaction  Zone  Size  Using 
Spectrum  Test  Data  -  The  overload  interaction  zone  size  selected 
for  the  analyses  was  based  on  the  retardation  measured  in  spectrum 
tests  in  dry  air.  Once  the  shut-off  ratios  were  determined  from 
the  overload  test  results,  overload  interaction  zone  sizes  were 
determined  by  correlating  analysis  results  with  results  obtained 
from  the  center  crack  panel  spectrum  tested  in  dry  air  to  the 
accelerated  stress  history.  Results  of  these  correlations  are 
discussed  in  Section  VI. 

6.  SEMI-ELLIPTICAL  SURFACE  FLAW  ANALYSIS 


a.  Stress  Intensity  Solutions  -  Stress  intensity  solutions 
for  elliptic  surface  flaws  were  determined  using  a  slice  synthesis 
technique  described  in  Reference  1.  A  comparison  of  these  solu¬ 
tions  with  others  is  presented  in  Figure  52.  In  Figures  B-l  and 
B-2  of  Appendix  B,  the  solutions  are  graphically  summarized. 

Closed  form  expressions  are  presented  in  Appendix  B  which  can  be 
readily  used  in  crack  growth  prediction  computer  routines. 


The  slice  synthesis  technique  described  in  Appendix  B  was  used 
to  obtain  stress  intensity  solutions.  In  this  approach,  the  flaw 
is  idealized  as  a  series  of  center  crack  "slices"  joined  through 
the  thickness  with  edge  crack  slices  to  account  for  the  shear 
coupling  between  the  center  crack  slices.  The  results  presented 
in  Figure  52  show  good  agreement  with  those  obtained  by  Raju  and 
Newman  (References  18  and  19)  using  fine  mesh  finite  element 
models  which  are  the  most  detailed  employed  thus  far  in  the 
analysis  of  the  surface  flaw.  The  slice  synthesis  solutions 
include  the  effects  of  finite  width  as  well  as  thickness. 

b.  Prediction  of  Flaw  Shape  Change  -  The  elliptic  surface 
flaw  requires  analysis  of  the  growth  in  both  the  surface  and  depth 
directions.  Because  the  crack  aspect  ratio  changes  as  the  crack 
grows,  the  stress  intensity  relationships  at  the  surface  and  depth 
change.  It  would  be  unnecessarily  complicated  in  spectrum  analy¬ 
sis  to  track  the  growth  at  both  locations  on  a  cycle-by-cycle 
basis.  Instead,  the  changing  aspect  ratio  can  be  predicted,  based 
on  constant  amplitude  loading  analysis.  With  the  aspect  ratio 
known,  growth  at  only  one  location  needs  to  be  tracked,  simplify¬ 
ing  spectrum  crack  growth  prediction.  With  this  data 
pre-determined ,  it  is  then  possible  to  predict  stress  intensity 
relationships  as  a  function  of  surface  length,  permitting  the  two 
dimensional  growth  to  be  characterized  by  surface  growth  alone. 

Appendix  B  describes  a  routine  used  to  predict  growth,  shape 
changes,  and  stress  intensity  factor  corrections  for  a  semi¬ 
elliptical  surface  flaw  under  constant  amplitude  loading.  This 
routine  was  used  to  determine  flaw  shape  changes  and  stress 
intensity  factors  for  both  constant  amplitude  and  spectrum  crack 
growth  analyses.  For  spectrum  crack  growth  analyses,  the  highest 
stress  level  in  the  spectrum  and  initial  flaw  size  and  shape  were 
used  for  analysis. 
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Computation  of  stress  intensity  factors  at  the  crack  depth  and 
plate  surface  is  based  on  results  from  the  slice  synthesis  tech¬ 
nique  described  in  Reference  1.  Stress  intensity  factors  and 
crack  growth  both  at  flaw  depth  and  surface  are  computed  on  a 
cycle-by-cycle  basis.  The  following  assumptions  are  inherent  in 
this  routine: 

(1)  Constant  amplitude,  R=0,  sinusoidal  loading  is 
applied. 

(2)  Crack  growth  rates  at  plate  surface  and  flaw  depth 
are  dependent  on  the  K  which  is  applied  at  those 
points . 

(3)  Crack  growth  rate  versus  K  relationship  obtained  from 
a  center  cracked  panel  can  be  used  to  predict  flaw 
growth  at  any  point  along  the  crack  front. 


(4)  Stress  intensity  factor  expressions  are  assumed  valid 
until  flaw  depth  exceeds  90%  of  plate  thickness. 
Values  of  Ka  and  Kb  at  a/t=0.8  and  a/t=0.9  are  used 
for  extrapolations  beyond  a/t=0.9. 

(5)  Fracture  is  assumed  to  occur  when  the  stress 
intensity  factor  for  the  input  constant  amplitude 
stress  level  is  equal  to  Kc. 


c.  Correlation  of  Analysis  and  Test  -  Figures  41  through  44 
compare  measured  and  predicted  crack  growth  in  the  surface  flaw 


specimens  under  constant  amplitude  loading.  The  test  data  was 
obtained  from  the  constant  amplitude  elliptic  flaw  testing 
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SECTION  IV 

LANDING  GEAR  STRESS  HISTORY 


The  flight-by-flight  stress  history  used  in  the  verification 
test  phase  was  derived  for  the  axle-piston  fork  interface  area  on 
the  F-4J  main  landing  gear  system.  It  was  derived  from  the  design 
fatigue  load  spectrum  for  F-4J  landing  gears  and  includes 
estimated  times  of  stress  application.  The  times  required  to 
attain  the  various  landing  gear  loading  conditions  can  vary 
widely,  especially  during  braking  where  the  pilot  establishes  load 
duration.  The  times  computed  were  the  maximum  possible  for  each 
braking  condition. 

1.  F-4J  LANDING  GEAR  DESCRIPTION  -  The  F-4J  main  landing  gear 
(Figure  53)  is  a  retractable,  single  wheel,  pneudraulic  strut 
mounted  in  the  wing  outboard  of  the  fuselage  between  the  main  and 
rear  spars.  Retraction  of  the  gear  is  accomplished  by  a  hydraulic 
side  brace  actuator  having  one  end  attached  to  the  gear  and  the 
other  end  mounted  in  the  wing.  Retraction  of  the  actuator  causes 
the  gear  to  rotate  about  its  trunnion  axis  inboard  and  upward  into 
the  wing  cavity.  While  in  the  extended  position,  the  actuator  is 
converted  into  a  rigid  brace  by  means  of  an  internal  locking 
device.  Hydraulic  pressure  unlocks  the  device  prior  to 
retraction. 

Energy  developed  during  landing  is  dissipated  by  metering 
fluid  through  a  series  of  orifices  in  a  tube.  Vertical  loads 
encountered  during  taxiing  and  ground  handling  are  supported  by 
air  pressure.  Fore  and  aft  loads,  side  loads,  and  moments  on  the 
piston-fork  applied  at  the  wheel  ate  taken  out  as  a  couple  by  the 
upper  and  lower  strut  main  bearings.  Torsional  moments  about  the 
strut  centerline  are  reacted  by  the  torque  arm  lugs.  Vertical 
loads  are  reacted  by  internal  fluid  and/or  air  pressure  acting 
against  an  orifice  fitting. 

All  the  components  discussed  above  except  the  axle-piston  fork 
are  fabricated  from  4330V  steel  heat  treated  to  220-240  ksi 
ultimate  tensile  strength.  The  axle-piston  fork  is  an  integral 
unit  fabricated  from  300M  steel  heat  treated  to  280-300  ksi 
ultimate  tensile  strength. 

The  four  areas  of  the  F-4  main  gear  predicted  to  sustain  the 
most  fatigue  damage  are  identified  in  Figure  54.  The  axle-piston 
fork  interface  area  was  predicted  by  crack  growth  analyses  to  have 
the  most  severe  of  the  stress  spectra  considered.  The  axle-piston 
•  >rk  area  of  the  F-4  main  gear  is  the  only  area  known  to  have 
.  r  -d  iced  -racks  in  service. 

•  "'tin  landing  gear  was  selected  as  a  representative  highly 
■  •  *  if lgue  critical  component.  A  fatigue  spectrum  was  derived 

i  < 1  •  -piston  fork  region  based  on  design  loads.  This 
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Figure  54.  Areas  of  F-4J  Main  Gear  Predicted  to  Sustain  Greatest  Fatigue  Damage 


region  experiences  representative  loads  in  each  of  the  major 
loading  conditions  for  the  gear;  i.e.,  landing,,  braking,  turning 
etc.,  and  is  subject  to  most  of  the  lg  sustained  loads  when  the 
aircraft  is  parked. 
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2.  LOADS  SPECTRUM  SUMMARY  -  The  F-4J  main  landing  gear  were 
designed  to  the  ultimate  strength  requirements  of  MIL-A-8862  and 
the  repeated  loads  requirements  of  MIL-A-8866,  except  that  the 
service  life  was  conservatively  selected  as  4,000  takeoffs  and 
landings  rather  than  the  3,000  specified  in  MIL-A-8862.  A  scatter 
factor  of  two  was  used,  resulting  in  8,000  takeoffs  and  landings 
in  the  design  fatigue  spectrum. 

The  landing  gear  repeated  load  conditions  for  the  F-4J  main 
gear  are  presented  in  Table  8.  These  design  conditions  are  in 
accordance  with  MIL-A-8866.  Side  drift  landings  as  defined  in 
MIL-A-8863  were  included  although  not  specifically  required  by 
MIL-A-8866.  Landing  condition  19  was  separated  into  three  con¬ 
ditions  to  better  represent  the  highest  loads  in  the  spectrum. 
Landing  conditions  not  defined  in  Table  8  are  assumed  to  produce 
lg  stress  levels.  Turning  conditions  at  design  landing  gross 
weight  were  predicted  to  produce  negligible  damage  and  were 
omitted.  The  remaining  ground  handling  loads  conditions  were 
conservatively  assumed  to  occur  at  design  takeoff  gross  weight. 

Loads  for  both  ground  handling  and  landing  load  conditions  are 
summarized  in  Table  9.  The  sequence  in  which  loads  are  applied  to 
the  landing  gear  are  defined  in  Figure  55.  The  maximum  stresses 
for  the  axle-piston  fork  interface  area  for  each  loading  condition 
are  summarized  in  Table  10.  These  stresses  are  based  on  design 
values  and  not  on  actual  measurements  on  field  hardware. 

3.  FLIGHT-BY-FLIGHT  STRESS  HISTORY  DEVELOPMENT  -  The  event 
sequence  used  to  develop  the  flight  stress  history  for  the  axle- 
piston  fork  interface  is  summarized  in  Table  11.  This  stress 
sequence  is  used  for  every  flight  in  the  spectrum,  except  that 
pivoting  cycles  are  only  applied  every  tenth  flight  and  hard 
braking  cycles  before  landing  are  applied  every  second  flight. 

A  high  stress  level  landing  load  cycle  is  applied  every  sixth 
or  seventh  flight  in  the  spectrum.  Stresses  for  each  landing  load 
condition  were  repeated  uniformly  throughout  the  spectrum,  the  two 
highest  stress  levels  were  placed  at  the  middle  and  end  of  the 
spectrum.  The  resulting  landing  stress  level  distribution  is 
shown  in  Figure  56. 

A  comparison  of  the  flight-by-flight  stress  spectrum  with  the 
stresses  for  the  repeated  loads  spectrum  presented  in  Reference  1 
is  shown  in  Figure  57.  The  only  significant  differences  occur  at 
the  largest  stresses  where  the  two  highest  landing  stress  levels 
are  separately  identified  in  the  flight-by-flight  history.  The 
occurrences  of  each  loading  condition  in  the  reduced  design  spec¬ 
trum  and  the  flight-by-flight  spectrum  are  compared  in  Table  12. 
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TABLE  8.  F-4J  LANDING  GEAR  REPEATED  LOAD  CONDITIONS 


Description 

Condition 

Sink  Rate 
(ft/sec) 

Roll 

Angle 

Number  of  /^\ 
Cycles 

Parameters  Defining  Loads 

1 

15.6 

0°- 

330 

Landing  Weights  Vary  from  36,000  to 

2 

15.6 

2°-Lo  Gear 

205 

38,000  lb 

3 

15.6 

2°-Hi  Gear 

153 

4 

15.6 

5°-Lo  Gear 

65 

5 

15.6 

5°-Hi  Gear 

27 

6 

18.5 

0° 

122 

7 

18.5 

2°-Lo  Gear 

76 

8 

18.5 

2°-Hi  Gear 

53 

9 

18.5 

5°-Lo  Gear 

25 

10 

18.5 

5°-Hi  Gear 

15 

Arrested 

11 

20.4 

0° 

25 

Landing 

12 

20.4 

2°-Lo  Gear 

15 

13 

20.4 

2°-Hi  Gear 

13 

14 

20.4 

5°-Lo  Gear 

7 

15 

20.4 

5°-Hi  Gear 

5 

16 

23.4 

0° 

11 

17 

22.8 

2°-Lo  Gear 

9 

18 

22.4 

2°-Hi  Gear 

11 

19a 

22.8 

5°-Lo  Gear 

1 

14°  Pitch  Angle 

19b 

22.4 

5°-Lo  Gear 

1 

14°  Pitch  Angle 

19c 

22.8 

5°-Lo  Gear 

7 

10.5  Pitch  Angle 

20 

22.8 

5°-Hi  Gear 

8 

Side  Drift 

21 

- 

- 

10 

Side  Reaction  =  0.8  Vertical  Reaction 

Landing 

22 

— 

“ 

10 

Side  Reaction  =  0.6  Vertical  Reaction 

Ground  Handling  Load  Conditions 


Description 

Condition 

Aircraft 

Weight 

Number  of 
Cycles 

Parameters  Defining  Loads 

Turning 

23a,  23b 

Takeoff 

15,000 

Side  Reaction  =  0.4  Vertical  Reaction 

Hd  Braking 

24 

Takeoff 

12,800 

Drag  Reaction  =  0.8  Vertical  Reaction 

Md  Braking 

25 

Takeoff 

32,000 

Drag  Reaction  =  0.4  Vertical  Reaction 

Pivoting 

26 

Takeoff 

800 

Torque  =  0.5  Limit  Torque  Load 

Note:/l\ Number  of  cycles  represent  8,000  takeoffs  and  landings 

GPO3-0a3»-«0 
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TABLE  9.  F-4J  MAIN  LANDING  GEAR  DESIGN  FATIGUE  SPECTRUM 


Gear  Stroke 
in.  from 
Extended 


52.970 
63,750 
51,960 
56,940 
68,050 
75,130 
86,850 
60,700 

71.970 
71,590 
96,790 

101,420 

68,150 

82,950 

78,600 

119,760 

119,860 

85,420 

103,690 

98,300 

105,170 

89,050 


49,810 

49,240 


30,370 

18,170 


16.85  23,000  15,190 


-21,190 

-23,370 

-14,700 

-20,990 

-20,280 

-30,540 

-35,700 

-17,320 

-30,490 

-34,590 

-32,570 

-45,030 

-32,330 

-39,630 

-41,810 

-43,220 

-48,450 

-40,490 

-53,780 

-51,360 

-44,290 

-51,300 


-  5,000 
-16,000 
10,000 
-32,280 
-30,760 
-25,020 
10,000 


-40,020 

-29,670 


Medium  25  16.85  23,950  7,070 

Braking 


Pivoting  26  16.85  24,890  -  2,305 


Notes:  1  All  loads  are  in  Strut  Coordinate  System.  Py  =  vertical  load, 

PD  =  drag  load,  Pg  ■=  side  load,  Mp  =  moment  about  vertical 

load  axis 

2  Parking  load  at  takeoff  weight  is  Py  =  24,500  lb,  at  landing 
weight,  Py  =  15,470  lb 
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Arrested  Landing  Side  Drift  Landing 


Time  -» 


Figure  55.  Loading  Diagrams  for  F-4J  Main  Gear  Fatigue  Spectrum 
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TABLE  10.  STRESS  SPECTRUM  FOR  AXLE-PISTON  FORK  AREA 


Description 

Condition 

Stress 

(Max/Min) 

■  ■  . . 

Percent  of 
Maximum  Stress 

Number  of 
Cycles  a\ 

1 

93.8/0 

32.5/0 

330 

2 

116.6/0 

40.4/0 

205 

3 

69.9/0 

24.2/0 

153 

4 

94.3/0 

39.1/0 

65 

5 

65.2/0 

34.3/0 

27 

6 

136.7/0 

47.4/0 

122 

7 

169.1/0 

58.7/0 

76 

8 

7.14/0 

24.8/0 

53 

9 

169.2/0 

58.7/0 

25 

10 

108.0/0 

37.5/0 

15 

Arrested 

11 

164.5/0 

57.1/0 

25 

Landing 

12 

212.6/0 

73.7/0 

15 

13 

84.7/0 

29.4/0 

13 

14 

217.1/0 

75.3/0 

7 

15 

113.6/0 

39.4/0 

5 

16 

207.3/0 

71.9/0 

11 

17 

249.2/0 

86.4/0 

9 

18 

106.2/0 

36.8/0 

11 

19a 

288.3/0 

100.0/0 

1 

19b 

274.1/0 

95.1/0 

1 

19c 

258.9/0 

89.8/0 

7 

20 

118.9/0 

41.2/0 

8 

Side  Drift 

21 

216.7/0 

75.2/0 

10 

Landing 

22 

0.0/— 35 

0.0/- 12.1 

10 

Turning 

23 

8.85/— 8.1 

30.7/— 2.8 

7,50u 

Hd  Braking 

24 

41.1/36.6 

14.3/12.7 

Md  Braking 

25 

37.2/36.6 

12.9/12.7 

32,000 

Pivoting 

26 

37.2/36.6 

12.9/12.7 

800 

Note:  /l\  Number  of  cycles  represent  8,000  takeoffs  and  landings  i 
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TABLE  11.  FLIGHT  STRESS  SEQUENCE  FOR  AXLE-PISTON  FORK  AREA 


Description 

Condition 

Cycles  per 
8,000  Flights 

Cycles  per 
Flight 

Pivoting 

26 

800 

A 

Turning 

23 

8,000 

A 

Hard  Braking 

24 

4,000 

Medium  Braking 

25 

16,000 

2 

Landing 

1  -  22 

8,000 

A 

Hard  Braking 

24 

8,000 

i 

Medium  Braking 

25 

16,000 

2 

Exceedances  per  8,000  Landings 


NADC-7 
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Percent  Maximum  Stress 
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Figure  57.  Comparison  of  Flight-by-Flight  Spectrum  and  Design 
Spectrum  for  F-4J  Main  Landing  Gear 


.urn 


80 


NADC-7 909 5-60 


TABLE  12.  COMPARISON  OF  DESIGN  AND  FLIGHT-BY-FLIGHT  STRESS  SPECTRA 
F-4J  MAIN  LANDING  GEAR  AXLE-PISTON  FORK  SPECTRA 
8,000  TAKEOFF  AND  LANDINGS 


Description 

- 

-# 

Condition 

Percent  of 
Spectrum 
Maximum  Stress 

A\  Design 
Number 
of  Cycles 

/2\  Flight-by- Flight 
Number  of  Cycles 

1 

32.5/0 

330 

330 

2 

40.5/0 

205 

205 

3 

24.2/0 

153 

153 

4 

39.1/0 

65 

65 

5 

34.3/0 

27 

27 

6 

47.4/0 

122 

122 

7 

58.7/0 

76 

76 

8 

24.8/0 

53 

53 

9 

58.7/0 

25 

25 

10 

37.5/0 

15 

15 

Arrested 

11 

57.1/0 

25 

25 

Landings 

12 

73.7/0 

15 

15 

13 

29.4/0 

13 

13 

14 

75.3/0 

7 

7 

15 

39.4/0 

5 

5 

16 

71.9/0 

11 

11 

17 

86.4/0 

9 

9 

18 

36.8/0 

11 

11 

19a 

100.0/0 

9 

1 

19b 

95.1/0 

— 

1 

19c 

89.8/0 

— 

7 

20 

41.2/0 

8 

8 

Side  Drift 

21 

75.2/0 

10 

10 

Landings 

22 

0.0/- 12.1 

10 

10 

Turning 

23 

30.7  -  2.8 

7,500 

8,000 

Hd  Braking 

24 

14.3/12.7 

12,800 

12,000 

Md  Braking 

25 

12.9/12.7 

32,000 

32,000 

Pivoting 

26 

12.9/12.7 

800 

800 

Notes:  /j\  MCAIR  design  repeated  loads  spectrum 
/2\  Flight-by-flight  spectrum 


QPO 3-06 36-46 


4.  TIMES  FOR  STRESS  APPLICATION  -  The  spectrum  tests  include 
simulation  of  the  times  of  stress  application  as  well  as  the  value 
of  stresses.  Therefore,  an  evaluation  of  the  times  required  to 
attain  the  landing  gear  loading  conditions  was  performed.  Load 
application  times  vary  widely,  especially  during  braking  where  the 
pilot  dictates  the  duration  of  the  load.  The  times  for  stress 
application  computed  herein  are  the  maximum  possible  times  for 
each  braking  condition.  Two  stress-time  histories  were  created  - 
one  having  the  greatest  possible  duration  of  stress  application, 
the  other  having  the  same  stress  levels  applied  as  quickly  as 
possible  without  sustained  loads.  The  time  estimates  are  design 
estimates  and  not  based  on  actual  measurements  on  field  hardware. 
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Pivoting  -  Pivoting  was  assumed  to  occur  by  herd  braking  of 
one  main  gear,  simultaneously  turning  the  nose  gear  such  that  no 
drag  force  is  developed  on  the  nose  gear.  The  turning  torque  due 
to  the  locked  wheel  is 


T 


pVmglm 


where  u  is  the  coefficient  of  friction  (specified  to  be  0.8  for 
hard  braking ) . 

Vmg  is  the  main  gear  vertical  load  (computed  to  be  0.44  times 
the  aircraft  weight,  W)  and  is  the  distance  between  the  main 
gear  tire  and  the  aircraft  centerline  (8.94  feet). 

T  =  (0.8)(0.44W) (8.94)  =  3.15W,  ft-lb 


The  rotational  acceleration  of  the  aircraft  is 

5)  =  T/I 


where  w  is 
and  I  is 


the  yaw  rotational  rate 
the  yaw  moment  of  inertia 


where  g  is  gravitational  acceleration  (32.2  ft/sec2) 

and  p  is  the  yaw  radius  of  gyration  (computed  to  be  10.74  ft) 


«  <3-19W><32.21  .  0.879  rad/sec2 
(10.74)^ 


The  relationship  of  velocities  and  accelerations  are 

v  =  Rd 

a  =  Rui2 

hence  u  =  a/v 

where  a  *  F  /M 

s 

Fs  is  the  side  force  acting  through  the  aircraft  mass 
centroid  *  uW 

where  u  is  specified  to  be  0.4 
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M  =  W/g 

hence  a  =  yg 

and  v  is  the  velocity  (assumed  to  be  15  knots,  25.32  ft/sec) 
thus  a)  =  (0.4)  (32.2)/(25.32)  =  0.509  rad/sec. 

Assuming  that  angular  acceleration  increases  linearly  with  time  to 
the  maximum  angular  velocity,  the  time  required  to  reach  this 
velocity  is 

t  =  «/•  =  ( .509)/(0.879)  =  0.58  sec. 

Turning  -  Turning  was  assumed  to  occur  by  turning  the  nose 
wheel  with  no  braking  of  either  main  gear.  It  was  assumed  that 
the  aircraft  yaws  such  that  its  centerline  is  tangent  to  the  path 
of  motion  at  the  center  of  gravity.  It  was  assumed  that  a 
constant  velocity  turn  is  made  and  that  the  turning  radius  varies 
such  that  the  total  side  force  linearly  increases  and  decreases 
with  time  to  and  from  its  maximum  value.  The  turning  torque 
provided  by  the  nose  gear  is: 


T  "  ^NG  LN 


where  y  is  the  coefficient  of  friction  (specified  to  be  0.4  for 
side  force) 

V  _  is  the  nose  gear  vertical  load  (computed  to  be  0.12  times 
°the  aircraft  weight,  W) 

and  is  the  distance  between  the  nose  gear  tire  and  the 
centerline  of  the  main  gear  (20.4  ft). 

T  =  (0.4) (0.12W) (20.4)  =  0.979W  ft- lbs 

The  rotational  acceleration  of  the  aircraft  is 
<L  =  T/I 

where  w  is  the  yaw  rotation  rate  and  I  is  the  yaw  moment  of 
inertia  (3.58W)  ft- lb-sec2) 
u  =  (0.979W) (3.58W)  =  0.273  rad/sec2 
a)  =  a/v 

where  a  =  yg  =  (0.4) (32.2)  =  12.88  ft/sec2 
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and  v  =  is  the  velocity  (again  assumed  to  be  15  knots,  25.32 
f t/sec^ ) 

u  =  (12. 88) /2 5. 32  =  0.509  ft/sec. 

The  time  required  to  reach  this  maximum  angular  velocity  is 

t  =  0 j/a  =  0.509/0.273  =  1.86  sec. 

Hard  and  Medium  Braking  -  Braking  times  can  vary  from  very 
short  pulses  to  the  time  required  to  stop  the  aircraft  during  high 
speed  landing  or  a  rejected  takeoff.  Assuming  the  effect  of 
aerodynamic  drag  is  small  in  comparison  to  braking  forces,  the 
deceleration  is  computed  as: 


2pV 

MG 

x  =  — w —  g 

where  x  is  the  stopping  distance 

v  is  the  coefficient  of  friction  (specified  to  be  0.4  for 
medium  braking  and  0.8  for  hard  braking) 

VMG  is  the  main  gear  vertical  load  (0.44  times  the  air¬ 
craft  weight,  W) 

g  is  gravitational  acceleration 

W  is  aircraft  weight 

The  total  main  gear  load  is 

2V_  =  0.88W 

Li 

Hence  the  deceleration  rate  is: 

2 

x  =  0.88ug  ft/sec 

The  time  to  stop,  t,  depends  on  the  initial  velocity,  which  can  be 
as  high  as  185  knots,  (312.2  ft/sec)  for  an  aborted  takeoff. 

t  =  v/x  =  312. 2/(0. 88yg)  =  11.02/u  sec. 

Stopping  time  for  medium  braking  is  27.6  sec. 

Stopping  time  for  hard  braking  is  13.8  sec. 

The  time  to  apply  or  remove  braking  loads  was  assumed  to  be  one 
second . 
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Landing  -  The  time  required  to  reach  the  peak  vertical  load 
for  the  main  gear  was  found  from  drop  test  data  to  be  approxi¬ 
mately  0.18  sec.  This  value  is  used  for  all  landing  loads. 

A  sustained  stress  representing  a  lg  static  loading  at  take¬ 
off  gross  weight  was  arbitrarily  added  in  equal  0.33  second  inter¬ 
vals  between  loading  events.  The  resulting  stress-time  profiles 
used  for  load  conditions  are  summarized  in  Figure  58. 

5.  CYCLE-BY-CYCLE  STRESS  SPECTRA  -  The  stress-time  profiles  of 
Figure  58  were  assembled  according  to  the  event  sequence  of  Table 
11  to  produce  a  eye le-by- cycle  stress  history.  The  resulting 
stress-time  history  for  the  landing  gear  is  made  up  of  flight 
histories  of  various  lengths,  depending  on  which  events  are 
included.  The  stress  levels  and  times  for  load  application  for 
the  longest  flight  are  presented  in  Table  13.  Time  increments 
represent  the  time  required  to  obtain  the  given  load  from  the 
previous  load  level.  When  successive  loads  have  equal  magnitudes, 
the  time  increment  represents  sustained  load  duration. 

The  stress-time  history  for  the  longest  flight  is  shown  in 
Figure  59;  that  for  the  shortest  flight  is  shown  in  Figure  60. 

The  total  time  required  to  apply  8,000  flights  of  this  spectrum 
including  sustained  loads  is  about  338  hours. 
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TABLE  13.  SUMMARY  OF  STRESS-TIME  HISTORY  FOR  BASELINE  FLIGHT 


Stress  Number 

Delta  Time 

Percent  of 
Maximum  Stress 

Condition 

Description 

1 

0.33 

12.7 

26 

Take-off 

2 

0.58 

12.9 

26 

Pivoting 

3 

0.58 

12.7 

26 

- 

4 

0.33 

12.7 

23 

Turning 

5 

1.86 

30.7 

23 

6 

3.72 

-2.8 

23 

7 

1.86 

12.7 

23 

8 

0.33 

12.7 

24 

Hard 

9 

1.00 

14.3 

24 

Braking 

10 

13.80 

14.3 

24 

11 

1.00 

12.7 

24 

12 

0.33 

12.7 

25 

Medium 

13 

1.00 

12.9 

25 

Braking 

14 

27.60 

12.9 

25 

15 

1.00 

12.7 

25 

16 

0.33 

12.7 

25 

17 

1.00 

12.9 

25 

18 

27.60 

12.9 

25 

19 

1.00 

12.7 

25 

- 

20 

0.33 

12.7 

1-22 

Landing  /2\ 

21 

0.03 

0. 

1-22 

22 

0.18 

12.7 

1-22 

23 

0.18 

12.7 

1  -22 

24 

0.33 

12.7 

24 

Hard 

25 

1.00 

14.3 

24 

Braking 

26 

13.80 

14.3 

24 

27 

1.00 

12.7 

24 

28 

0.33 

12.7 

25 

Medium 

29 

1.00 

12.9 

25 

Braking 

30 

27.60 

12.9 

25 

31 

1.00 

12.7 

25 

32 

0.33 

12.7 

25 

33 

1.00 

12.9 

25 

34 

27.60 

12.9 

25 

35 

1.00 

12.7 

25 

Park 

Notes:  1  Stress  numbers  1  -  35  are  repeated  8,000  times  to  represent 

A  two  design  lifetimes 

/ 2\  Stress  level  22  is  replaced  by  a  higher  landing  stress 
in  every  sixth  or  seventh  flight 
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SECTION  V 

VERIFICATION  TEST  PLAN 


1.  TEST  PROGRAM  SUMMARY  -  The  verification  test  plan  is  sum¬ 
marized  in  Table  14.  This  test  program  was  performed  to  provide 
the  data  required  to  verify  the  crack  growth  prediction  methodol¬ 
ogy.  Spectrum  fatigue  tests  on  center  cracked  and  surface  flawed 
panels  were  performed  in  dry  air  and  synthetic  sea  water  using 
both  accelerated  and  sustained  load  stress  histories.  Five  tests 
in  each  steel  are  divided  into  three  series. 

o  Test  Series  1  -  One  surface  flawed  panel  in  the  dry  air 
environment  using  the  accelerated  stress  history  to 
evaluate  predictions  of  surface  flaw  growth  in  an  inert 
environment. 

o  Test  Series  2  -  One  center  cracked  panel  and  one  surface 
flawed  panel  in  sea  water  using  the  accelerated  stress 
history  to  evaluate  the  predictions  of  crack  growth  in  an 
aggressive  environment  minimizing  sustained  load  effects. 

o  Test  Series  3  -  One  center  cracked  panel  and  one  surface 
flawed  panel  in  sea  water  using  the  sustained  load  stress 
history  to  evaluate  predictions  of  crack  growth  under 
maximum  conditions  of  sustained  load  and  aggressive 
environment . 

TABLE  14.  VERIFICATION  TEST  PROGRAM 


Specimen 

Number 

Specimen  Type 

Environment 

Load  Frequency 
(cps) 

Objective 

1 

Surface  Flawed  Panel 

<10%  RH  Air 

'v  10 

Evaluate  crack  growth  in  an 
inert  environment 

2 

3 

Center  Cracked  Panel 
Surface  Flawed  Panel 

Alternate 
Immersion 
in  Sea  Water 

o  o 

r*  v— 

P  P 

Evaluate  crack  growth  in  an 
aggressive  environment 

4 

5 

Center  Cracked  Panel 
Surface  Flawed  Panel 

e 

o  p 
— » — » 

Evaluate  crack  growth  with 
sustained  load  in  an 
agressive  environment 

Te«t  series  is  identical  for  each  material  qp03-om»-«9 


2.  TEST  SPECIMENS  AND  TEST  CONDITIONS 

a.  Test  Specimens  -  The  center  cracked  panel  and  surface 
flawed  panel  specimens,  Figures  1  and  2,  were  used  in  this  test 
program.  The  reduced  sections  1.1  the  surface  flawed  panel  speci¬ 
mens  were  required  in  order  to  properly  simulate  the  stresses 
experienced  in  landing  gear  components  and  maintain  load  levels 
within  the  capacity  of  available  fatigue  test  equipment.  For  both 
steel  alloys,  there  was  minimal  specimen  width  effect  on  predicted 
critical  flaw  sizes,  or  in  the  predicted  crack  growth  rates  at  crack 
lengths  near  those  critical  flaw  sizes. 
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b.  Spectrum  Stress  Levels  -  The  spectrum  used  in  the  testing 
is  defined  in  Section  IV.  The  cycle-by-cycle  stress  history  was 
derived  for  the  axle-piston  fork  interface  of  the  P-4J  main 
landing  gear  system.  In  order  to  achieve  the  desired  life  (8,000 
take-offs  and  landings  in  the  dry  air  environment),  consistent 
stress  levels  had  to  be  selected.  The  maximum  spectrum  stress 
level  used  in  the  Air  Force  landing  gear  program  was  approximately 
85%  of  the  ultimate  strength  of  the  HP  9-4-. 30  material.  This  per 
centage  is  representative  of  the  maximum  stress  level  in  a  landing 
gear  component.  To  permit  direct  comparison  between  the  crack 
growth  rates  of  Air  Force  and  Navy  landing  gear  spectra,  as  well 
as  maintaining  direct  comparison  between  the  materials,  it  is 
desirable  to  use  the  same  maximum  stress  in  both  steel  materials. 
Therefore  185  ksi  was  selected  as  the  maximum  spectrum  stress 
level  for  surface  flawed  panels  in  both  steels. 

Stress  levels  and  precrack  lengths  for  the  center  cracked 
panels  were  selected  to  result  in  8,000  take-offs  and  landings  in 
the  dry  air  environment.  The  combination  of  precrack  length  and 
stress  levels  provided  initial  stress  intensity  factor  levels 
similar  to  those  of  the  surface  flawed  panels. 

c.  Initial  Flaw  Sizes  -  Initial  flaw  sizes  were  selected  to 
produce  the  maximum  possible  crack  growth  in  8,000  spectrum 
landings.  Predicted  crack  growth  under  the  Navy  landing  gear 
spectrum  (Section  IV)  was  very  small  because  of  the  small  magni¬ 
tude  of  frequently  applied  loads  and  the  infrequent  application  of 
the  highest  loads.  In  order  to  obtain  the  greatest  possible 
growth  in  8,000  spectrum  landings,  initial  flaw  sizes  were 
selected  which  result  in  failure  at  the  end  of  the  8,000  landing 
spectrum.  Because  of  the  difference  in  critical  stress  intensity 
factor,  Kc,  between  300M  and  HP  9-4-. 30  steel,  the  same  initial 
flaw  size  could  not  be  used  for  both  materials.  In  fact,  initial 
flaw  sizes  used  for  HP  9-4-. 30  steel  were  greater  than  critical 
flaw  sizes  in  300M  steel.  In  center  cracked  panels,  initial  flaw 
lengths  were  1.23  inches  in  300M  and  2.52  inches  in  HP  9-4-. 30 
steel. 

In  surface  flawed  panels  initial  flaw  geometries  were  devel¬ 
oped  from  small  semi-circular  EDM  notches,  0.05  inch  long  by  0.025 
inch  deep.  Semi-circular  flaw  geometries  were  used  because 
surface  flaws  under  uniform  tensile  stress  tend  to  grow  in  a 
roughly  semi-circular  shape.  Crack  growth  in  the  depth  is  accompa¬ 
nied  by  growth  at  the  surface,  allowing  surface  measurements  to 
characterize  the  overall  growth  of  the  flaw.  EDM  notch  and 
initial  flaw  geometries  used  for  verification  tests  are  shown  in 
Figure  61. 

3.  PRECRACKING  PROCEDURES  -  An  Electrical  Discharge  Machined 
(EDM)  notch  o£  size  and  shape  expected  to  produce  the  desired 
initial  flaw  was  introduced  on  the  specimen  surface.  Specimens 
were  then  precracked  to  obtain  the  desired  surface  length.  The 
target  EDM  and  precrack  sizes  and  shapes  are  indicated  in  Figure 
61.  The  precracking  stress  levels  were  100  ksi. 
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HP  9-4-0.30  Steel 


Dimensions  in  inches 


Figure  61.  Development  of  Initial  Flaws 


4.  TEST  PROCEDURES,  CRACK  GROWTH  MONITORING,  INSTRUMENTATION,  AND 
ENVIRONMENTAL  CONTROL  -  Cyclic  testing  was  performed  in  a 
Materials  Testing  System  (MTS)  test  system.  The  specimen  was 
loaded  through  self-aligning  hydraulic  grips.  The  load  spectrum 
application  was  controlled  through  a  mini-computer;  the  time- 
sequence  load  spectrum  was  defined  by  the  load  levels  and  the  time 
point  at  which  each  load  was  to  be  applied.  A  haversine  wave 
shape  was  used  between  consecutive  loads.  Sustained  load  was 
defined  by  two  consecutive  loads  of  equal  magnitude  applied  at  the 
beginning  and  end  of  the  sustained  load  period. 
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During  the  testing,  crack  lengths  were  optically  monitored 
using  a  linear  displacement  transducer  with  a  30X  microscope. 

Crack  growth  measurements  were  made  after  every  100  landings  for 
the  first  2,000  landings  and  after  every  500  landings  thereafter. 

The  chambers  used  to  maintain  environmental  conditions  are 
shown  in  Figures  3  and  4,  Pages  8  and  9.  The  alternate  immersion 
cycle  was  10  minutes  immersion  followed  by  50  minutes  of  blown  air 
drying.  Low  humidity  air  was  obtained  using  silica  jell  dessicant 
contained  in  a  sealed  mylar  case  enclosing  the  crack  (Figure  3). 
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SECTION  VI 

COMPARISON  OF  SPECTRUM  LIFE  PREDICTIONS  AND  TEST  RESULTS 


1.  SUMMARY  -  The  crack  growth  algorithms  described  in  Section  III 
were  calibrated  to  correlate  with  the  spectrum  crack  growth  data 
obtained  from  center  crack  panels  tested  in  dry  air.  The  cali¬ 
brated  model  was  then  used  to  predict  the  crack  growth  in  the 
center  crack  panels  tested  to  both  the  accelerated  and  real  time 
stress  histories  in  the  alternate  immersion  environment.  The  real 
time  stress  history  is  the  history  described  in  Section  IV  includ¬ 
ing  sustained  stress  times.  The  accelerated  stress  history  uses 
the  same  stress  levels  and  sequence  but  includes  no  sustained 
stress  periods.  In  the  accelerated  stress  history,  stress  levels 
were  applied  at  about  10  cps.  Predictions  were  based  on  cycle-by- 
cycle  analyses  of  the  stress  spectrum  of  the  F-4J  main  landing 
gear  axle-piston  fork  interface.  This  spectrum  involves  small 
peak  stresses  sustained  for  up  to  27  seconds,  as  well  as  high  peak 
stresses  rapidly  applied.  Good  correlation  was  obtained  between 
the  predicted  and  measured  growth  behavior  in  the  aggressive 
(alternate  immersion  in  sea  water)  environment. 

The  calibrated  model  was  also  used  to  predict  crack  growth  for 
surface  flawed  panel  tests.  Verification  tests  included  applica¬ 
tion  of  the  stress  history  both  in  real  time  and  at  accelerated 
rates  to  specimens  containing  surface  flaws.  Tests  were  conducted 
both  in  lab  air  and  in  alternate  immersion  in  sea  water.  Verifi¬ 
cation  test  results  confirm  the  ability  of  the  model  to  predict 
the  growth  of  both  center  cracks  and  semi-elliptical  surface  flaws 
under  spectrum  loading  both  in  inert  and  aggressive  environments. 
The  relative  severity  of  Navy  and  Air  Force  landing  gear  spectra 
on  crack  growth  was  investigated. 

2.  CENTER  CRACK  PANEL  SPECIMEN  RESULTS 


a.  Model  Calibration  -  Model  calibration  is  discussed  in 
greater  detail  in  Section  III.  The  load-environment  interaction 
algorithms  were  calibrated  through  correlation  with  constant 
amplitude  test  results  from  the  development  test  program.  The 
Willenborg  spectrum  crack  growth  analysis  parameters  were  deter¬ 
mined  from  single  overloads  tests  in  dry  air  and  sea  water  and 
through  correlation  with  the  results  of  tests  of  center  cracked 
panels  tested  to  the  accelerated  stress  history  in  the  dry  air 
environment.  The  load  interaction  zone  size  parameter  was  used  to 
provide  correlation  with  center  cracked  panel  results. 

b.  Spectrum  Crack  Growth  Analyses  And  Test  Results  -  Spectrum 
crack  growth  was  predicted  for  center  crack  panels  in  the  alter¬ 
nate  immersion  environment  under  the  accelerated  and  sustained 
load  stress  histories  using  the  calibrated  crack  growth  algor¬ 
ithms.  Table  15  summarizes  crack  growth  life  results  for  all 
verification  tests.  Figure  62  summarizes  crack  growth  life 
results  for  the  center  cracked  panels.  Results  of  analyses  and 
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tests  in  each  material  are  presented  in  Figures  63  and  64. 
Correlation  with  dry  air  test  results  is  good  because  the  model 
was  calibrated  to  correlate  with  those  results. 


TABLE  15.  SUMMARY  OF  SPECTRUM  TEST  LIVES 

Results  of  Verification  Tests 


Load 

Frequency 

(cps) 

Life-Landings  to  Failure 

Specimen  Type 

Environment 

300IU1 

A 

HP  9-4-0.30  A 

Test 

Analysis 

Test 

Analysis 

Center  Cracked  Panel 

<  10%  R  H  Air 

10 

12,100 

11,780 

14,500 

10,870 

Surface  Flawed  Panel 

10 

11,780 

11,780 

11,300 

11,780 

Center  Cracked  Panel 

Surface  Flawed  Panel 

Alternate 
Immersion 
in  Sea  Water 

10 

10 

3,870 

11,000 

11,780 

11,200 

8,500 

9,370 

10,870 

11,780 

Center  Cracked  Panel 

0.1 

3,500 

3,780 

4,380 

10,870 

Surface  Flawed  Panel 

0.1 

4,000 

3,780 

7,990 

11,780 

GP03-0838-82 

Notes: 


A  Initial  flaw  lengths  in  300M  steal  were  1.23  in.  for  center  cracked  panels  (48  ksi  maximum  stress)  and 
approximately  0.1145  in.  in  surface  flawed  panels  (185  ksi  maximum  stress) 

A  Initial  flaw  lengths  in  HP  9-4-0.30  steel  were  2.52  in.  for  center  cracked  panel  (48  ksi  maximum  stress) 
and  approximately  0.278  in.  in  surface  flawed  panels  (185  ksi  maximum  stress). 
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10  CPS  0.1  CPS  10  CPS  0.1  CPS 

GP03  0838  79 


Figure  62.  Comparison  of  Analysis  and  Verification  Test  Lives  for  Center  Cracked  Panels 
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Figure  63.  Verification  Test  Results  for  300M  Center  Cracked  Panels 
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Figure  64.  Verification  Test  Results  for 
HP  9-4-0.30  Center  Cracked  Panels 


Spectrum  crack  growth  predictions  for  center  cracked  panels  in 
the  alternate  immersion  environment  showed  inconsistent  correla¬ 
tion  with  test  behavior.  Generally,  good  correlation  was  obtained 
in  all  tests  prior  to  application  of  the  maximum  load  in  the  spec¬ 
trum.  Two  specimens,  one  300M  specimen  tested  to  the  accelerated 
time  history  and  one  HP  9-4-. 30  specimen  tested  to  the  real  time 
stress  history,  failed  during,  or  shortly  after,  application  of 
the  maximum  stress  level  in  the  spectrum  at  about  3,870  landings. 
Both  300M  and  HP  9-4-. 30  steel  center  cracked  panels  showed 
considerable  tunneling  in  alternate  immersion  tests,  increasing 
stress  intensity  factor  during  application  of  the  maximum  stress 
level.  Also,  both  materials  showed  some  variation  in  failure 
crack  lengths  in  the  constant  amplitude  tests.  Some  tunneling 
coupled  with  low  Kc  for  individual  specimens  probably  contributed 
to  the  early  failure  of  these  specimens. 
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The  300M  specimen  tested  to  the  real  time  stress  history 
showed  erratic  behavior.  Between  2500  and  3000  landings  crack 
growth  jumped  from  below  that  predicted  to  somewhat  above  that 
predicted.  The  crack  continued  to  grow  slowly  until  after 
application  of  the  maximum  stress  level  at  3,870  landings. 
Thereafter  the  crack  growth  stopped  (at  least  at  the  surface). 
After  8,000  landings  were  applied  under  real  time  loading,  the 
spectrum  was  accelerated.  The  alternate  immersion  cycle  was 
maintained  throughout  the  test.  The  crack  remained  essentially 
stationery  until  27,500  landings  had  been  applied.  At  that  point 
the  crack  began  growing  and  failure  occurred  at  28,380  landings. 

In  this  case  it  is  expected  that  the  increased  stress  intensity 
factor  produced  by  tunneling  may  produced  sufficient  plasticity  at 
the  crack  tip  to  halt  subsequent  growth.  Several  authors 
(References  9-11)  have  shown  that  overloads  can  arrest  a  growing 
crack  if  the  overload  is  of  sufficient  magnitude.  For  purposes  of 
the  comparisons  shown  herein  (Figure  62-64),  failure  of  this 
specimen  was  said  to  occur  at  the  application  of  the  maximum 
spectrum  stress  level  at  3,870  landings. 

3.  SURFACE  FLAW  SPECIMEN  TEST  RESULTS  -  Using  the  analysis  algor- 
ithms  described  in  Section  III,  spectrum  crack  growth  predictions 
were  made  for  the  surface  flaw  geometries. 

Predictions  of  the  shape  change  occurring  during  growth  from 
the  given  initial  flaw  geometries  are  obtained  from  constant 
amplitude  analysis,  as  described  in  Section  III.  At  any  surface 
crack  length,  the  crack  aspect  ratio  is  assumed  equal  to  that  ob¬ 
tained  by  constant  amplitude  analysis. 

Table  15  and  Figure  65  summarize  the  predicted  and  measured 
crack  growth  lives  for  each  surface  flaw  test.  The  crack  growth 
lives  shown  in  Figure  65  are  presented  for  growth  to  fracture  from 

the  largest  initial  flaw  length  in  each  material  type.  In  300M 
steels,  this  flaw  length  was  0.114  inch  and  in  HP-9-4-.30,  it  was 
0.278  inch. 

Comparisons  of  predicted  and  measured  crack  growth  for  each 
surface  flaw  test  are  presented  in  Figures  66  and  67.  Correlation 
of  crack  growth  predictions  with  the  surface  flaw  test  results  in 
30 0M  as  shown  in  Figure  66  is  poor  except  for  the  alternate  immer¬ 
sion  test  using  the  10  cps  spectrum.  While  correlation  with  crack 
growth  appears  poor,  correlation  with  test  lives  is  good  as  indi¬ 
cated  in  Table  15  and  Figure  65.  This  discrepancy  was  explained 
by  examination  of  the  fracture  surfaces  of  each  specimen.  These 
examinations  showed  that  for  the  dry  air  test  and  the  slow  cycle 
rate  alternate  immersion  test,  crack  growth  at  the  surface  (where 
measurements  were  taken)  was  arrested,  while  growth  beneath  the 
surface  continued.  Crack  sizes  at  failure,  beneath  the  surface, 
approached  those  predicted.  Unfortunately  corrosion  of  the  crack 
surfaces  during  storage  made  it  impossible  to  determine  subsurface 
crack  lengths  accurately. 
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Figure  65.  Comparison  of  Analysis  and  Verification  Test  Lives  for  Surface  Flawed  Panels 
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Figure  66.  Verification  Test  Results  for 
300 M  Surface  Flawed  Panels 
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Figure  67.  Verification  Test  Results  for  HP  9-4-0.30  Surface  Flawed  Panels 


4.  CRACK  GROWTH  UNDER  NAVY  LANDING  GEAR  SPECTRUM  -  This  test 
series  verified  the  slow  crack  growth  rates  predicted  for  cracks 
growing  under  the  Navy  landing  gear  spectrum.  As  shown  in  Figure 
68,  crack  lengths  used  in  verification  tests  were  significantly 
larger  than  those  found  in  the  initial  flaw  survey  of  References  1 
and  2.  Yet  test  lives,  for  both  center  cracked  and  surface  flawed 
panels,  under  the  Navy  landing  gear  spectrum  are  longer  than  the 
8,000  landings  desired  in  the  verification  tests. 
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Figure  68.  Predicted  Growth  of  Small  Surface  Flaws  Under  Navy  Landing  Gear  Spectrum 


The  Navy  landing  gear  spectrum  is  characterized  by  a  few  occur¬ 
rences  of  high  stress  levels  caused  by  high  sink  rate  (carrier) 
landings  and  many  occurrences  of  very  low  stress  levels  caused  by 
ground  handling  loadings  (taxiing,  turning,  pivoting,  and  brak¬ 
ing).  As  shown  in  Figure  69,  all  of  the  high  peak  stresses  in  the 
Navy  gear  spectrum  are  caused  by  carrier  landing  loads.  In  fact, 
the  Navy  landing  gear  spectrum  without  carrier  landing  stresses  is 
similar  to  Air  Force  landing  gear  spectrum  used  in  Reference  1,  as 
shown  in  Figure  70.  In  this  comparison,  both  stress  spectra  have 
been  normalized  to  the  parking  load  stress  level  so  that  differ¬ 
ences  in  aircraft  weights  and  design  criteria  are  minimized. 
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Figure  69.  Navy  Main  Landing  Gear  Stress  Spectra 


Ratio  of  Stress  Level  to  1  g  Stress  Level 


GP03-10S2-6 


Figure  70.  Comparison  of  Navy  Landing  Gear  Spectrum  Without 
Carrier  Landings,  and  Air  Force  Landing  Gear  Spectrum 
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When  carrier  landing  stress  levels  are  included  in  the  Navy 
gear  spectrum,  as  they  were  in  this  program,  the  Navy  stress 
spectrum  differs  significantly  from  the  Air  Force  gear  spectrum 
used  in  the  program  of  Reference  1.  Exceedance  curves  for  these 
spectra  are  compared  in  Figure  71.  Ground  handling  stress  levels 
in  the  Navy  spectrum  are  no  more  than  35  percent  of  the  maximum 
stress  level  in  the  spectrum,  while  ground  handling  stress  levels 
in  the  Air  Force  spectrum  exceed  70  percent  of  the  maximum  stress 
level.  This  difference  in  relative  stress  level  due  to  ground 
handling  loads  makes  the  Navy  spectrum  much  less  severe,  from  a 
crack  growth  life  point  of  view,  than  the  Air  Force  spectrum. 

This  difference  in  severity  is  reflected  in  Table  16  which 
compares  crack  growth  lives  for  surface  flawed  panel  tests  using 
the  Air  Force  and  Navy  spectra. 


Percent  Maximum  Stress 
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Figure  71.  Comparison  of  Air  Force  and  Navy  Main  Landing  Gear  Stress  Spectra 
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TABLE  16.  COMPARISON  OF  TEST  LIVES  UNDER  AIR  FORCE  AND  NAVY 
LANDING  GEAR  SPECTRA 


Steel 

Material 

Environment  A 

Sustained 

Loads 

Life-Landings  /$\ 

Air  Force 
Spectrum 

Navy 

Spectrum 

Dry  Air 

No 

60 

m 

300M 

Aggressive 

No 

20 

mm 

Aggressive 

Yes 

12 

EZI 

Dry  Air 

No 

MM 

HP  9-4-0.30 

Aggressive 

No 

wM 

Aggressive 

Yes 

500 

Notes:  QP0M#3M9 

A  Aggressive  environment  for  the  Air  Force  spectrum  is  continuous  immersion  in 

3.5%  NaCI  salt  water,  for  Navy  spectrum  is  alternate  immersion  in  synthetic  sea  water. 

/^\  Initial  flaw  lengths  are  0.1 14  in.  in  300M  steel  and  0.278  in.  in  HP  9-4-0.30. 

Maximum  stress  level  for  both  spectra  is  185  ksi. 


While  the  implication  of  this  data  is  that  cracks  will  not 
grow  in  landing  gear  of  Navy  aircraft,  service  experience  shows 
that  this  is  not  the  case.  Corrosion  cracking  is  a  primary  cause 
for  gear  replacement  in  the  field.  These  flaws  are  caused  by  high 
residual  tensile  stresses  in  local  regions  or  by  actuator  loadings 
or  other  loadings  independent  of  the  primary  loading  of  the  gear 
during  landings.  To  minimize  crack  growth  in  Navy  gear,  residual 
stresses,  actuator  loading,  or  redundant  loads  may  require  more 
attention  during  design  than  stresses  caused  by  primary  loading. 
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SECTION  VII 
CONCLUSIONS 


Immersion  in  sea  water  accelerates  crack  growth  in  300M  steel. 
The  acceleration  decays  with  drying  time  after  immersion,  indicat¬ 
ing  that  continuous  immersion  is  a  more  aggressive  environment  for 
crack  growth  in  300M  steel  than  is  alternate  immersion. 

Immersion  in  sea  water  has  little  impact  on  crack  growth  rate 
in  HP-9-4-.30  steel  except  at  stress  intensity  factor  ranges  below 
30  ksi  «'/in.  Alternate  immersion  at  low  cycle  rates  tends  to 
increase  threshold  stress  intensity  factor  from  below  7  ksi  >^Tn  at 
10  cps  to  25  ksi  /"in  at  0.1  cps.  When  the  maximum  stress 
intensity  factor  is  below  these  threshold  values,  crack  growth 
arrests  during  the  drying  portion  of  the  immersion  cycle  and  does 
not  reinitiate  during  subsequent  immersions. 

Currently,  we  conjecture  that  the  drying  cycle  of  the  alter¬ 
nate  immersion  environment  served  to  dry  corrosion  products  on  the 
crack  faces,  increasing  crack  closure  to  the  point  that  the  crack 
could  not  grow  in  the  dry  air  environment.  Subsequent  ten  minute 
immersion  times  were  not  sufficient  to  dissolve  the  product  such 
that  the  crack  could  grow.  Once  the  specimen  was  precracked 
beyond  the  apparent  threshold,  the  crack  growth  rates  closely 
approached  those  found  in  the  dry  air  tests. 

Crack  growth  in  alternate  immersion  environments,  such  as  that 
applied  in  this  program,  can  be  accurately  predicted  using  linear 
summation  of  sustained  load  crack  growth  rates  derived  from  cyclic 
tests  in  an  aggressive  environment  and  cyclic  crack  growth  rates 
determined  from  tests  in  an  inert  environment.  The  decay  in  crack 
growth  rate  following  removal  from  an  aggressive  environment  is 
accurately  predicted  by  a  function  which  decays  exponentially  with 
time  after  removal. 

The  main  landing  gear  stress  history  used  in  this  program  is 
based  on  Navy  specified  sink  rates  for  carrier-based  fighter 
aircraft  and  so  is  typical  of  that  for  landing  gear  of  many 
carrier-based  fighter  aircraft.  This  stress  history  is  character¬ 
ized  by  few  applications  of  high  stress  levels  caused  by  high  sink 
rate  landings  and  numerous  applications  of  low  stress  levels 
produced  by  ground  handling  (taxiing,  turning,  pivoting,  braking) 
load  cycles.  This  stress  history  is  typical  of  many  landing  gear 
components.  However,  nose  gear,  actuators,  and  some  redundant 
components  can  be  subjected  to  significantly  different  stress 
histories. 
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This  stress  history  produces  very  little  crack  growth  over  the 
aircraft  service  liftime,  except  for  crack  lengths  near  critical 
flaw  size.  Large  flaw  sizes  were  used  in  the  verification  tests 
reported  herein.  So  little  growth  occurs  under  sustained  loads  in 
this  stress  history  that  crack  growth  lives  are  not  as  sensitive 
to  environment  as  were  those  reported  in  previous  work  with  Air 
Force  fighter  aircraft  landing  gear  (Reference  1). 
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APPENDIX  A 

CRACK  GROWTH  ANALYSIS  ROUTINE  INCLUDING 
ENVIRONMENT-LOAD  INTERACTION  EFFECTS 


1.  INPUT  DEFINITIONS  -  The  computer  program  has  the  following 
inputs : 

Analysis  inputs  include 

a.  title  identifying  problem  (TITLE) 

b.  threshold  stress  intensity  factor,  critical  stress 
intensity  factor  for  material  thickness  (DELKTH,  KSUBC) 

c.  initial  crack  length  (AZERO) 

d.  stress  ratio  correction  parameter,  plastic  zone  size 
correction  factor,  and  shut-off  overload  ratio  (ALPHA, 
ROOT,  OLMAX) 

e.  stress  intensity  correction  factor  type  (ICOR) 

f.  for  single  and  double  cracks  from  holes  (ICOR  =  1  or  2, 
respectively),  radius  of  the  hole  and  shortest  distance 
perpendicular  to  the  load  from  hole  center  to  the  plate 
edge  (RADIUS,  ECCEN)  are  input;  for  constant  shape 
semi-elliptical  surface  flaws  (ICOR  =  3)  the  flaw  aspect 
ratio,  plate  thickness,  and  shortest  distance  from  crack 
center  to  plate  edge,  (A,  C,  THICK,  ECCEN)  are  input;  and 
for  a  through  crack  (ICOR  =4),  the  shortest  distance  from 
the  crack  center  to  plate  edge  (ECCEN)  is  input. 

g.  if  a  stress  intensity  factor  correction  table  is  to  be 
used,  the  number  of  points  in  the  table  (NPTS),  and  the 
crack  length  and  K  correction  factors  for  each  point 
(TABLEA ( I ) ,  TABLEB ( I ) ) ,  are  input. 

h.  number  of  spectrum  repetitions  to  be  analyzed  (NBLKS) 

i.  design  limit  stress  (DLS) 

3.  number  of  flights  between  printouts  of  accumulated  time, 
crack  length,  and  reference  K  (PRTFLTS) 

k.  index  for  consideration  of  environmental  effects,  0  if  no 
load-environment  interaction  is  to  be  considered,  1  if 
load-environment  interaction  will  be  considered. 
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Material  data  inputs  include  one  data  set  if  load-environment 
interaction  is  not  considered,  two  data  sets  if  load-environment 
interaction  is  considered.  The  two  data  sets  describe  high 
frequency,  sinusoidal  wave  crack  growth  rate  data  from  tests  in  an 
inert  environment  and  low  frequency,  trapezoidal  wave  data  from 
tests  in  an  aggressive  environment.  The  single  data  set  used  for 
anaysis  without  load-environment  interaction  is  the  same  as  the 
first  of  the  following  sets. 

High  frequency,  sinusoidal  wave  data  (used  with  or  without 
load-environment  interaction  effects) 

a.  title  identifying  material,  product  form,  and  loading 
(TITLE) 

b.  monotonic  tensile  yield  stress,  cyclic  tensile  yield 
stress  (TYLD,  CYLD) 

c.  number  of  points  in  da/dN  table  and  sine  wave  frequency 
(NDADNA,  FA) 

d.  AK  and  da/dN  values  of  points  in  table  (DKA(I),  DADNA(I)). 

Low  frequency,  trapezoidal  wave  data  (used  only  for 
load-environment  interaction  effects) 

a.  title  identifying  material,  product  form,  and  loading 
(TITLE  2) 

b.  monotonic  tensile  yield  stress,  cyclic  tensile  yield 
stress  (TYLD,  CYLD) 

c.  number  of  points  in  da/dN  table  and  trapezoidal  wave 
frequency  (NDADNB,  FB) 

d.  AK  and  da/dN  values  of  points  in  table  (DKB(I),  DADNB ( I ) ) . 

Spectrum  inputs  include: 

a.  title  identifying  spectrum  (TITLE) 

b.  number  of  stress  levels  in  spectrum  and  number  of  flights 
(landings)  represented  by  spectrum  (NLYR,  FLTS ) 

c.  stress  level  and  time  increment  for  each  half  cycle 
(SMAX(I),  TIME  (I)). 

2.  DESCRIPTION  OF  OUTPUT  -  The  program  outputs  include  titles  of 
analysis,  material,  and  spectrum  input  data  sets,  design  limit 
stress  (DLS),  initial  flaw  length  (AZERO),  and  a  table  of  elapsed 
flights  (landings),  crack  length,  and  reference  stress  intensity 
factor  printed  at  the  interval  specified  by  PRTFLTS.  When  Kmax 
exceeds  Kc,  fracture  is  predicted  to  occur  and  the  flight  in  the 
spectrum,  crack  length  at  fracture,  and  Kmax  is  printed. 
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LISTING  OF  PROGRAM 
c 

C  THIS  PROGRAM  COMPUTES  SPECTRUM  CRACK  GROWTH 

C  INCLUDING  ACCELERATION  DUE  TO  SUSTAINED  LOADINC 

C  IN  AN  AGGRESSIVE  ENVIRONMENT.  THE  UII.LENBORG 

C  MODEL  IS  USED  TO  PREDICT  SPECTRUM  LOAD 

C  INTERACTION  EFFECTS.  A  LINEAR  SUMMATION  APPROACH 

C  IS  USED  TO  PREDICT  ENVIRONMENTAL  ACCELERATION. 

CYCLE-BY-CYCLE  ANALYSIS  IS  USED.  PROCRAM  WAS 
DEVELOPED  BY  C.R.SAFF  OF  MCDONNELL  AIRCRAFT  CO., 

C  ST.  LOUIS,  MISSOURI,  TELEPHONE  (314)  232-3356. 

C 

DIMENSION  S MAX (1000) .TIME (1000) ,CARRAY(100) ,SNARAY(100) 
DIMENSION  S  LAND ( t  50  0) 

INTEGER  RETARD, TITLE (20 ) 

REAL  KM AX , KM  IN , KS UB C , KM XE F F , KMNEF F 

COMMON/ A/ A, RADIUS , ECCEN , AOC, THICK, NPTS , TAB LEA ( 5 0 ) , TABLEB ( 5 0 ) 
COMMON/ B/DKC ( 32 )  , DKD ( 3  2 )  ,  DADN  D ( 3  2 ) , DADT (3  2)  ,NDADT , NEN, DELKTH 
PI-3.  141  59265 
RETARD-0 
C 

C  READ  AND  ECHO  INPUT  DATA 

C 

C  ANALYSIS  INPUTS  INCLUDE- 

C 

C  TITLE  WHICH  IDENTIFIES  PROBLEM 

C  DELKTH-THRESHOLD  STRESS  INTENSITY  FACTOR 

C  KSUBC-CRITICAL  ST°ESS  INTENSITY  FACTOR  FOR  MATERIAL  THICKNESS 

C  AZERO-INITIAL  CRACK  LENGTH 

C  RC-STRF.SS  RATIO  CORRECTION 

C  ROOT  2 -PL A  ST I C  ZONE  SIZE  FACTOR 

C  OLMAX-SHUT-OFF  OVERLOAD  RATIO 

C  ICOR-STRF.SS  INTENSITY  FACTOR  CORRECTION  TYPE 

C  1-SINGLF.  CRACK  FROM  HOLE 

C  2 -DOUBLE  CRACK  FROM  HOLE 

C  3-SEMI-ELLIPTICAL  SURFACE  FLAW  (FIXED  A/C) 

C  A -THRO UGH  CRACK 

C  RADIUS-RADIUS  OF  CRACKED  HOLE 

C  F.CCEN-SHORTEST  DISTANCE  PERPENDICULAR  TO  LOAD 

C  FROM  HOLE  CENTER  TO  PLATE  EDGE 

C  AOC-RATIO  OF  SURFACE  CRACK  DEPTH  TO  HALF  LENGTH 

C  THICK-PLATE  THICKNESS 

C  ECCEN-SHORTEST  DISTANCE  FROM  CRACK  CENTER  TO  PLATE  EDGE 

C  NPTS-NO.  OF  POINTS  IN  STRESS  INTENSITY  FACTOR.  CORRECTION  TABLE 

C  TABLE A(I ) ,TARLEB(I )-CRACK  LENGTH  AND  K-CORRECTION 

C  FACTOR  VALUE  FOP.  TABLE 

C  NBLKS-NO.  OF  SPECTRUM  REPETITIONS  TO  BE  ANALYZED 

C  DLS-DES IGN  LIMIT  STRESS 

C  PRT  FLTS-NO .  OF  FLIGHTS  BETWEEN  PRINTOUTS  OF  ACCUMULATED 

C  TIME, CRACK  LENGTH, AND  REFERENCE  K 

C  NEN-INDF.X  FOR  CONSIDERATION  OF  ENVIRONMENTAL  EFFECTS 

C  0-SINCLF  DA/DN  CURVE  IS  TO  BE  USED 

C  NO  FREQUENCY  EFFECT  CONSIDERED 

C  1-ENVIRON  MENTAL  EFFECT  WILL  B  F.  CONSIDERED 

C  INPUT  HIGH  FREO.  SINF.  WAVE  DATA  & 

C  LOW  FREQ.  TRAPEZOIDAL  WAVE  DATA 

C  ALFA-RATE  OF  DECAY  OF  ENVIRONMENTAL  ACCELERATION 

C  FREF-REFERF. NCF.  FREOUENCY  FOR  ENVIRONMENTAL  ACCELERATION 

C  P-EXPOtlF.NT  OF  FREQUENCY  RATIO  USED  TO  DETERMINE 

C  ENVIRONMENTAL  ACCELERATION 

C 
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READ(1,2200)  TITLE 
READ (1,1 000 )  DELKTH.KSUBC 
READ (1,1 000 )  AZERO 
REA D  ( 1  ,  1000 )  RC , ROOT  2 , OLMAX 
IF(RC.LT.-1.85)  RC  — 1-85 
REA D(l, 2 000)  ICOR 
GO  TO  (10,  10,20, 30)  .ICOR 
10  REA  D ( 1  ,  1 000 )  RADIUS, ECCEN 
GO  TO  40 

20  READ(l.lOOO)  AOC  ,  Tl!  IC  K  ,  ECC  EN 
GO  TO  40 

30  READ(1  ,  1000)  ECCEN 
40  R EA 0(  1,2000)  NPTS 

IF(NPTS.LE.O)  GO  TO  60 
DO  50  1*1, NPTS 

50  READ(l.lOOO)  TA fi LEA ( I ) , TAB LEB { I ) 

60  RCAD(1 ,2000)  N3LKS 
READ ( 1 , 1000)  DLS 
REA  P ( 1  ,  1000)  PRTFLTS 
READ ( 1,2000)  WEN 

IF(NEN.NE.O)  READ(l.lOOO)  ALP.FREF.P 
WRITE(6, 2300)  TITLE 
C 

C  READ  MATERIAL  GROWTH  RATE  DATA 

C 

C  COMPUTE  DA/DT  AND  DA/DN 

C 

CALI,  ENVIR  (CARR  AY,  SNA  RAY,  NDADN  ,  KS  L'BC  ,  TYLD  ,  CYLD  ) 

DO  70  I-l.NPADN 

SNA  RAY ( I )  ■ ALOG 10 (SNARAY ( I )  * (KSUBC-CAR  RAY ( I ) ) /KSUBC ) 

70  CARRAY(I)“AI.0C10  (  C  ARRAY  (  I  )  ) 

C 

C  READ  SPECTRUM  INPUT  DATA 

C 

C  SPECTRUM  DATA  INPUTS  1NCLUDE- 

C 

C  TITLE  WHICH  IDENTIFIES  SPECTRUM  INPUT 

C  HLYR-NO.  OF  STRESS  LEVELS  IN  SPECTRUM 

C  FLTS-NO.  OF  FI.  IGHTS  (LANDINGS)  REPRESENTED  BY  SPECTRUM 

C  SM AX ( I ) "STR  ESS  LEVEL  FOR  EACH  HALF  CYCLE 

C  TIME(l)-TIME  INCREMENT  FOR  EACH  HALF  CYCLE 

C 

100  READ(4,2200)  TITLE 

READ (4 ,20  00)  HLYR,FI.TS 

DO  110  I-l.NLYK 

READ (4 , 1000)  SMAX(I)  ,TIME(I  ) 

SMAX(I )«SMAX(I ) *DLS/ 1 00. 

110  CONTINUE 

R EAD (5 ,  5000  )  (S LAND ( I ) , I- 1  ,  1  204 ) 

5000  F ORMAT ( 1  OF 6 • 1) 

URITE(6, 2300)  TITLE 
WRITE(6, 3460)  NBLKS.DLS 
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SET  UP  FOR  CYCLIC  GROWTH  ANALYSIS 

7.K-0 . 3  304  5+0.  1  5  1  64  *  RC -0 . 0  l  4  7  6*  RC  *  *2 
FR-CYLO/TYLD-l . 

ZK-ZK*( 1 . +0 . 6*FR-0 . 1 56*FR**2 ) 

SIGMAY-CYLD 

A-AZ  ERO 

ASUBP-AZERO 

CALL  BETA (ICOR, BETAT ) 

WRITE (6,4300)  A 
PFLTS-O. 

WRITF.(6,3500) 

THESE  LINES  ARE  USED  TO  SET  UP  COUNTERS  FOR  ENVIRONMENT 
AND  RANDOM.  LANDINC  LOADS 


SMIN-O. 

NCYC-0 

KEY-1 

TIMEl-O. 

TYME-0. 

TYII 1-3600  . 

TYM 2-600. 

LSTRS-0 
KONT  3A-6 
J3A-2568 
J3B-8000 
KONT 1-0 
KONT 2-0 
KONT 3-0 

DO  230  Jl-l.NBLKS 
KONT 1-KONT 1+1 
KONT2-KONT2+1 
KoNT  3-KO  NT  3  +  1 
IF(K0NT1.0T.  10)  KO  NT  1*1 
IF( KONT1.LT. 10)  KEY  1-0 
IF (KONT l.EQ. 10)  KEY  1 » 1 
I F ( KO  NT2.GT.2)  KONT 2-1 
I F (  KONT  2 .  LT  .  2  )  KF.Y2-0 
IFCKONT2.EO. 2)  KEY2-1 
I F ( J 1 . GT .  J3A)  KO  NT  3A - 7 
IF(J  I.LF..J3B)  GO  TO  1  IS 
KONT 3 A -6 
J3A-J  3A+80O0 
J3B-J  38+8000 

119  I F ( KO  NT  3 • OT . KO  NT  3 A  )  KONT 3-1 
IFCKONT 3. LT . KONT 3A  )  KEY3-0 
I F  (  KO  NT  3 .  F.  0  •  KO  NT  3  A  )  KF.Y3-1 
DO  220  J4-1.NLYR 

TYME-T  YME+T IME ( J  4 ) 
IF(KFY.EO.l)  GO  TO  120 
IF(TYME.LT.TYMl)  GO  TO  125 
KEY-1 

TYM 1 -TYM 1+3 600. 

TYII  2 -TYM  2  +  360  0 . 

GO  TO  125 

120  I  F (TYilE.  LT  .  TYM2  )  CO  TO  125 
KF.Y-0 

125  CONTINUE 
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CYCLIC  CROWTH  ANALYSIS 
DADN-0. 

SIGMAX-SMAX(J4) 

IF  (  J  4 .  LT  .  4 .  AND  •  KE  Y  1 .  EQ .  0  )  GO  TO  220 
IF(J4.LT. 8. OR. J4.GT.il)  GO  TO  126 
IF ( KEY2 • EO. 0 )  GO  TO  220 
126  IF( J 4. NE . 22 . OR • KEY  3. NE •  1 )  GO  TO  129 
LSTRS-LSTRS+1 
IF(LSTRS.CT. 1204)  LSTRS-1 
SIGMAX-S LAND (LSTRS)*DLS/ 100. 

COMPUTE  KM AX  &  KM  IN 


129  CONTINUE 

130  I F(SICMAX. LT. SMIN )  GO  TO  210 

KMIN-SMIN*SQRT(PI*A)*BETAT 

R-SMIN/SICMAX 

IF(KMIN.LT.SKMIN)  S  KM  IN -KM  I N 
KMAX«S I CM A X ‘SORT (PI *A) *  BET AT 
RETARD  -  1 

IF(A.GE.ASUBP)  GO  TO  160 
COMPUTE  EFFECTIVE  KMAX  6  KM  IN 

SIGREF-(SIGMAY*SQRT(2.*( ASUB  P-A )/A))/BETAT/SQRT(ROOT2) 

SIGRED-SICREF-SIGMAX 

IF(OLMAX.NE.O.)  GO  TO  140 

PHI-1 . 

GO  TO  150 

140  THRSLD-DELKTH*(  1  .  -P  “2  ) 

PH  I- ( 1 .-THRSLD/XMAX) /(OLMAX-1. ) 

SIGR  ED-PH I ‘S IGRED 
150  IF(SICRED.LT.O. )  SIGRED-O. 

SIGMAX-S ICMAX-S IGRED 
IF(SIGMAX. LE . 0. )  GO  TO  210 
KMNEFF-KM  IN-SIGREO*SQRT(PI *A)*BETAT 
IF(KMNEFF.LT .0. )  KMNEFF-0. 

KMXEFF-SIGMAX*SQRT(PI*A) *  BETAT 
RF.FF-KMNEFF/KMXEF  F 

COMPARE  EFFECTIVE  PLASTIC  ZONE  INTERFACE  WITH  PREVIOUS 
INTERFACE 

RSUBYl«(KMAX“2-3./32.*SKMIN“2)/2  .  /P  I  ‘ROOT  2  /  SIGMA  Y“2 
IF(A+RSUBY1.LT.  ASUBP)  GO  TO  170 

THIS  CYCLE  IS  NOT  RETARDED 


160  RETARD-0 

RSUBY-(KMAX“2-3./32.*SKMI 

asubp-a+rsuby 
S  KM  IN  - 1  .  E  20 
KMXEFF-KMAX 
KMNEFF-KM  IN 
REFF-KMIN/KMAX 


N  “  2  )  /  2  •  /PI  *ROOT2/SIGMAY“2 


COMPUTE  DA/DN 


1  70 


I  F  ( KM  XF.FF.LE.  1 . 00 1  ‘KMNEFF  )  GO  TO  t  RO 
IF(KMAX.CE.KSUBC)  GO  TO  240 
THRSLD-DELKTH*(1  .-R“2  ) 

DELTAK-KMXEFF* ( 1 . -REFF) 

DELK-O  . -REFF)  *(  l  .-( 1  • -RF.FF  )  *ZK) /(  1  . -ZK) ‘KMXEFF 

IF ( REFF • LT . 0. )  DE LK- ( 1 . -ZK*EXP ( 0 . 1 *REF F) ) ‘KMXEFF/ ( 1 . -ZK 
DELTAK-ALOG lO(DELK) 

CALL  TLU<C ARRAY, SNARA Y , ND A DN , DELT A K , DADN) 
DADN-10.“DADN‘KSUBC/  (KSUBC-KMAX) 
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COMPUTE  ENVIRONMENTAL  ACCELERATION 

180  IF(NEN.EQ.O)  GO  TO  210 

IF(KMXEFF.LT.DELKTH)  GO  TO  210 
IF(TIME(J4)  .  EQ . 0 . )  GO  TO  210 
F-0.5/TIME(J4) 

XKMX-ALOGIO(KMXEFF) 

IF (R • GT .0.9)  GO  TO  190 

CALL  TLU (DKD , DADND, NDADT , XKMX, DADTL) 

GO  TO  200 

190  CALL  TLU(DKC, DADT , NDADT , XKMX, DADTL) 

F“l . /TIME ( J  A ) 

200  IF(DADTL.LT.-99.  )  DADTL  — 99. 

DADT l-l 0. **DADTL 

IF ( KE  Y . NE . 0 )  GO  TO  205 

T1-TIME1 

T2-T  IMF.  1+TIM  E  (  J  4  ) 

DADT  I -DADT  1*  (EXP  (-AI.P*T  1  )-EXP  (-ALP*T2)  )  /  ALP/(T2-T  1  ) 

GO  TO  206 
20  5  TIME  1-0. 

206  EDADN -DADT1/F*(FREF/F)**P 
DA DM  »D ADN+E  DADN 

INCREMENT  A  AND  COMPUTE  STRESS  INTENSITY  CORRECTION 

210  A-A+DADN 

TIME1-TIME  1+T  IM F.  (  J  4  ) 

SMIN-SMAX(J4) 

IF (DADN • LE . 0. )  CO  TO  220 
IF(RDADN . LT. DADN )  R  DADN -DADN 
CALL  BETA( ICOR, BETAT) 

REFK-DLS  *SQRT (P1*A)*BET,\T 
IF (RETARD. EQ.O)  AS UB P-A S UB P+D ADN 
220  CONTINBE 

PRINT  RESULTS 

CFLTS-FLTS  *J  1 

IF (CFLTS-PFLTS . LT . PRTFLTS)  GO  TO  230 
PFLTS-PFLTS+rRTFLTS 
WRITE  (6  ,  4400)  C  FLTS  ,  A  ,  P  F.F  K 
230  CONTINUE 
GO  TO  250 
240  WRITE (6, 41 00) 

CFLTS-FLTS*J 1 

WRITE(6 , 4200)  CFLTS,A,REFK 
250  STOP  1776 
1000  FORMAT ( 5F  1 0. 0 ) 

2000  F OR M AT (I  10, 5F 10.0) 

2200  F ORMAT (2  0A4 ) 

2300  F0RMAT(/1X, 20A4) 

3460  FORMAT( /T2, 18, *  BLOCKS  AT  '.E12.5,*  PSI  DESIGN  LIMIT  STRESS') 
3500  FORMAT  ( '  FLIGHTS  A  RF.FF'/) 

4100  FORMAT (/6 8(1H») /I X,' KMAX  EXCEEDS  KSUBC.  PROBLEM 

♦'TERMINATED' /68(1H*) /, IX, 'LAST  CALCULATED  VALUES  ARE'/) 

4200  F0RMAT(1X, 'FLIGHT  IN  SPECTRUM  '.F16.2/ 

+  IX, 'CRACK  LENGTH  ', E16.fi/ 

+  IX, 'REFK  '  ,E16. 8) 

4300  F  ORMAT  ( /  / '  BEGIN  SPF.CTRUM  CRACK  GROWTH  ANALYSIS,  A  -*,F10.5/) 
4400  FORMAT(F8.0,F10.5,F10.5) 

4500  F  ORMAT ( / ) 

END 
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SUBROUTINE  TI.U  (X  ,  Y  ,  N  ,  XVA  L  ,  Y  VAL  ) 

* 

TABLE  LOOK  UP  ROUTINE 

DIMENSION  X(N),Y(N) 

IF(X(N) .CE.XVAL)  GO  TO  10 
I-N 

CO  TO  30 
10  DO  20  I«1,N 

IF ( X ( I  )  .GE.XVAL)  CO  TO  30 
20  CONTINUE 
30  IF(I.EQ.l)  1-2 

YVAL- (Y ( I )-Y (I-1))/(X(I)-X(I-1) )*(XVAL-X(I-1  ) )+Y (1-1  ) 

RETURN 

END 


SUBROUTINE  B  ETA  (  ICOR  ,  B  ETAT  ) 

STRESS  INTENSITY  FACTOR  CORRECTIONS 

COMMON /A/ A, RADIUS , ECC E N , AOC , TH IC K , NPTS , TAB LE A ( 50 ) , TAB LEB ( 5 0 ) 
DATA  PI /3  .  141  5926  5/ 

CO  TO  (10,20, 30,40) , ICOR 

SINGLE  THR0UCH  CRACK  FROM  OPEN  HOLE 

10  A A  1 -1 . /COS (PI  * ( A+2 .* RADIUS )/2 . / ( 2 . * ECC EN-A  )  ) 

IF(AA1.LE.0. )  AA1-1.E20 

BOWIE-Q. 67620624(0.8733015/(0. 324 544 2+A / RADIUS ) ) 
BETAT-B0WIE*SQRT(AA1) 

CO  TO  50 

DOUBLE  THROUGH  CRACK  FROM  OPEN  HOLE 

20  AAl-1. /COS(PI*(A+RADIUS )/2./ECCEN) 

IF(AA1 .LE. 0. )  AAI-1.E20 

BOWIE-0. 94 3851 0+(0. 6805078/(0. 27 7196 5+A/ RADIUS) ) 

BETAT -B  OW IE  *S  QRT ( A A  1  ) 

GO  TO  50 

NEWMAN'S  SEMI-ELLIPTICAL  SURFACE  FLAW  CORRECTION 

30  AT-A/THICK 

IF ( AT . GT  .0.95)  GO  TO  35 

XMW-1 . /S QRT (COS ( P I  * A/ 2 . / AOC / EC C EN  ) ) 

XM1-1. 1 3-0 . 09*AOC 
Q-l  .  4-1 . 4 7*  AOC**  1.64 

XME-XM1+S QRT (Q / AOC -XM l)*AT**(2.+8.*  AOC  *  *3 ) 
BETAT-XMW*XME/SQRT(Q) 

GO  TO  50 

35  A- (S  QRT ( AOC )+ 1 . ) **  2  *  B  E  TAT  **  2  *  A  /  4 . 

ICOR -4 

CENTER  CRACK  -FINITE  WIDTH 


40  AAl-1 ./COS (PI*A/2 ./ECCEN) 

IF ( AA1 . LT . 0. )  AAl-1. E20 
BETAT-SQRT(AAl) 

TABULAR  CORRECTION 

50  I F  (  N  PTS  .  EO .  0  )  GO  TO  60 

CALL  TLU (TAB LEA, TAB LEB, NPTS, A, CORl) 
BETAT-BETAT*COR 1 
60  RETURN 
END 
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SUB ROUT  INF  ENVIR ( D KA , DADNA , NDADN A , KS UBC , T YL D , CYL D ) 

COMPUTATION  OF  SUSTAINED  LOAD  GROWTH  DATA 
FOR  AGGRESSIVE  ENVIRONMENT 

DIMENSION  AK(32 ) ,DADNA(50) ,DADN3(50) ,DKA(50) ,DKB (50) 

INTEGER  TITLE1 (20) ,TITLE2(20) 

REAL  KSUBC 

COMMON /B/DEC( 32  ) ,DKD (32)  ,DADND(32) , DADT (32), NDADT , NEN , DE LKT  11 
DATA  AK/l*|2«|3*|4t|3*|6>|7i|8<|9i|10*)12«|I4»|16> .18*t20.)22.5» 
42  5.  , 30.  , 35. ,  AO.,  A 5. .  50. ,60. , 70. ,80. ,90. ,  10  0. ,  1  10. ,  120. , i 30- , 
4140.  ,  1  50. / 


READ  DA/DN  DATA  FOR  10  CPS  ,  R-0,  IN  DRY  AIR 
INPUT  DATA  INCLUDES  - 

TITLF.l  WHICH  IDENTIFIES  MATERIAL  4  PRODUCT  FORM  FOR  REFERENCE 
TYLD “MONO TONIC  TENSILE  YIELD  STRESS 
CYLD-CYCL1C  TENSILE  YIELD  STRESS 

NDADN  A-NO .  OF  POINTS  IN  DA/DN  TABLE  (SINE  WAVE) 

FA-LOAD  CYCLE  FREQUENCY 

DKA (1 )-D  ELT  A  K  VALUES  OF  POINTS  IN  DA/DN  TABLE 
DADNA(I )-DA/DN  VALUES  OF  POINTS  IN  DA/DN  TABLE 


REA  D  (  2  ,  1000)  TITLF.l 
1000  FORMAT ( 2  OA  4) 

R  EAD (2,  3000)  TYLD.CYLD 
REA  D  (  2  ,  2  000 )  NDADNA , FA 
2000  FORMAT (I  10, 5F 10. 0 ) 

NR  -0 

DO  10  I - 1 , NDADNA 

READ(2  ,  3000)  D  KA ( I )  ,  DADN A ( I  ) 

3000  FORMAT (8F10.4) 

IF(0KA( I ) .GE. KSUBC  )  GO  TO  5 
GO  TO  10 
5  NR-NR+1 
10  CONTINUE 

NDADN A -NDADN A -NR 
IF(NEN.NE.O)  GO  TO  15 
12  WRITE  (6,  1050)  TITLE  1 
1050  FORMAT(/2X,2044) 

GO  TO  999 

15  WP.ITE(6,  1  100)  TITLE  1 

1100  F0RMAT(/1X,'  DA/DT  IS  DETERMINED  FROM  '/1X.20A4) 
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c 

C  READ  DA/DN  DATA  FOR  0.1  CPS,  P-0,  IN  AGGRESSIVE  ENVIRONMENT 

C 

C  INPUT  DATA  INCLUDES  - 

C 

C  TITLE2  WHICH  IDENTIFIES  MATERIAL  &  PRODUCT  FORM  FOR  REFERENCE 

C  NDADNB-NO.  OF  POINTS  IN  DA/DN  TABLE  (TRAPEZOIDAL  WAVE) 

C  FB-LOAD  CYCLE  FREQUENCY 

C  DKB ( I )-DELTA  K  VALUES  OF  POINTS  IN  DA/DN  TABLE 

C  DADNB(I )-DA/DN  VALUES  OF  POINTS  IN  DA/DN  TABLE 

C 
C 

READ (3 , 1000 )  TITLE2 

RE AD  (  3 , 2  000)  NDADNB.FB 

DO  20  I-l.NDADNB 

READ ( 3 , 3000)  DKB (I ) ,DADNB(I ) 

DKB (I  )-ALOG10(DK3 (I )  ) 

20  DADNB(I  )-ALOG10(DADNB(I )) 

WR IT E ( 6 ,  1200)  TITLE2 
1200  FORMAT ('  AND ' / 1 X , 20A 4 ) 

DO  25  I-l.NDADNA 
DKA ( I )-AL0G10(DKA( I ) ) 

25  D ADN A (I ) - ALOG 10(DADNA(I)) 

DETERMINE  DA/DT 

DKLIM-ALOG 10(1 . 01 *DELKTH) 

UPLIM-ALOGIO(KSUBC) 

DO  40  J-l , 32 
DO  30  1-1,32 
B  K-AI.OC  10(AK(I)) 

IF(BK.LT.ONLIM)  CO  TO  30 
I F ( BK.OT . UPLIM)  GO  TO  45 
CALL  TLU(DKA,DADNA,NDADNA,BK,DADN1) 

CALL  TLU(DKB,DADNB, ND ADN B , B K , DA DN 2 ) 

DADN 1-10. **DADN 1 
DADN  2-10.**DADN2 
IF(DADN2.LE.DADN1)  GO  TO  30 
GO  TO  35 
30  CONTINUE 

35  DADT (J )- (DADN2-DADN 1  )*FB 
DKC ( J ) -B  K 

DN  LIM-ALOG10(AK(I  +  1 ) ) 

40  CONTINUE 

45  N DADT “J-l 
DO  50  I - 1 , NDADT 
DADT (I )-AL0C 10(DADT(I ) ) 

50  CONTINUE 

I F (NDADT . CT . 0 )  GO  TO  46 
NEN-0 
CO  TO  999 

INTEGRATE  DA/DT  FOR  R-0,  1  CPS  SINE  WAVE 
C 

46  DO  130  1-1, NDADT 
DKX-10.**DKC(I) 

CALL  GAUSS(DKX.CC) 

DAD!!D(I)-ALOC10(CC) 

DKD (I )-DKC (I  ) 

130  CONTINUE 

DO  150  I • 1 , NDADN A 
DKA(I)-10.**DKA(I) 

DA DM A(I )-10.**DADNA(I  ) 

150  CONTINUE 
999  RETURN 
END 
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SUBROUTINE  0AUSS(X4,CC) 

10  POINT  CAUSS'  RULE  INTEGRATION  ROUTINE 
DIMENSION  U(5),R(5) 

COMMON / B/ DKC ( 32 )  ,DKD(32) ,  DADN D ( 3 2 ) , DADT ( 3 2 ) , ND A DT , NE N , DEL KTH 
DATA  U/0.0744371695,0.216697697,0.339704784,0.432531683, 
60.486953264/ 

DATA  R/0. 147762 11 2, 0.13463336, 0.1 09 5431 81, 0.074 72 56746, 
60.0333356722/ 

DATA  PI/3.14159265/ 

A“0 . 

B-0. 55 

C1-0.5MB+A) 

C2-B-A 

S-O. 

DO  20  1-1,5 
W-C2*U  (t ) 

X-Cl+W 

J-l 

10  XK-X4/2.*(l.-COS(2.*PI*X>) 

XKL-ALOG 1 0 (XK ) 

IF(XKL.CT .0. )  GO  TO  12 

Y-l.E-99 

GO  TO  14 

12  CALL  TLU(DKC,DADT,NDADT,XKL,YL) 

Y«10.**YL 
14  S-S4R(I)*Y 

IF ( J . EQ. 2 )  GO  TO  20 

X-Cl-W 

J-2 

CO  TO  10 
20  CONTINUE 
CC-S*C2 
RETURN 
END 
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SAMPLE  INPUT 


ANALYSIS  FILE 


SURFACE  FLAWED  PANEL  ANALYSIS 


15.5 

60. 

0.05725 

0. 

0.028 

4 

0.5 

7 

.  0545 

.  6696 

.0614 

.  683  6 

.  070  5 

.  6945 

.083 

.  705 

.  0998 

•  7172 

.  1243 

.  72  92 

.1457 

.  72  95 

100000 

135  . 
500. 

1 

0.00175 

0. 1 

HIGH  FREQUENCY  DA/DN  DATA  FILE 


DA/DN  DATA  -  300M  STEEL  -  10  CPS  -  DRY  AIR 
248.  228. 

13  10. 


12. 

.  63E-6 

14. 

.  90E-6 

16. 

1.  31F.-6 

13. 

1. 85E-5 

20. 

2. 55E-6 

22.5 

3. 6E-6 

25. 

4. 75E-6 

30. 

7.5E-6 

35. 

12. 8E-6 

40. 

18.2E-6 

45. 

25. 5E-5 

50. 

40.0E-6 

60. 

600. E-6 
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LOW  FREQUENCY  DA/DN  DATA  FILE 

DA/DM  PAT  \  -  300»1  STEEL  -  0.1  CPS  -  ALTERNATE  IMMERSION  IN  SEA  WATER 


10 

0. 1 

18. 

200. E-6 

20. 

330. E-6 

22.  5 

400. E-6 

25. 

430 . E-6 

30. 

475. E-6 

35. 

525. E-6 

40. 

5  70  . E-6 

45. 

600. F - 6 

50. 

650. E-6 

60. 

710. C-6 

SPECTRUM 

FILE 

F-4J  MAIN  G EAR  AXLE-PISTON  FORK  STRESS  HISTORY  -  FAST 


35 

1 . 

12.  7 

0.001 

12.9 

0.05 

12.7 

0.05 

12.7 

0.001 

30.  7 

0.  05 

0. 23 

0.05 

12.  7 

0.  05 

12.  7 

0.  001 

14.3 

0.05 

14.3 

0.001 

12.7 

0.05 

12.7 

0.001 

12.9 

0.  05 

12.9 

0.001 

12.7 

0.05 

12.7 

0.  001 

12.9 

0.05 

12.9 

0.001 

12.7 

0.05 

12.7 

0. 001 

0. 

0.05 

12.7 

0.05 

12.  7 

0.05 

12.  7 

0.001 

14.3 

0.05 

14.3 

0.001 

12.7 

0.05 

12.7 

0.001 

12.9 

0.05 

12.9 

0.001 

12.7 

0.05 

12.7 

0.001 

12.9 

0.05 

12.9 

0.001 

12.  7 

0.05 
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RANDOM  LANDING  STRESS  LEVELS 


32.  5 

32.5 

32.5 

32.5 

40.4 

3  2.5 

24.2 

32.5 

47.4 

32.5 

40.4 

32.5 

32.5 

24  .  2 

58. 7 

32.5 

40.4 

39.1 

4  7.4 

32  .  5 

40.4 

24.3 

24.2 

32.  5 

32.5 

32  .  5 

32.  5 

40.4 

47.4 

2  4.2 

58. 7 

32.5 

32.  5 

32.5 

40.4 

3  2.5 

39.1 

24.2 

47.4 

32.5 

40.4 

40.4 

24.2 

34  .  3 

24.8 

58.7 

5  7.1 

58.  7 

47.4 

32.5 

32.  5 

4  0.4 

32.5 

24.  2 

39.  1 

32.5 

3  2.5 

4  0.4 

4  7.4 

32.5 

32.  5 

24.2 

58.  7 

40.  4 

3  2.5 

32.  5 

4  0.4 

2  4.8 

47.4 

24.2 

32.5 

32.5 

32.  5 

39.  1 

40.4 

24.2 

47.4 

5  8.7 

37.  5 

73  .  7 

32.  5 

40.4 

32.5 

32.5 

40.4 

2  4.2 

40.4 

47.4 

34.3 

24.8 

39.  1 

29.4 

24.2 

58.  7 

5  7.1 

58.7 

32  .  5 

47.4 

40.4 

32.5 

32.  5 

24.2 

32.5 

40.4 

40.4 

2  4.2 

4  7.4 

36.8 

71.9 

58.7 

39.  1 

32.5 

24.8 

32.5 

32.5 

4  0.4 

24.2 

47.4 

-12.1 

75  .2 

32.  5 

32.5 

40.4 

32.5 

24.2 

58.  7 

40.4 

47.4 

39.1 

32.5 

24.2 

34.3 

86.4 

32.5 

40.4 

2  4.8 

32.  5 

4  7.4 

32.5 

40.4 

CM 

• 

<NI 

58.  7 

5  7.1 

58.  7 

32.  5 

4  0.4 

47.4 

39.1 

24.2 

41.2 

32.5 

40.4 

32.5 

40.4 

2  4.2 

47.4 

5  3.7 

24 . 8 

37.5 

73  .  7 

32.5 

32.5 

32.5 

40.4 

2  4.2 

39.  1 

47.4 

32.5 

32.5 

4  0.4 

39.8 

75.3 

24.2 

*58.7 

32.5 

40.4 

47.4 

34.  3 

32.5 

24.2 

24.8 

40.4 

32.5 

39.  1 

29.4 

40.4 

47.4 

24.2 

32.5 

58.7 

57.  1 

58.  7 

40 . 4 

32.5 

32.5 

24.2 

47.4 

32.5 

40.4 

32.5 

40.4 

24.2 

39.  1 

24.8 

58.7 

32.5 

47.4 

32.5 

32.5 

32.5 

40.4 

24.2 

32.5 

32.5 

40.4 

47.4 

36.8 

71.9 

24.2 

58.  7 

39.  1 

34.  3 

40.4 

40.4 

40.4 

47.4 

24.8 

24.2 

40.4 

2  4.2 

47.4 

58.  7 

39.1 

57.1 

58.7 

37.  5 

73.7 

-12.  1 

75.2 

39.4 

32.5 

32.  5 

24.2 

32.5 

40.4 

47.4 

32.5 

32.5 

24.8 

32.5 

24.2 

40.  4 

58.  7 

32.5 

32.5 

47.4 

40.4 

24.2 

39.1 

32.5 

32.5 

32.5 

4  0.4 

24.2 

47.4 

34.3 

86.4 

32.5 

58.7 

40.4 

32.5 

24.8 

40.4 

24.2 

47.4 

39.1 

29.4 

32.5 

32.5 

32.5 

40.4 

32.5 

24.2 

4  0.4 

58.  7 

47.4 

57.  1 

58  .  7 

32.5 

32.5 

24.2 

32.  5 

4  0.4 

47.4 

24.8 

39.  1 

32.5 

4  0.4 

2  4.2 

58.  7 

41.2 

32.5 

32.5 

40.4 

47.4 

2  4.2 

32.5 

32.5 

32.5 

32.5 

40.4 

34.  3 

24.2 

39.  1 

47.4 

58.7 

40.4 

24.8 

32.5 

37.5 

73.  7 

2  4.2 

4  0.4 

32.5 

47.4 

40.4 

36.  8 

71.9 

32.5 

24.2 

32.  5 

58.7 

39.1 

40.4 

47.4 

57.  1 

58.7 

24.2 

40.4 

24.8 

32.5 

40.4 

89.8 

75.3 

47.4 

2  4.2 

32.5 

58.7 

32.5 

32.5 

39.  1 

40.4 

32 . 5 

24.2 

47 . 4 

34.3 

3  2.5 

4  0.4 

24.2 

-12.1 

75.  2 

4  0.4 

24.8 

58.7 

4  7.4 

32.5 

4  0.4 

2  4.2 

3  9.1 

2  9.4 

32.5 

32.5 

32.5 

40.4 

4  7.4 

32.5 

24.  2 

32  .  5 

32.5 

58.7 

40.4 

24.2 

47. 4 

5  7.1 

58.  7 

24.  8 

40.4 

39.  1 

3  2.5 

32.5 

32.  5 

40.4 

2  4.2 

hi .  4 

40.4 

58.7 

24.2 

3  4.3 

37.5 

73  .  7 

86.4 

32.5 

32.5 

47.4 

40.4 

32.5 

39.1 

24.8 

24 . 2 

4  0.4 

58.  7 

32.5 

32.5 

47.4 

32.5 

40.4 

2  4.2 

3  2.5 

3  2.5 

32.5 

32.  5 

4  0.4 

24.2 

47.4 

32  .  5 

3  9.  1 

5  8.7 

4  0.4 

40.4 

24.2 

24.8 

57.1 

58. 7 

47.4 

32.5 

36.8 

71.9 

32.5 

4  0.4 

2  4.2 

32  .  5 

47.4 

58.7 

39.  1 

34.3 

4  0.4 

32.5 

24.2 

32  .  5 

32.5 

41  .  2 

40.4 

47.4 

24.8 

32.5 

24.2 

3  2.5 

40.4 

58.7 

40.4 

47.4 

39.  1 

29.4 

24.2 

32  .  5 

40.4 

40.4 

4  0.4 

24.2 

2  4.2 

47.4 

58. 7 

39.  1 

24.8 

57.  1 

5  8.7 

37.  5 

73.  7 

-12.1 

75  .  2 

39.4 

32.5 

47.4 

32.  5 

3  2.5 

4  0.4 

2  4.2 

32.5 

3  2.5 

34.3 

58.  7 

4  0.4 

47.4 

32 . 5 

24.2 

32.5 

32.5 

4  0.4 

2  4.8 

39.1 
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32  .  5 

24.2 

4  7.4 

3  2.5 

4  . 4 

5'  .  7 

7  2.9 

.1  .  5 

3  2.5 

4  0 . 4 

?4.  2 

32.5 

4  7.4 

4  0 . 4 

3  0.8 

7  5.7 

3  2.9 

3  .  1 

2  4.2 

4  0.4 

5 8 .  7 

2  4.3 

4  7.4 

3°  .  5 

3  2.5 

4  0.4 

2  4.2 

5  7.1 

58.7 

3  2  .  . 

4  *  *  *  4 

4  7.4 

°  4.  2 

j  4  •  j 

’*  >  •  4 

3  .  1 

9  "  .  7 

3  2.:> 

4  "  .  4 

24.2 

47.4 

3  2.5 

4  o .  4 

24.8 

36.8 

7  1.9 

32.5 

32.5 

24.2 

4u .  4 

4  7.4 

58.7 

39.1 

2  9.4 

32.5 

4  •.) .  4 

24.2 

3  2.5- 

37.5 

73.7 

47.4 

4  0.4 

3  2.5 

32  .  5 

24.2 

24.8 

32.5 

40.4 

58.7 

32  .  5 

4  7.4 

2  4.2 

39.1 

4  0.4 

5  7.1 

5  8.7 

3  2.5 

34.  3 

40.4 

24.2 

47.4 

32.5 

32.5 

32.5 

58.7 

4  0.4 

3  2.5 

24.2 

2  4.8 

47.4 

3  9.1 

40.4 

32  .  5 

4  0.  4 

24.2 

47.4 

58.  7 

-12.1 

75.2 

41.2 

100.0 

32.5 

40.4 

24.2 

32.5 

3  2.5 

3  2.5 

40.4 

47.4 

39.  1 

24.2 

24.8 

32.5 

32.  5 

4  0.4 

58.7 

32.5 

32.5 

24.2 

47.4 

40.4 

32  .  5 

34.3 

57.1 

58.7 

40.4 

2  4.2 

39.1 

32.5 

47.4 

32.5 

58.  7 

4  0.4 

32.5 

24.  8 

24.  2 

32.5 

4  0.4 

4  7.4 

37.5 

73.  7 

32.5 

32.5 

2  4.2 

40.4 

39.  1 

29.4 

58.7 

40.4 

47.4 

32  .  5 

24.2 

32.5 

36.  8 

7  1.9 

40.  4 

2  4.8 

32.5 

24.2 

47.4 

32  .  5 

40 . 4 

24.2 

58.7 

39.  1 

34.3 

86.4 

40 . 4 

32.5 

47.4 

32.5 

57.1 

58.7 

40.4 

24.2 

32.5 

4  0.4 

4  7.4 

56.  7 

24.8 

32.5 

32.5 

24. 2 

39.1 

40 . 4 

89.8 

75.3 

32.5 

4  7.4 

32.5 

24.2 

4  0.4 

32.5 

32  .  5 

32.5 

58.7 

40.4 

24 . 2 

47.4 

32.  5 

40.4 

39.  1 

24.8 

32.  5 

4  0.4 

4  0.4 

24.2 

4  0.4 

47.4 

24.2 

58.  7 

34.3 

47.4 

39.1 

57.1 

58.  7 

37.5 

73  .  7 

-12.1 

75.2 

39.4 

24.  2 

32.5 

32.5 

24.8 

40.4 

58.  7 

32.5 

47.4 

24.2 

32.5 

32.5 

40.4 

32.5 

40.4 

2  4.2 

47.4 

39.  1 

29.4 

32.5 

32.5 

32.5 

58.  7 

40.4 

32.5 

24.  2 

32.5 

24.8 

47.4 

40.  4 

41.2 

32.5 

32.5 

24.2 

4  0.4 

47.4 

34.3 

39.  1 

58.  7 

32.5 

4  0.4 

24.2 

32.5 

36.  8 

71.9 

4  0.4 

47.4 

57.1 

58.  7 

24.2 

24.8 

32.5 

32.5 

40.4 

58.7 

39.1 

24.2 

47.4 

32.  5 

40. 4 

32.5 

32.  5 

32.  5 

32.5 

2  4.2 

40.4 

32.5 

4  7.4 

32.5 

40.4 

58.  7 

40.4 

24.2 

24.8 

39.  1 

24.2 

47.4 

34.3 

37.5 

73.  7 

86.4 

32.5 

4  0.4 

32.5 

32 . 5 

58.7 

40.4 

4  7.4 

24.2 

32.  5 

32.5 

40.  4 

2  4.2 

39.1 

24.8 

5  7.1 

58.  7 

47.4 

32.5 

32.  5 

40.4 

58 . 7 

32.5 

32.5 

24.2 

4  0.4 

4  7.4 

32.5 

32.5 

40. 4 

39.1 

29.4 

24.2 

32.  5 

40.4 

2  4.2 

4  7.4 

58.7 

24.8 

-12.1 

75  .  2 

32.5 

32.5 

4  0.4 

3  2.5 

34  .  3 

4  7.4 

2  4.2 

32.5 

40.  4 

39.  1 

32.  5 

32.5 

58.  7 

40.4 

2  4.2 

47.4 

89.8 

75.  3 

32.  5 

40.4 

24.8 

24.2 

57.1 

58.  7 

32.5 

47.4 

4  0.4 

39.1 

58.  7 

4  0.4 

24.2 

36.9 

7  1.9 

32.5 

32.5 

4  7.4 

40.4 

24.2 

37.5 

73.7 

32.5 

4  0.4 

2  4.8 

24.2 

58  .  7 

4  7.4 

39.1 

32.5 

34.3 

40.4 

32.5 

32  .  5 

32.5 

40.4 

2  4.2 

4  7.4 

3  2.5 

32.5 

40.4 

58.7 

41.2 

24.2 

32.5 

4  7.4 

39.1 

24.8 

32.5 

4  0.4 

32  .  5 

24.2 

32.5 

57.1 

58  .  7 

4  0.4 

4  7.4 

58.  7 

32.5 

24.2 

32  .  5 

40.4 

32.5 

40.4 

39.  1 

2  9.4 

47.4 

2  4.2 

32.5 

32.5 

24  .  8 

40.4 

24.2 

58.  7 

3  4.3 

9  6.4 

4  7.4 

32.5 

40.4 

32.5 

32  .  5 

32.5 

24.2 

39.  1 

4  (J.  4 

3  2.5 

4  7.4 

4  0.4 

40.4 

58.7 

40.4 

24.2 

24.2 

24.8 

47  .  4 

39.1 

57.  1 

58.7 

37.  5 

73.  7 

-12.1 

75.2 

39.4 

32.5 

24 . 2 

58  .  7 

47.4 

32.5 

40.4 

32.5 

32.5 

32.5 

32.5 

40.4 

2  4.2 

24.  P 

47  .  4 

40.4 

4  0.4 

39.1 

34. 3 

58.  7 

24.2 

36.  8 

7  1.9 

47.4 

32.  5 

32.5 

32.5 

32.5 

24.2 

40.4 

32.5 

32.5 

40.4 

47.4 

2  4.2 

58  .  7 

24.8 

39.  1 

32.  5 

32.5 

32.5 

4  0.4 

32.5 

4  7.4 

2  4.2 

32  •  5 

40.4 

57.  1 

58.  7 

32.5 

58. 7 

40.4 

24.2 

47.4 

39.1 

29.4 

32.5 

32.  5 
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40.4 

2  4.8 

2  4.2 

3  2.5 

34.3 

47.4 

40.4 

32.5 

58.7 

24.2 

89.8 

75.3 

4  0.4 

32.5 

32.5 

47.4 

39.1 

24.2 

4  0.4 

40.4 

24.2 

37.5 

73  .  7 

24.8 

5  8.  7 

47.4 

32.5 

32.  5 

32.5 

32.5 

40.4 

24.2 

41.2 

47.4 

39.  1 

32.5 

40.4 

5  7.1 

58.  7 

32.5 

58.7 

24.2 

40.4 

32.5 

47.4 

40.4 

24.8 

24.2 

34.3 

86 . 4 

32.5 

3  2.5 

40.4 

39.1 

47.4 

32.5 

58.7 

24.2 

32.5 

40.4 

32.5 

-12.1 

75.2 

24.2 

47.4 

4  0*4 

32.5 

40.4 

24.2 

2  4.8 

58.  7 

39.1 

36.  8 

71.9 

47.4 

32.5 

32.5 

40.4 

32.5 

32.5 

24.  2 

32.5 

40.4 

4  0.4 

47.4 

57.1 

58.  7 

58.7 

24.2 

39.  1 

29.4 

4  0.4 

24.8 

34.3 

47.4 

32.5 

24.2 

32.  5 

32.5 

40.4 

24.2 

37.5 

73  .  7 

58.  7 

47.4 

32.5 

40.4 

32.5 

39.1 

32.5 

32.5 

40.  4 

24.2 

47.4 

2  4.8 

32.5 

32.5 

32.5 

4  0.4 

58.7 

24.2 

32.5 

32.5 

47.4 

4  0.4 

32.5 

24.2 

39.  1 

32.5 

32.5 

40.4 

24.2 

47.4 

57.  1 

58.  7 

58.  7 

40.4 

2  4.8 

34.3 

32.5 

32.5 

40.4 

24.2 

47.4 

40.4 

39.  1 

40.4 

40.4 

40.4 

40.4 

24.2 

58.  7 

4  0.4 

4  7.4 

4  0.4 

2  4.2 

24.2 

24.2 

24.2 

47.4 

24.8 

47.4 

47.4 

58.7 

39.  1 

58.  7 

39.  1 

24.8 

34.  3 

57.1 

58.7 

37.5 

73.  7 

29.4 

36.8 

71.9 

-12.  1 

75.2 

86. 4 

41.2 

89.8 

75.3 

39.4 

95.  1 

I 

I 
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OUTPUT  FOR  SAMPLE  INPUT 


SURFACE  FLAWED  PANEL  ANALYSIS 
DA/DT  IS  DETERMINED  FROM 

DA/Dtt  DATA  -  3COM  STEEL  -  10  CPS  -  DRY  AIR 
AND 

DA/DN  Data  -  3COM  STEEL  -  0.1  CPS  -  ALTERNATE  IMMERSION  IN  SEA  WATER 
F-4J  1 1 A  IN  GEAR  AXLE-PISTON  FORK  STRESS  HISTORY  -  FAST 
100000  BLOCKS  AT  .18500E+03  PSI  DESIGN  LIMIT  STRESS 


BEGIN  SPECTRUM  CRACK  GROWTH  ANALYSIS,  A  =  .05725 


FLICHTS 

A 

REFK 

500. 

.05801 

53 . 89398 

1000  . 

.05883 

54.42016 

1500. 

.05969 

54 . 97040 

2000. 

.06052 

55.49924 

2500. 

.06141 

56.07418 

3000. 

.06225 

56.54943 

3500. 

.0631  5 

57.06480 

4  000. 

.06416 

57.64067 

4  5  00  . 

.  065  1  3 

58. 19100 

5000. 

.06607 

58.72416 

5500. 

.06701 

59.251 24 

6000. 

. 06803 

59.82860 

6500  . 

.06902 

60.38621 

7000. 

.07008 

60.98072 

7500. 

.07116 

61 .56888 

8000  . 

.07301 

62.54464 

85  00. 

.07409 

63  .  1  1  053 

9000. 

.07532 

63. 75800 

9500. 

.07670 

64.47853 

10000. 

.07797 

65.14455 

10500. 

. 0  7  Q6  5 

66.01905 

1 1000. 

.08098 

66 . 70664 

******************************************************************** 

KMAX  EXCEEDS  KSUBC.  PROBLEM  TERMINATED 

********************************  ************************************ 

LAST  CALCULATED  VALUES  ARF. 

FLIGHT  IN  SPEC  TR  U.'f  11171.00 

CRACK  LENGTH  . 8 1 4 1 2 1 8  1  E -0  1 

REFK  . 66932 839E+02 
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APPENDIX  B 

STRESS  INTENSITY  SOLUTIONS  FOR  ELLIPTIC  SURFACE  FLAWS 


1.  SUMMARY  OF  RESULTS  -  Stress  intensity  solutions  computed  from 
a  slice  synthesis  model  and  NASTRAN  three  dimensional  finite 
element  results  were  compared  with  those  obtained  from  three 
dimensional  finite  element  analyses  of  Raju  and  Newman,  References 
16  and  17  to  substantiate  the  model.  Results  of  these  comparisons 
are  shown  in  Figures  B-l  and  B-2.  Good  agreement  is  obtained  for 
the  stress  intensity  at  the  flaw  depth.  Figure  B-l,  and  reasonable 
agreement  is  obtained  for  the  stress  intensity  at  the  surface  for 
all  but  deep  flaws,  Figure  B-2.  The  results  reported  in  Reference 
17b  are  based  on  the  most  detailed  model  employed  thus  far  in  the 
analysis  of  the  surface  flaw. 

Because  good  agreement  is  shown  between  the  results  of  Raju 
and  Newman,  and  the  slice  synthesis  model,  as  shown  in  Figures  B-l 
and  B-2,  the  slice  synthesis  model  results  are  expected  to  be 
reasonably  accurate.  The  slice  synthesis  model  has  the  advantage 
that  finite  width  as  well  as  thickness  can  be  included  in  the 
analysis . 

In  order  to  readily  use  these  results  in  crack  growth  predic¬ 
tion  computer  routines,  equations  relating  stress  intensity  to 
crack  shape,  and  depth  to  plate  thickness  ratio,  have  been 
empirically  established  based  on  the  slice  synthesis  model 
results: 


"a  ■  6a  <”'”i 

kb  =  6B 


where: 


Ka  =  stress  intensity  at  the  depth 
Kg  =  stress  intensity  at  the  surface 
a  =  crack  depth 
c  =  half  surface  length 
t  =  thickness 
w  =  width 


eA  =  Ml^Ml  +  M2a)/>/Q 
eB  =  M1B*(1  +  M2B)/'/Q 


a  =  applied  stress 
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Figure  B-1.  Surface  Flaw  Stress  Intensity  Solution  of  Depth 
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and:  for  (a/c  <_  1.0) 

Q  =  1.0  +  1.47  (a/c)1*64 
M1a  =  1.14  -  0.1  (a/c) 

M1b  =  (a/c)  (1.36  -  0.17 2 (a/c) ) 

< 

5 

M2a  =  (a/t)2  [-0.166  +  0 . 258 (a/t ) 2  +  0.358(c/a)  -  0.242(c/a) (a/t)2] 


M2b  =  (c/t)  [0.0935  +  0.162  (a/t)  -0.0097  (c/a)] 


The  form  of  these  equations  is  the  same  as  presented  in  j 

Reference  20. 

2.  SURFACE  FLAW  GROWTH  ANALYSIS  ROUTINE  j 

Input  Definitions  -  The  computer  program  has  the  following 
inputs: 

a.  title  describing  material  and  product  form  for  reference 
(TITLE) 

b.  number  of  points  in  da/dN  table  ( NDADN ) 

c.  stress  intensity  factor  range  and  crack  growth  rate  for 
each  point  in  da/dN  table  (DK(I),  DADNA  (I)). 

d.  initial  surface  flaw  depth  and  length,  material  thickness, 
width  and  maximum  stress  in  spectrum  of  interest  (A,  C,  T, 

W,  SMAX,  KSUBC) 

e.  number  of  cycles  analyzed  between  printout  of  stress 
intensity  factor  and  crack  growth  (NPRNT) . 

Description  of  Output  -  The  program  outputs  include  title,  A, 

C,  T,  W,  SMAX,  and  a  summary  of  a,  6 a*  kA»  and  c'  3 B '  KB  printed 
at  the  end  of  each  specified  print  interval.  The  print  interval 
is  reduced  by  a  factor  of  ten  where  a  exceeds  0.8t  and  again  when 
transition  to  a  through  crack  occurs.  The  stress  intensity 
correction  factors,  and  SB  are  defined  as 


(a/c  >1.0) 

1.0  +  1.47  (c/a) 1-64 
1.04  (c/a) 

1.0  +  0.188  (c/a) 
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LISTING  OF  PROGRAM 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  PROGRAM  COMPUTES  CRACK  GROWTH  AT  SURFACE  AND  THROUGH  DEPTH 
FOR  A  SEMIELLIPTICAL  SURFACE  FLAW.  STRESS  INTENSITY  FACTORS 
ARE  BASED  ON  MCA  I R  SLICE  SYNTHESIS  MODEL  RESULTS.  CRACK 
FRONT  PENETRATION  OF  SACK  FACE  IS  PREDICTED  BY  LINEAR  EXTRAPOLATION 
OF  STRESS  INTENSITY  FACTOR  RESULTS.  ASSUMED  LOADING  IS  CONSTANT 
AMPLITUDE,  R-0.  ,  SINUSIODAL  LOADPNG.  PROCRAM  WAS  DEVELOPED  BY 
C.R.SAFF  OF  MCDONNELL  AIRCRAFT  CO.,  ST.  LOUIS,  MISSOURI, 

TELEPHONE  (314)  232-3356. 

INTEGER  T ITLE (20) 

REAL  DADN  A ( 2  5 ) ,DK(25)  ,KSUBC 
PI-3. 14159265 

INPUTS  ARE  -  TITLE 

FTY.FCY 

NDADN , FREQ 

DK(2  5 ) , DADN  A ( 2  5  ) 

A, C,T,W, SMAX, KSUBC 
NPRNT 

MATERIAL  PARAMETERS 

TITLE  DESCRIBES  MATERIAL  &  PRODUCT  FORM  FOR  REFERENCE 

FT Y-MO NOTONIC  TENSILE  YIELD  STRENGTH 

FCY-CYCLIC  TENSILE  YIELD  STRENGTH 

NDA Dt«  -NO .  OF  POINTS  IN  DA/DN  TABLE 

FREQ-LOAD  CYCLE  FRFOUENCY 

DK(I)-DELTA  K  VALUES  OF  POINTS  IN  DA/DN  TABLE 
DADN A (I  )-PA/DN  VALUES  OF  POINTS  IN  DA/DN  TARLE 

GEOMETRIC  PARAMETERS 

A-INITIAL  SURFACE  FLAW  DEPTH 
C-INITIAL  SURFACE  FLAW  LENCTH 

T -MATER IAL  THICKNESS  -  CRACK  DEPTH  DIRECTION 
W-MATERIAL  WIDTH  -  CRACK  SURFACE  DIRECTION 
SMAX-MAXIMUM  SPECTRUM  STRESS  APPL'F.D 

KS  UBC-CRIT ICAL  STRESS  INTENSITY  FACTOR  FOR  MATERIAL  THICKNESS 

NPRNT-NO.  OF  CYCLES  BETWEEN  PRINTOUT  OT  K  DATA  AND  CROWTH 

READ ( 5 , 1000)  TITLE 
V'RITE(6 , 1010)  TITLE 
READ(5, 1020)  FTY.FCY 
REAO(5, 1030)  NDADN, FREQ 
DO  10  I- l, NDADN 
READ( 5,1020  )  D  K ( I ) , 0 ADNA ( I ) 

DK ( I )-AL0G10(DK(I) ) 

10  DADNA( I )-ALOGlO(PADNA(I ) ) 

OUTPUT .'INPUT  A, C.T.W, SMAX, KSUBC' 

REA 0(5,  1020)  A, C.T.W, SMAX, KSUBC 

OUTPUT, 'INPUT  NO.  OF  CYCLES  BETWEEN  OUTPUT' 

REAl)(5,  1030)  NPRNT 

WRITE(6, 1040)  A, C.T.W, SMAX 

ICIIK-O 

IBRK-0 

CYC-O. 

NP-0 

WR I TE ( 6 , 1050) 


nonon  oooo 
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DO  LOOP  80  COMPUTES  STRESS  INTENSITY  FACTORS  AND  GROWTH  AT 
DEPTH  AND  SURFACE  OF  FLAW 

DO  80  1-1,100000 

IF(I.EQ.l)  N INC-1 

IF  ( I.  EQ.  2)  N  INC-tlPRNT/  1000-V 

IF( I • GT . 2 )  NINC-NPRNT/1000 

IF ( N INC • LT .  1  )  N INC - 1 

CYC -C YC+N INC 

ICYC-CYC 

IF ( ICHK. EQ. 1 )  GO  TO  AO 

AOC-A/C 

COA-l . / AOC 

AOT-A/T 

IF ( AOC •  GT  •  1 .  )  GO  TO  20 
0-l.+1.47*<A0C**1.64> 

XM 1 A- 1 • 1 3-0 .09*AOC 
XM1B-U  .  36-0. 185* AOC)* AOC 
GO  TO  25 

20  Q-l .+1 . 47*COA**l . 64 
XM 1A-1 • 04*COA 
XM  IB  - 1  *  +0  •  1  88  *COA 

2  5  XM 2A- < -0 . 16  6+0. 2 58 *AOT*AOT+0 . 358 *C0A-0 . 2 4 2 *COA* AOT *AOT ) * ( AOT *AOT ) 
XM2B-(0. 09  3  5+0. 1 6 2* AOT-0 . 009 7*COA) * ( COA* AOT ) 

BETAA-XM 1 A* ( 1 . +XM 2A ) /SORT (0 ) 

BETAB-XM1B* ( 1 . +XM 2B ) / SQRT (O ) 

STRESS  INTENSITY  FACTOR  EXPRESSIONS  ARE  ASSUMED  GOOD  UNTIL 
FLAW  DEPTH  EXCEEDS  90 Z  OF  THICKNESS.  VALUES  AT  A/T-0.8  AND 
A/T-0.9  ARE  USED  FOR  EXTRAPOLATIONS  BEYOND  A/T-O. 9. 

IF ( A . LT . 0 . 8*T )  GO  TO  50 
N PRNT -N PENT /  I  0 
IF(NPRNT.LT.l)  NPRNT-1 
ICHK-1 
A I -0 . 8  *T 
A2-0 • 9*T 

IF ( AOC . GT . I . )  CO  TO  30 
Q=l .+1 . 47* (AOC**l . 64) 

XM1A-1. 1 3—0 .09  *AOC 

XII  IB-  (1  •  36-0 .  1 8  5  *  AOC  )  *  AOC 

GO  TO  35 

30  Q-l .+1 . 47*COA**l . 64 
XH1A-1. 04*COA 
XM IB-1 . +0. I 38*COA 

35  XM2A 1  —  (— .  1 66+. 2 58*A I /T*A 1 /T+. 358 *COA  — .  2 42 *COA* A 1/T*A1/T)*(A1/T)**2 
XM2A2-(-. 166+. 258*A2/T*A2/T+. 358*COA-. 2 42* COA* A2 / T* A2 / T ) * ( A 2 /T ) * *2 
XM2B I -(0.093 5+0.  1 6  2  * A  1 /T-0 . 009  7*C0A ) * ( COA* A 1 /T ) 

XM  2B2- (0 . 0935+0 . 1 62* A2 /T-0 . 009  7 *COA ) * ( COA* A2 /T) 

BAI-XM1A*( I.+XM2A1)/S0RT(Q) 

BA2-XM1A* (1. +XM2A2) /SQRT(q) 

BB 1 -XM IB  * ( 1 . +X  1 2B 1) /SORT(q) 

BB2-XM18*(l. +XM2B  2)/SQRT(Q) 

40  BETAB-(A-A1 ) /( A2-A1 )*( BB2-BB1 )  +  BB  1 
BETA  A- ( A -A  1 ) / ( A  2  —  A  1 ) *( BA2-BA1 )+RA 1 
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50  AK-RHTAA*SMAX*SQRT(PI *A) 

BK-BETA8*SMAX*SQRT (PI *C/C0S (PI*C/W) ) 

IF (IRAK. ME. 1)  GO  TO  60 

WRITE (6, 1070)  CYC,A,BETAA,AK,C,BETAB,BK 
WR IT  E ( 6 ,  1100) 

IBRK-I BRK+1 
CO  TO  70 

60  I F( IC YC. LT. NP)  GO  TO  70 
NP-NP+NPRNT 

WRITE (6 , 1070 )  CYC, A , BETAA , AK , C , B  ETAB , BK 

FRACTURE  IS  ASSUMED  TO  OCCUR  WHENEVER  THE  STRESS  INTENSITY 
FACTOR  AT  EITHER  THE  SURFACE  OR  THE  DEPTH  EXCEEDS  KSUBC. 

70  1F(AK.CT .KSUBC)  GO  TO  120 
IF(8K.CT. KSUBC)  GO  TO  130 

DA/ON  LOOK  UP  IS  LOC (DA/DN )  VS.  LOG (DELTA  K).  CRACK  GROWTH  RATE 
CURVE,  is  THE  DEPTH  DIRECTION  IS  ASSUMED  SAME  AS  ALONG  SURFACE. 

C 

1F( BK.GT . 0. S*KS1IBC. OR . AK.CT.O. 8*KSUBC )  N INC-1 
AKI-ALOCIO(AK) 

BKI-ALOCIO(BK) 

CALL  TLU(DK, DADNA, NDADN ,AK1, ADADN ) 

CALL  TLU(DK, DADNA, NDADN , RK1 , BDADN) 

ADADN” 1 0 . ** ADADN 
BDADN” 1 0 •**  BDADN 
A”A+ADADN*NINC 
C-C+BDADN  *N INC 
C 

C  CRACK  IS  ASSUMED  TO  CROW  AT  A  FIXED  A/C  PATIO  AFTER  CRACK 

0  FRONT  BREAKTHROUGH  UNTIL  BETAR  EXCEEDS  1.  THEN  FLAW  IS  TREATED 

C  AS  A  CENTER  CRACK  Of  LF.NCTH  2C. 

C 

IF ( BETAB . CE  l .  )  CO  TO  90 
IF(A.LT.T)  CO  TO  80 
IF(IBRK.CE.l)  CO  TO  80 
URITE(6, 1060) 

IBRK-I BRK+1 
80  CONTINUE 
90  UKITF,(6,  1080  ) 

WRIT  E ( 6 ,  1070)  CYC,A,BETAA,AK,C,BETAB,BK 

A-C 

WRITE (6 ,  1090  ) 

URITE(6, 1070  )  CYC, A 
DO  110  1-1,100000 
CYC-CYC+N INC 
IC YC-CYC 

CK-SMAX*SQRT (PI* A/COS (PI*A/W) ) 

IF(ICYC.LT.NP)  CO  TO  100 
WRITE(6, 1070)  CYC, A 
NP-NP+NPRNT-1 
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100  IF(CK.GT.KSUBC)  GO  TO  140 
IF(CK. CT .0. 9*KSUBC )  NINC-1 
CK-ALOGIO(CK) 

CALL  TLU!DK, DADNA, NDADN , CK.DADN ) 

DAi)N-10.**DADN 
A-A+0ADN*N INC 
110  CONTINUE 
STOP 

120  WRIT  E ! 6 , 11  10) 

WRITF.  (6  ,  1070  )  CYC,A,BETAA,AK,C,BETAB,BK 
STOP 

130  URITE!6,  1  120) 

WRIT  F.(6  ,  1070  )  CYC  ,  A,  BETA  A,  AK,  C,  BE  TAB  ,  BE 
STOP 

140  WRIT E ( 6 , 1130) 

WRITE (6 , 1070)  CYC, A 
STOP 

1000  F ORMAT ( 2  OA  4 ) 

1010  FORMAT (1H1./20A4) 

1020  FORMAT ( 1  OF  10.0) 

1030  FORMAT (I  10) 

1040  FORMAT!/'  AO  -  ',FS.3,'  CO  -  '.F6.3, 

4/'  T  -  ' , F8 • 3, '  W  -  ',F6.3 ,/'  SMAX  -  ',F8.3) 

10  50  FORMAT (/T 4, ' C YC ' , T 1 6 , ' A ' ,T24,  'BETA A' , T34 , ' K A ' , T4 4 ,  '  C ' , 
4'BETAB ' ,T42 , ' KB' /) 

1060  FORMAT!'  CRACK  FRONT  PENETRATES  BACK  FACE  AT') 

1070  FORMAT! F8.0, 2F10.4,F8.2,2F10.4,F8.2) 

1080  FORMAT!'  SURFACE  CRACK  BECOMES  THROUGH  CRACK  AT') 

1090  FORMAT!//'  EQUIVALENT  THROUGH  CRACK  LENGTH  IS'/, 

&  T4  , '  C  YC  ' ,T16,'A'/) 

1100  FORMAT ( 1 H  ) 

1110  FORMAT!/'  FRACTURE  AT  A') 

1120  FORMAT!/'  FRACTURE  AT  C') 

1130  FORMAT!/'  THROUCU  CRACK  FRACTURE') 

F.NU 


SUBROUTINE  TLU ( X , Y , N , X VAL, Y VA L ) 

DIMENSION  X!N)  ,Y!N) 

IF! X !N ) .OE . XV AL )  GO  TO  10 
I-N 

GO  TO  30 
10  00  20  1-1 , N 

IF ( X  < I ) . GE.XVAL)  GO  TO  3  0 
20  CONTINUE 
30  IF! I . FQ. 1)  1-2 

YVAL-(Y !I ) -Y ! I - 1 ) ) / <  X ! I ) -X ! I -1 ))*(XVAI,-X(I-1 ) )+T(I-l> 

RETURN 

END 
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SAMPLE  INPUT 


DA/ON  DATA  -  300M  STEEL  -  10  CPS  -  DRY  AIR 


248. 

228  . 

13 

10. 

12. 

.035-6 

14. 

.  90E-6 

16. 

1. 31E-6 

18. 

1 . 85  E-6 

20. 

2. 55E-6 

22.5 

3. 6E-6 

25. 

4. 75E-6 

30. 

7. 5E-6 

35. 

1 2. 8E-6 

40. 

18.2E-6 

45. 

25. 5E-6 

50. 

40. OE-6 

60. 

600. E-6 

0.04 

0.0545 

500 
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OUTPUT  FOR 

SAMPLE 

INPUT 

PA/DN  DATA 

-  3 00M 

STEEL  -  10 

CPS  -  DRY 

AIR 

AO  - 

.  040 

CO  - 

055 

T  = 

.250 

W  -  1  . 

000 

SMAX  - 

100.000 

CYC 

A 

BETAA 

KA 

C 

BETAB 

KB 

1. 

.0400 

.  7«1  3 

27.69 

.0545 

.6696 

2  7.91 

500. 

.0432 

.  77  58 

2  8.57 

.0578 

.6769 

29.08 

1000. 

.0466 

.7712 

29.52 

.0614 

.6836 

30. 32 

1500. 

.0504 

.  76  81 

30.5  7 

.0656 

.6893 

31.64 

2000. 

.054  7 

.  7661 

31.76 

.  0705 

.6945 

33.10 

2500. 

.0596 

.  7651 

33.11 

.0762 

.6996 

34.74 

3000. 

.0653 

.  7648 

34.65 

.0830 

.  7050 

36.61 

3500. 

.0720 

.  7651 

36.38 

.0907 

.7111 

38.  75 

4000. 

.0795 

.  76  69 

38.34 

.0998 

.7172 

41.18 

4500. 

.0883 

.7712 

40.62 

.1107 

.  72  33 

43.99 

5000. 

.0988 

.  77  93 

43.41 

.1243 

.  72  92 

47.39 

5500. 

.  1120 

.  8025 

4  7.60 

.  1457 

.  72  95 

52.10 

FRACTURE 

AT 

C 

5  767. 

.1284 

.  882  5 

56.05 

.  1  934 

.  6986 

60.11 
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